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ANALYSIS  AND  OPTIMIZATION  OF  THE  MULTIEFFECT 
MULTISTAGE  FLASH  DISTILLATION  PROCESS 


CHAPTER  1 
INTRODUCTION 

The  present  study  is  directed  to  the  System  analysis  and  optimiza- 
tion of  a  multieffect  multistage  (MEMS)  flash  distillation  process. 
The  MEMS  flash  distillation  process  is  a  rather  recent  development  in 
the  flash  distillation  technology  and  offers  the  most  promise  in  the 
foreseeable  future  for  producing  large  quantities  of  potable  water 
economically  from  seawater  (1). 

A  better  understanding  of  the  MEMS  system  is  obtained  by  following 
the  developments  of  the  process.   Regular  distillation  is  a  familiar 
water  purification  process.   Flash  distillation  was  introduced  because 
of  better  control  of  scale  formation  (2).   In  the  flash  distillation 
process,  heated  saline  water  is  released  into  a  closed  vessel  which  is 
maintained  at  a  lower  pressure  than  the  vapor  pressure  of  the  solution. 
Since  the  vapor  simply  flashes  off  the  warm  liquid,  the  resulting  pre- 
cipitates form  in  the  liquid  and  not  on  the  heat  transfer  surface  (3). 

The  brine  concentration  in  a  flashing  chamber  is  nearly  uniform 
due  to  the  vigorous  mixing  resulting  from  the  flashing  and  is,  there- 
fore, equal  to  that  of  the  discharge  stream.   In  a  single  stage  oper- 
ation, the  feed  brine  with  low  concentration  is  mixed  with  flashing 
brine  with  high  concentration.  This  causes  a  large  amount  of  free 
energy  loss  due  to  the  irreversible  mixing  of  two  solutions  with  con- 
siderable concentration  difference.   However,  the  concentration  differ- 
ence between  stages  in  a  multistage  operation  is  considerably  reduced. 
Therefore,  the  thermodynamic  efficiency  for  this  operation  is  greatly 
improved.   This  multistage  operation  is  the  so-called  "single-effect 
multistage  (SEMS)"  flash  distillation  process  (4). 
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A  recirculated  SEMS  is  a  single-effect  multistage  flash  distilla- 
tion with  a  recycle  operation.   The  main  reasons  for  using  a  recycle 
stream  are  to  increase  the  total  heat  capacity  of  the  flashing  brine 
and  the  percentage  conversion  of  the  brine  feed  into  fresh  water  (S,  6). 
Due  to  the  high  latent  heat  of  vaporization  of  water  and  the  low  heat 
capacity  of  the  aqueous  solution,  the  solution  cools  off  considerably 
when  only  a  small  fraction  of  the  solution  is  flash  evaporated.   In 
the  flash  distillation  process,  the  highest  flashing  temperature  is 
limited  by  the  scale  formation  problem  and  the  lowest  temperature  is 
limited  by  the  temperature  of  the  seawater  which  is  used  as  the  coolant. 
The  percentage  conversion  of  a  flash  distillation  process  without  re- 
cycle, which  is  operated  within  the  temperature  range  mentioned  above, 
is  less  than  20%  (7,  8).   Since  the  saline  water  feed  stream  has  to 
be  pumped  and  pretreated,  a  low  percentage  conversion  of  feed  water 
into  fresh  water  will  result  in  a  poor  overall  economy  for  the  process. 

A  MEMS  process  consists  of  several  SEMS  plants  with  recycles 
connected  in  series.   Each  SEMS  system  is  considered  as  an  effect. 
W.  R.  Williamson  et.  al.  (7)  have  summarized  the  advantages  of  the  MEMS 
system  as  compared  to  the  SEMS  system  as  follows: 

(a).   Reduces  total  feed  treatment  cost  by  50%. 

(b).   Reduces  heat  transfer  surface  by  at  least  20%. 

(c).   Allows  for  more  stages  at  the  hot  end  of  the  plant  and 

fewer  stages  at  the  cold  end  of  the  plant. 
Thus,  a  MEMS  system  lends  itself  to  better  control  of  the  operating 
variables  so  that  a  lower  water  cost  can  be  achieved. 

The  MEMS  process  is  described  in  detail  in  Chapter  2.   In  Chapter 
3  a  mathematical  model  which  fairly  accurately  describes  the  process  is 
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developed.   Each  effect  is  first  assumed  to  consist  of  an  infinite 
number  of  stages  (infinite  stage  operation)  and  differential  equations 
are  set  up  to  obtain  the  idealized  performance  equations.   Since  an 
actual  plant  consists  of  a  finite  number  of  stages  (finite  stage 
operation) ,  appropriate  correction  terms  are  added  to  the  idealized 
performance  equations.   The  capital  and  operating  cost  equations  are 
set  up  in  Chapter  4. 

The  discrete  maximum  principle  combined  with  a  search  technique 
is  used  to  optimize  the  MEMS  system.   Two  search  techniques  are  used 
in  this  study:   the  parametric  search  and  the  Simplex  method.   The 
optimization  procedure  and  numerical  results  are  illustrated  and  the 
comparison  of  the  two  methods  is  discussed  in  Chapter  5.   The  final 
optimal  policy  of  the  system  and  the  capital  and  operating  costs 
allocation  are  given  in  Chapter  6.   The  computer  program  of  each  method 
and  the  sample  results  are  listed  in  the  Appendix. 
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CHAPTER  2 
PROCESS  DESCRIPTION 

Figure  1  presents  a  simplified  process  flow  diagram  of  a  three- 
effect  multistage  flash  system  and  Figure  2  depicts  a  typical  effect, 
the  n-th  effect,  of  the  system.   The  critical  locations  of  the  system 
are  denoted  by  numbers,  n',  n",  and  n,  which  divide  the  system  into 
various  sections,  namely  sections  H-n,  MR-n,  HR-n,  and  R-n.   The  n-th 
effect  consists  of  preheater  H-n,  mixing  section  MR-n,  heat  recovery 
section  HR-n,  and  heat  rejection  section  R-n.   From  Figure  2  it  is 
easily  seen  that  the  preheater  of  the  n-th  effect,  H-n,  coincides  with 
the  heat  rejection  section,  R-(n-l),  of  the  (n-l)-th  effect. 

F,  L,  and  Rn  represent  respectively  the  flow  rate  of  the  feed 
brine,  flashing  brine,  and  recycle  brine  in  the  n-th  effect.   Wn 
represents  the  condensate  produced  in  the  n-th  effect.  The  feed  brine 
and  the  recycle  brine  together  are  referred  to  as  the  non-flashing 
brine  stream.   Tf,  Tj ,  and  Tc  represent  respectively  the  temperature 
of  the  flashing  brine,  non-flashing  brine,  and  condensate.   The  sub- 
script is  used  to  indicate  the  location.   For  example,  (T-),  and 
(Tj)nt  are  respectively  the  temperature  of  the  flashing  brine  at  lo- 
cation 1  and  temperature  of  the  non-flashing  brine  at  location  a'. 

In  Figure  1,  the  seawater  feed  is  heated  in  section  R-3  and  then 
degasified  to  remove.  C02  and  other  dissolved  gases.   After  being 
heated  successively  in  sections  HR-3,  H-3,  HR-2,  H-2,  and  HR-1 ,  it  is 
mixed  with  the  recycle  brine  R],  to  form  a  brine  stream  which  is  heated 
in  the  brine  heater,  H-l ,  and  then  introduced  into  the  first  effect  as 
the  flashing  brine  [L)it, 
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As  is  shown  in  Figure  2,  the  flashing  brine  at 

location 

a  is 

divided  ir.to  two  streams:   one  stream,  (l)  ,  is  fed 

into  the 

(n+l-)-th 

effect  and  the  other  stream,  R  ,  is  recirculated  by 

the  recyc 

:le  pump, 

■V 

heated  in  sections  HE-n  and  H-n,  and  then  mixed 

with,  the 

brine 
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A  stage  within  each  effect  consists  of  a  flashing  chamber  arc 
a  condensing  chamber  and  a  demister  which  separates  the  two  chambers. 
Each  stage  is  maintained  at  a  lower  pressure  than  the  preceding  one. 
Brine  flows  from  stage  to  stage,  giving  up  _additional  vapor  as  the 
pressure  drops;  the  vapor  then  passes  through  the  demister  to  the 
condensing  chamber,  where  it  is,  condensed  to  heat  the  nor.-f lashing 
brine . 
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CHAPTER  3 
PROCESS  ANALYSIS 

3-1.   Outline  of  Process  Analysis 

Quantitative  relations  among  the  operating  variables  are  derived 
in  the  following  sections.   The  performance  of  a  MEMS'system  is 
characterized  by  the  temperature  -  composition  diagram  for  the  (n-l)th 
and  the  n-th  effects  in  Figure  3.   The  general  approach  is  to  obtain 
idealized  performance  equations  by  assuming  infinite  stage  operation 
in  each  effect  and  applying  correction  terms  for  the  finiteness  of  the 
number  of  stages. 

The  lines,  a-b-c,  and  d-e-f,  show  how  the  temperature  of  the 
flashing  brine,  T^,  decreases  as  the  concentration  of  the  brine,  Cf, 
increases  in  an  infinite  stage  operation  in  the  (n-l)-th  effect  and 
n-th  effect,  respectively.   The  relations  representing  these  lines 
are  derived  in  section  3-5.   The  concentration  gaps,  between  (n-1) 
and  n',  n  and  (n+1 ) ,  are  caused  by  mixing  of  brines  due  to  recirculation 
in  the  n-th  and  (n+l)-th  effects,  respectively.   The  stepped  lines 
along  the  lines,  a-b-c  and  d-e-f,  represent  the  temperature  of  the 
flashing  brine  in  the  (n-l)-th,  and  n-th  effects  respectively  in  an 
actual  process  where  the  number  of  stages  in  each  effect  is  finite. 

The  lines,  a'-b'-c1,  and  d'-e'-f,  show  the  relations  between 
the  condensate  temperature,  T  ,  and  the  flash  brine  composition,  C^, 
at  various  locations  in  the  system  for  an  infinite  stage  operation. 
The  stepped  lines  along  a'-b'-c'  and  d'-e'-f  again  represent  an 
actual  finite  stage  operation. 
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3.     Temperatures    Tf    ,  Tc    ,   Ti      vs      concentration      cf 
flashing       brine      Cf      (  schematic  )  . 
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The  vertical  distance  between  the  two  3Cts  of  linos  described  in 

the  last  two  paragraphs  represents  (T,  -  T  )  at  various  locations  in  the 

system.  This  difference  is  due  to  the  boiling  point  elevation  of  the 

flashing  brine  and  the  pressure  difference  across  the  demister  at  each 

location.  This  difference  will  be  denoted  by   a.  The  magnitude  of  * 

varies  throughout  an  effect;  this,  is  mainly  due  to  the  Varying  composition, 

and  consequently  the  varying  boiling  point  elevation.   The  average 

value  of  a  in  the  n-th  effect  is  denoted  by  a   and  is  derived  in  section  3-6* 

n 

Similarly,  the  lines  a"-b"-c"  and  d"-e"-f"  in  Figure  3  show  the 
relation  between  non-flashing  brine  temperature,  T  . , ■ and  flashing  brine 
composition,  C^ .    AT  is  used  to  represent  the  temperature  difference 
(T.-T.).   As  illustrated  in  the  figure,   AT,  in  an  infinite  stage  operation 
is  nearly  constant  within  a  heat  recovery  section.   However,  it  is  no- 
constant  within  a'  heat  rejection  section.   For  example,    AT  varies  gradually 

from  (  AT)  „  to  (  AT)   in  the  heat  rejection  section  R-n.  These  items 
\    'n"    v    'n 

are  explained  further  in  section  3-8.  The  average  &    in  the  n-th  effect, 

AT  ,  is  derived  in  section  3-9 • 

The  heat  loads  in  the  brine  heater  and  the  n-th  effect  are  derived 
in  section  3-^*-   A t  is  used  to  denote  the  temperature  difference  required 
for  the  heat  transfer  in  the  various  sections.  The  average  values  of 

A t  in  the  brine  heater  and  the  n-th  effect  are  calculated  by  the 
relations  derived  in  section  3-10.   From  these,  the  relations  for  the 
heat  transfer  area  requirements  are  derived  in  section  3-11 •  The  pumping 
head  required  for  each  of  the  circulating  pumps  is  derived  in  section 
3-12. 
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3-2.   Flow  Rates  of  Flashing  Brine,  Kecycle  Brine  and  Condensate  Streams- 

The  notation  representing  the  various  fluid  streams  has  been 
defined  in  Chapter  II;  The  recycle  ratio  in  the  n-th  effect,  r  , 
is  defined  by 

R 


*  'n-1 


=  1,  2,  i,  (J) 


where 

(L)f 


The  flow  rate  of  the  brine  stream  leaving  the  n-th  effect, 
(L)   is  related  to  its  concentration,  (C,)  ,  by  the  following  equation, 

C_  .  . 

(L)„  =   F  —2 ,    a  -  1,  2,  3,  (M 

(Cf>n 

where  C-  is  the  salt  concentration  in  the  feed.  Therefore,  by  combining 
equations  (3)  and  C+),  the  flow  rate  of  the  recycle  brine  stream  can  be 
written  as  .  . 


cv 

R 
n 

- 

r  F     - 
n 

<Cf)n-l 

Note 

that 

n  -  1,  2,  3-  (5) 


(cf)0    -  <y. 

The  flow  rate  of  the  brine  stream  entering  the  n-th  effect  is  - 

given  by 

C-i 

(-)  ,   -  (L)„  ,  +  B    -  F   —  (Ur  ),'  n  =  1,  2,  3-   (6) 

v  'n '         r.-l     n  ,   s    v    r.  k  ' 

^  f "r.-l 
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From  the  overall  material  balance  over  the  n-th  effect,  we  obtain 
the  following  equation  for  the  flow  rate  of  the  condensate  produced 
in  the  n  th  effect. 


Wn   =  F  j   *-F     -    <-F     \.      ,      n  =  1,  2,  3.  (7) 

Jcf>n-1    (Cf)„ 


Therefore,  the  total  water  production,  2wn,  is 
SW   =   ZWn=  F   1  -    CF   1    .  (8) 

a-J      I     (cf)3 


3-3.   Mixing  of  the  Recycle  Brine  Stream  with  the  Flashing  Brine  Stream. 

As  mentioned  previously,  the  brine  streams  having  different  com- 
positions are  mixed  isothermally  at  mixing  points.   By  making  a  salt 
material  balance  at  the  mixing  point,  Mn,  as  shown  in  Figure  2,  we  obtain 

<L>n-l   (cf)n-l   +  «n  <cf>n  =  <L>n''  <Cf>n' 

Substituting  equations  (4),  (5),  and  (6)  into  the  foregoing  equation 
yields 

(Cf)n,  =  (Cf>n-1  +  fnjffj  n      ,  n  =  1,  2,  3.  (9) 

1  +  rn 

3-4.   The  Heat  Loads  in  the  Brine  Heater  and  the  n-th  Effect. 

The  heat  load,  qs ,  in  the  brine  heater,  is  a  very  important  • 
operating  variable.   A  large  value  of  qs  gives  rise  to  an  increased 
steam  cost.   But  it  also  gives  rise  to  a  large  temperature  difference 
A  t  for  heat  transfer  and  consequently  a  low  plant  cost. 
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If  the  heat  load  due  to  flashing  of  the  condensa 

to  stream  : 

.s 

neglected,  the  heat  load  in  the  n-th  effect  can  be  ap 

proxima' 

ied 

by  the  latent  heat  required  for  condensate  production 

in  the 

a-1 

;  h 

effect.   Therefore,  we  have 

-  f  CP        CF     ") 

-  i,  a 

1        S  •. 

(10) 

V°*W    (Cf)n  J   ' 

where   X  is  the  latent  heat  of  flashing  brine. 

J»5"   Temperature  and  Composition  of  the  Flashing  Brine,  T_„  vs. 

cf 

In  Figure  4,  which  represents  a  flashing ■ chamber 

of  an  : 

;  ■*,■&< 

.nite 

stage  system,  L,  C^.,  T,  and  h„  are  respectively  the  quantity. 

I 

concentration,  temperature  and  unit  e.nthalpy  of  the  f' 

Lashing 

b?i 

,ne  • 

Let  dV  be  the  quantity  of  water  vapor  evaporated,  and 

let  :-: 

V 

be 

the 

unit  enthalpy  of  the  vapor.  .A  total  material  balance 

gives 

L  =  L  +  dL  +  dV. 

or 

dL  =   -  dV. 

(11) 

A  salt  balance  gives 

L  Cf  =  (L  +  dl)   (Cf  +  dCf), 

Neglecting  the  tern  dLdC»  in  this  equation  yields 

dL_    m        _  ^f_ 

(12) 

L             Cf 

The  enthalpy  balance  is 

L  h-  .  (L  +  dL)(hf  +  dhf)  +  IIvdV. 
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stage  system,  L,  C^.,  T,  and  h„  are  respectively  the  quantity. 

I 

concentration,  temperature  and  unit  e.nthalpy  of  the  f' 

Lashing 

b?i 

,ne  • 

Let  dV  be  the  quantity  of  water  vapor  evaporated,  and 

let  :-: 

V 

be 

the 

unit  enthalpy  of  the  vapor.  .A  total  material  balance 

gives 

L  =  L  +  dL  +  dV. 

or 

dL  =   -  dV. 

(11) 

A  salt  balance  gives 

L  Cf  =  (L  +  dl)   (Cf  +  dCf), 

Neglecting  the  tern  dLdC»  in  this  equation  yields 

dL_    m        _  ^f_ 

(12) 

L             Cf 

The  enthalpy  balance  is 

L  h-  .  (L  +  dL)(hf  +  dhf)  +  IIvdV. 
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Neglecting  the  terra  dLdh,  and  substituting  equation  (11)  into 

the  above  equation  yields 

L  dh.  =  (H   -  h,)  dL. 
I    v  v     t' 

Since  dh.  «  C   dT  ,  and  H   -  h,.  can  he  approximated  by  the  latent 

heat  of  vaporization,  X   ,    the  above  equation  becoc.es 

C              dL 

x       iTt    - 

i 

By  substituting  equation  (12)  into  the  foregoing  equation,  we  obtain 

■c                        dCf 

-—  dif  .......                              (13J 

Lf 

Assuming  that  C  /x        is  constant  and  integrating  the  above  equation 
P 

between  locations  n'  and  n,  we  have 

(cO      c   (                             ~) 

m        r"            p     )  ^Tf>-    -   (*A   y    ■               ("0 

In   c    /    I  Ti-          I  n   f 

(cf)n,     x   L_           J  . 

Substituting  equation  (9)  into  this  equation  and  noting  that 

^fh>      *   ^A.-l' 

wo  obtain                                                — \ 
In  (C.)   .  In   *  n  X — *-fi   t  — 2-    4  ^f'n-l    l'fV  f 

r+  rB           ■ ' ' »        L                  J 

■ 

n  -  1,  2,  5.             (15) 

Ihia  is  the  equation  of  the  straight  line,  d-e-f,  in  Figure  3. 
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3-6.   Average  <x  in  the  n-th  Effect,  d-n    . 

oL  is  defined  as  the  difference  in  temperature  of  the  flashing 
brine  and  the  condensate  in  a  state.   dn  is'  the  average  value  of  OL 
in  the  n-th  effect.   The  value  of  CL  in  a  stage  depends  on  the  boiling 
point  elevation  of  the  brine  and  the  drop  in  condensation  temperature 
due  to  the  demister  pressure  drop.   Due  to  lack  of  information,  the 
drop  in  condensing  temperature  due  to  the  demister  pressure  drop  was 
assumed  to  be  1°F  in  each  effect. 

Figure  5  shows  the  boiling  point  elevation  of  the  brine  solution 
as  a  function  of  its  composition  at  constant  temperature.   The  average 
value  of  the  boiling  point  elevation  in  each  effect  was  evaluated  at 
the  average  temperature  in  the  corresponding  effect.   The  average 
temperatures  in  the  first,  second  and  third  effects  are  assumed  to  be 
225  F,  175  F  and  125  F,  respectively.   For  further  simplification,  the 
average  value  of  the  boiling  point  elevation  in  each  effect  was  calcu- 
lated from  the  slope  of  each  curve  at  the  appropriate  concentration 
ranges  instead  of  the  curve  itself. 

Therefore,   OL  n  can  be  expressed  as  functions  of  the  average  brine 
composition  in  the  n-th  effect  by  the  following  equations, 
,        (C-).,   +  (C  ) 

oi      =  i.oi  +  — —    — £-± ±-±_ 

1  0.03  2 


=  1.01  +  l 

0.03 


CF+rl<Cf>l  } 

1+r,      +  <Cf>l 
*—  (16) 
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aa  =  1.0075 

, 

(cf)?1  t  (cf)2 

-    0.0547 

■■-■-. _  — ■  +  (.C^K 

.  1.0075 

t 

1  +  r,          r'z 

^                                  '17s 

0.0347 

\l/; 
2 

a3  .  O.32 

i 

Co,),,    ♦    (cf), 

2 

i.0375 

-  0.^2 

+ 

1  T  '5               .   (ie) 

' 

O.0375 

2 

Jt7.  The  Flow  Rate  of  the 

Cooling  Water,  R. . 

An  entha] 

py  balance  a 

round  the  whole  system  gives 

q£  -  CM 

^VVy 

=  R,  C  (T.),„   +  W,C  (5  ),   +  (L),0  CO*  . 

4  pv  yy         f  pv  o'3   s  '3  pv  f'3 

Since 

C^),,,  - 

(*f)j  ,  (?c) 

j  -  (?f)j  -  »3  .and  P  .  Wf  +  (L)j  , 

the  above  equation  oar.  be 

solved  for  the  cooling  water  flow  rate  as 

h  - 

q 

C       '  n 
p 

v1?  J 

(*f)j  - 

<V3 
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From   the   foregoing   equation   we   obtain 


<3s  1 


F         Cp        lT-       CF  1  .  (20) 
Co 


(Cf)3 




(Tf)3      "    (Tj)3  *    "    (Cf)3 

3-8.   The  Temperature  Difference  Between  the  Flashing  Brine  and  Non- 
Non-Flashing  Brine,  AT,  in  Sections  HR-n,  and  R-n. 

A)    A  T  in  Section  HR-n. 

In  Figure  6  a  whole  stage  in  the  infinite  stage  system  is  taken 
for  analysis.   The  brine  feed  stream,  F,  and  the  recycle  brine 
stream,  R^,  are  introduced  to  the  stage  from  the  right  and  leave  from 
the  left.   The  condensate,  W,  and  the  flashing  brine  stream  enter  the 
stage  from  the  left  and  leave  from  the  right.   The  temperature,  unit 
enthalpy,  and  quantity  of  each  stream  per  hour  are  denoted  in  the  fig- 
ure. 

By  making  an  energy  balance  around  this  stage,  the  following  re- 
lation was  obtained. 

(F  +  Rnifhj   +  dhj)  +  W  hc   +  L  hf 

=  (F  +  Rn)h-i   +  (W  +  dW)(hc  +  dhc)  +  (L   +  dL)(hf  +  dhf). 
Since  dhc  =  dhf,  -  dL  =  dW  and  L  +  W  =  F  +  Rn,  the  above  equation  becomes 
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(F  +  It  )dh.   =  W  dh  f  h  dV  +  Ldh..  -r  h.  dL 
v     n/   j         co        S    .     t 

=   (L  +  W)ah„  -  (h  ,  -  h  )dW 
:       i     i    c' 

-  '  (p  h  aa)dhf  -  (hf  -  ho)dw. 

(21) 

' 

Since  the  feed  brine  and  the  recycle  brine  receive  the  heat  of 

condensation  of  the  water  vapor,  we  have 

-  (?  +  B  )dh.  =  X  dW 

or             (F  +  R  )dh. 

-  dW  .      "   J  .   . 

(22) 

X 

Substituting  equation  (22)  into  equation  (21)  yields 

(hf  -  h  )(F  1-  H  )dh. 

/—l          T,      \   IV                1     —1                         7^      \   -1.                          ~                '-' 

(23) 

tr  +ri  Jan.  =  Ij  -r-n  Jan,  -r                         . 
v      r/   j    >      r/   i              X 

Cr.  rearranging,  we  obtain 

(1  -    -  ^      )  dh.  .  dhf. 

(24), 

Since 

h,  -  h 
1     °   c  c   - 

<  <  i, 

A 

equation  (2V)  can  be  approximated  by 

dh.  •  dh. 

(25) 

or 

dT  .  -  dT«, 

j 

(26) 
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Integrating  this  between  location  n'  and  n"  and  rearranging  gives 

(AT)n,   =  (Tf)n,  -   (Tj)n,   =  (Tf)n„   -  (Tj)n„   =  (AT)n„    (27) 

This  derivation  leads  to  the  conclusion  that  AT  is  nearly  constant  in 
the  heat  recovery  s.ections. 

B)  A  T  in  Section  R-n. 

In  the  heat  rejection  section,  a  derivation  similar  to  that  de- 
scribed above  leads  to  the  following  result, 

(F  +  Rn)dTf  =  (F  +  Rn+iJdXj.  (23) 

Since  both  F  +  Rn  and  F  +  Rn+1   are  constant,  the  equation  can  be 
integrated  between  location  n"  and  n  to  give 


(F  +  Rn)     (Tf)n.,   -  (Tf)n    -  (P  ♦  RN+1.)  |  (TjOn"   -  (Tj  )n{ 
Since  (Tf)n  =  (Tj)nn  ,  the  above  equation  becomes 

(F  -  Rn)  {(Tf)n„   -  (TjJn-i  =  (F*Rn+l)  j  <Tf>n   "  (Tj)nl   (29> 


or 

(F  +  Rn)(AT)n..  =  (F  +  Rn+1)   (  AT)n  .  (30) 

Therefore  the  temperature  difference  varies  gradually  from 

(AT)nn   to  (  AT)n  ,  as  is  shown  in  Figure  3. 
3-9.   Average  AT  in  the  n-th  Effect,  ATn 
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It  is  known  from  the  last  section  that  AT  is  constant  in  the 
heat  recovery  section  of  the  n7th  effect,  that  is,  (b.J)n,    =  (A  T)n„. 
In  the  heat  rejection  section,  however,  AX  changes  slightly  from 
(  AT)n.i  to  (  AT)n  according  to  equation  (30).   However,  since  the  heat 
recovery  section  contains  as  large  a  number  of  stages  as  the  heat  re- 
jection section  does,  we  may  reasonably  assume  that  the  AX  in  the  heat 
recovery  section  is  the  average  AT  in  the  effect.   Xherefore,  we  have 

AXn  =  (  AT)n,  -  (AT)n„  (31) 

By  applying  an  enthalpy  balance  between  locations  0  and  n',  we  obtain 

n=l 

%    +    <F  +  Rn>  CD<Ti>n.  =  '^   W,)C  (T  L.  +  <L)   C  (T  ) 
*       .n   pjn*    ■  -    i  p  c  n '      n'  p  i  n' 

Substituting  equations  (5),  (6),  and  (8)  into  the  above  equation  and 

rearranging  yields 

r      C  q  C 

(T,)n,    -    (T,)n,    (1    ♦  -JLJ-)    -  JLi  ♦    (X    - 


lf'n'         ^j'n' 


<Cf>n-l  F    Cp  <Cf>n-l 


[  (X,)  ,-11    )  ]. 

l,_v  f'n'  v  c'nM 


Since  (Tf)n,  -  (Xc)n,  =  <*n,  and  (Tf)n,  -  (T A,    »  AT^  the  foregoing 

equation  becomes 

V_l  +  <1  -    CF    )     d 

ATn=— J2 *^  ,   »-l.a.  3.  (3.) 

1  +   r"   F 
(Cf)n-1 

3-10.   The  Effective  At  for  Heat  Transfer  in  the  Brine  Heater, AtQ, 

and  in  the  n-th  Effect,  At  . 

A)  At  in  the  brine  heater,   AtQ 

If  we  let  Xs  be  the  steam  temperature,  then 

At  at  the  inlet  =  T   -  (T.)i  ■ 
s    v  j  '1 ' 

At  at  the  outlet  =  T£  -  (T-),, 
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(  AT)n.i  to  (  AT)n  according  to  equation  (30).   However,  since  the  heat 
recovery  section  contains  as  large  a  number  of  stages  as  the  heat  re- 
jection section  does,  we  may  reasonably  assume  that  the  AX  in  the  heat 
recovery  section  is  the  average  AT  in  the  effect.   Xherefore,  we  have 

AXn  =  (  AT)n,  -  (AT)n„  (31) 

By  applying  an  enthalpy  balance  between  locations  0  and  n',  we  obtain 

n=l 

%    +    <F  +  Rn>  CD<Ti>n.  =  '^   W,)C  (T  L.  +  <L)   C  (T  ) 
*       .n   pjn*    ■  -    i  p  c  n '      n'  p  i  n' 

Substituting  equations  (5),  (6),  and  (8)  into  the  above  equation  and 

rearranging  yields 

r      C  q  C 

(T,)n,    -    (T,)n,    (1    ♦  -JLJ-)    -  JLi  ♦    (X    - 


lf'n'         ^j'n' 


<Cf>n-l  F    Cp  <Cf>n-l 


[  (X,)  ,-11    )  ]. 

l,_v  f'n'  v  c'nM 


Since  (Tf)n,  -  (Xc)n,  =  <*n,  and  (Tf)n,  -  (T A,    »  AT^  the  foregoing 

equation  becomes 

V_l  +  <1  -    CF    )     d 

ATn=— J2 *^  ,   »-l.a.  3.  (3.) 

1  +   r"   F 
(Cf)n-1 

3-10.   The  Effective  At  for  Heat  Transfer  in  the  Brine  Heater, AtQ, 

and  in  the  n-th  Effect,  At  . 

A)  At  in  the  brine  heater,   AtQ 

If  we  let  Xs  be  the  steam  temperature,  then 

At  at  the  inlet  =  T   -  (T.)i  ■ 
s    v  j  '1 ' 

At  at  the  outlet  =  T£  -  (T-),, 
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and  the  average  At  for  the  heat  transfer  in  the  brine  heater  is 

AtQ  -,  T,   -  —   [(T.)^   ♦   (Tf)lt]  . 

Since 

A*l     uf'l.<      Uj'V   '•       ... 

the  above  equation  can  be  rearranged  to  give 

1 

1       C-/F                            r«i 

where 

l*f'0   *•  f;l' 

The  maximum  value  of  (T,)n  must  generally  be  limited  in  order  to 
2    0 

control  scale  formation.   In  this  study,  the  values  of  T   ar.c  (T.)- 

s     *  t   0 

are  assumed  to  be  fixed. 

B) .    At  in  the  r.-th  Effect,  it 

Referring  to  a  stage  ir.  section  .HR-n  of  Figure  5 ,  it  car.  be  seen 

that  the  effective  it  for  heat  transfer  in  the  infinite  stage 

operation  is  given  by          . 

uw  t  xy 
' —  =  uw  «  xy 

2 

?or  S  stage  operation,  the  effective  it  becomes 

vv  *-  xy 

2 
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because  the  condensation  temperature  is  constant  within  each  stage. 
Therefore,  the  loss  in  the  effective  it  for  heat  transfer  is 


uw  ^  xy      vv  +  xy 


In  other  words,  the  loss  in   At  for  heat  transfer  is  one-haHf  of 
the  temperature  drop  from  stage  to  stage.   Since  the  average 
temperature  drop  from  stage  to  stage  in  the  n-th  effect  is  given 
by 


(TA-i  -  ("A 


(3*) 


where  N  is  the  number  of  staged  in  the  n-th  effect,  the  average  loss 
in  effective  A  t   for  heat  transfer  is 


At 


n,  loss 


(T-)   -   -   (-0 
v  f'r.-i     v  fa 


2K 


(35) 


Therefore,    the   effective     At      is   given  by 


At        =    A  T        -     a  -      At        - 

a  n  •   n  n,    loss 


Hs  +      (1   - 

F  C 


ra~ 


2'S. 


(c-;. 


a  .  1,   2,  3. 


(36) 
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5-11.   lieat  Transfer  Area  Requirements  in  the  Brine  Heater  and 
the  n-th  Effect. 

Equations  for  the  heat  loads  and  temperature  difference  for  heat 
transfer  in  the  brine  heater  At   and  that  in  the  n-th  effect,   At   , 
have  been  developed  in  sections  i-k   and  i-10,    respectively.  These  are 
used  to  calculate  the  heat  transfer  areas  by  the  equation 


(37) 


u  U  t) 

By  assuming  U  is  constant  in  the  brine  heater  and  each  effect, 
and  substituting  equation  (53)  into  equation  (37),  the  heat  transfer 
area  in  the  brine  heater  A0  ,  can  be  obtained  as 


(3&) 


1    Ut         ] 


U  \     T   -  (T,)n   +  — 

Substituting  equations  (ll)  and  (36)  into  equation  (37)  yields 

the  heat  transfer  area  in  the  n-th  effect,  A   ,  as 

n 


(c„),_.    (cfL  J 


X 


a     1  C, 


-«„   "  . > 

r  0  2H    '  :    ) 

n   =   1,  2,  3-  (39) 
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3-12.   Pumping  Head  for  the  Recirculation  Pump,  J  ■ . 

n 

The  recirculation  pump,  J  ,  takes  in  recycle  brine  8   of 

concentration  (0f)   at  temperature  (I.)  ,  pressurizes  it  to  a  pressure 

sufficiently  high  so  that  the  recycle  brine  does  not  soil  within  the 

heating  tubes.   The  highest  temperature  to  which  this  brine  stream  is 

heated  is  (T,)   .  •  Thus,  the  Dumping  head  (  *p  )/D 

-  n.-i  *    n"  M  recuirod  for  J 

is 

can  be  evaluated  as, 


U  P)„     1 


U 


^  ff> 


^  ;n-l  v  ;r. 


J' 


(40) 


where  r}-  is  the  fractional  excess  pumping  head  required  due  to  fi 

losses  and  (P)   is  the  vapor  pressure  of  the  brine  at  concentrat: 

(C-)   and  temperature  (T,)  . 
*  r '.n  I  'n 

The  vapor  pressure  of  brine  at  a  given  temperature  is  less 
vapor  pressure  of  pure  water  at  the  same  temperature  due  to  the 
pressure  depression'  of  the  solution.   Thus, 


lotion 

on 


tr.ar.  trie 
vaoor 


(*)B  -  (P°)B  -(b), 


»nd 


'a-1 


(P°)r 


I  »  ) 


n-1  ' 


where  (p  )n  and  (P  )   ,  are  the  vapor  pressures  of  pure  water  at 
temperatures  (Tf)n  and  (Tf)n_1  respectively  and  (  g  )   and  (  g )  _  ar 
vapor  pressure  depressions  of  the  brine  streams 
and  (lf)B_j  an°  concentrations  (Cf)   and  (C,)      .    ,    respectively 


temperatures  (T„,) 
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If  wo  assume  that 


wg  obtain 


and  equation  ('t-0)  becomes 


i  +  rtf)  [(P°)n_i   -  (p\] 


Within  the  operating  temperature  range  of  the  MEMS  system,  the 
vapor  pressure  of  water  may  be  represented  by 


i-i  P 


RT 


where  \    is  the  latent  heat  of  vaporization  and  D  is  an  integration 

constant.   On  rearranging,  we  obtain 

_   "  RT  RT 

r   =  e   e       ■  I  I 

where  E  ■=  eJ.   By  assuming  a  constant  value  of  1000  Btu/lb   for  A 

and  from  the  steam  table  at  two  temperatures,  3  is  evaluated  to  have 

q 
a  value  of  1.523  x  10   lb„/f^2  . 
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Substituting  equation  (^2)  into  equation  (41)  g 


gives 


A  P 


(1  tn,)   )  axp(-  T^r-Y )   -  -  exp( 

/        V  f;n-l 


v  in 


X i 


n  =   1,  2,    $. 


C'5) 
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CHAPTER  4 
ECONOMIC  ANALYSIS 

The  water  production  cost  depends  primarily  on  the  capital  cost 
and  operating  cost  of  the  plant.   Only  those  costs  which  are  most 
significantly  affected  by  changes  in  the  design  variables  are  con- 
sidered in  this  study.   The  cost  of  the  plant  site,  labor  cost,  over- 
head costs  and  insurance  costs  are  not  considered  here  as  they  are 
little  affected  by  changes  in  the  design  variables. 
The  capital  cost  consists  of  three  items: 
(a).  The  heat  transfer  area  cost, 
(b).   The  recirculation  pump  cost, 
(c).   The  outer  shell  cost. 
The  operating  cost  consists  of  four  items: 
(a).   Feed  brine  cost, 
(b).   Cooling  water  cost, 
(c) .   Steam  cost, 

(d).   Power  cost  for  recirculation  pumping  in  each  effect. 
Each  cost  item  is  expressed  as  the  cost  per  1000  gallons  of  fresh 
water  produced  in  the  whole  plant.  .The  following  notation  is  used 
to  represent  the  various  cost  items: 
Ex  =  Steam  cost, 

E2  =  Capital  cost  of  brine  heater, 

E3  =  Capital  cost  of  the  heat  transfer  area  in  the  n-th  effect, 
E4  =  Power  cost  for  recirculation  pumping  in  the  n-th  effect, 
£5  "  Capital  cost  for  recirculation  pump  in  the  n-th  effect, 
E&  =  Capital  cost  of  the  outer  shell, 
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27  ' 

Feed  brine  cost, 

^8  = 

Cooling  water  cost. 

4-1.  Ih« 

Capital  Costs. 

"he 

capital  costs  of  the  principal  items  of  equipment  are  evaluated 

per  unit 

of  production  in  a  unit  time.  The. annual  capitalization  charges 

for  these 

equipment  items  are  calculated  at  0.074  of  the  initial  cost 

per  year, 

as  recommended  in  the  Office  of  Saline  Water  procedures  (10). 

A  load  fa 

ctor  of  JjO  on-streara  days  per  year  will  be  assumed.   Therefore, 

the  oapit 

alization  charge,  ii  ,   is  J. 4x10"   of  the  initial  cost  per 

hour  on-s 

tream. 

(a)   Brine  Heater  Cost,  E„ 

The 

brine  heater  cost  is  assumed  to  be  proportional  to  the  brine 

heater  ar 

ea,  A  ,  which  is  given  by  equation  (38).  Therefore,  the 

capital  cost  per  100C  gallons  of  water  production  per  hour,  Z    ,    is 

given  by 

32 

*CA„W 
3  u 

W 
n 

"  °ht   ?                      ,  /,                     (<*) 

TT              f   n   >                                         .:,.    1     1 

0,(1+^).         (Cf)3 

where 

C- 

m     the  capital  cost  per  unit  of  heat  transfer  area  for 

the  brine  heater, 

Ght 

=   *CB  ' 

w 

mass  equivalent  to  1000  gallons  of  water, 

a 

>  T   ,  i.e.,  the  steam   temperature. 

This  COS1 

;  includes  both  the  heat  transfer  area  and  the  outer  shell  costs. 
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(b)   Heat  Transfer  Area  Cost  in  the  n  th  Effect,  E-. 

The  heat  transfer  area  for  the  n-th  effect,  A   ,  is  given  by 

n 

equation  (39)-  Therefore,  the  equation  for  the  capital  cost  per 
1000  gallons  of  water  produced  can  be  represented  by 


"3 


c„ 


<Cf)n-l    "   ^P7 


— 7T"   +  a  _  *.  1 


C_ 


n  ^    rc7T^j 


(*,)„_,   -   (*,), 


+  ^7v7 


1  Tc3 


l,  2,  3,      (*5) 


f'3 


C..  -  the  capital  cost  per  unit  of  heat  transfer  area, 
a 

Ccd  "  *CH  • 
(c)  Recirculation  Pump  Cost  in  the  n-th  Effect,  E   . 


The  pump. cost  is  assumed  to  be  proportional  to  its  power  rate.' 
Since  the  pumping  head  in  the  n-th  effect  is  given  by  equation  (V5), 
the  cost  equation  for  the  pump  J^  can  be  written  as 
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(c)  Recirculation  Pump  Cost  in  the  n-th  Effect,  E   . 
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*S 

=   vf  C,    A  Pn    Rn     W 

.,.  .    B  (1  *  Tf )   rnCFW   f           A                 ^     ) 
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CHAPTER  5. 

OPTIMIZATION 

The  performance  equations  of  the  XEXS  process  are  described  in  chapter 
3   and  the  cost  equations  of  the  water  production  are  derived  in  chapter  k. 
Armed  with  these  equations  we  can  proceed  by— ueiftg— t4-K>— o^iirri-z-art-.ro-ii--w^c4^*4<H^— 
to  determine  the  optimal  conditionsof  the  process.  Since  a  discrete  form  of 
the  maximum  principle  is  effective  for  seeking  the  optimal  condition; of  a 
sequential  multistage  multidecision  process,  we  shall  use  it  in  conduction 
with  search  techniques  to  optimize  the  process  (11).   Two  search  techniques 
are  used  here:   one  is  the  parametric  search  and  the  other  is  the  simplex 
method.   The  results  of  the  numerical  solution  and  the  comparision  of  the  two 
approaches  are  given  in  sections  5-3,  5-4,  and  5-5.  The  computer  programs. 
for  the  two  methods  are  presented  in  the  Appendix. 
5-1.   The  Performance  Equations 

From  equation  (51)  it  is  known  that  the  water  cost  S  is  a  function  of 

thirteen  variables:   qs/F,  C?,    (Cf)1>  (Cf),,  (Cf)„,  r^  r2>  ly    (Tf)0>  (-.-);, 

(T  ).,  (T  ),,  and  (T.)  .   But  C  ,  the  brine  concentration  of  sea  water  feed, 
L    £  Z   3  J  J         r 

is  fixed  and  assumed  to  be  3.5%  (5),  and  (T.)    the  temperature  of  the  brine 
stream  leaving  the  brine  heater,  is  fixed  because  of  the  need  for  the  controling 
of  the  scale  formation.  The  sea  water  temperature  (T.)„  is  rarely  changed 
and  is  assumed  constant.  Therefore  we  have 

S  =  S(qs/F,  (Cf)r  (Cf)2>  (Cf)3,  (Tf)1,  (Tf)2,  (Tf),,  r...  r,,  r.)  .     (52) 

However,  these  ten  variables  are  not  all  independent;  equation  (15)  gives 
three  relations  between  these  variables.  We  have  then  seven  independent  vari- 
ables.  According  to  the  maximum  principle  we  classified  the  variables  into 
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given  in  advance  before  we  proceed  to  optimize  the  cost  function. 
The  optimization  problem  we  have  imposed  is  as  follows : 
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x2  with  X]_   and  qs/F  preassigned. 
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Once  the  values  of  x±   and  qs/F  are  known,  the  optimal  value  of  Q   is 

sought  by  the  algorithm  of  the  maximum  principle,  but  the  optimal 

3 
values  of  x1  and  qs/F  must  be  found  by  one  of  the  two  search  techniques 

mentioned  before. 


5-2.   Search  for  Optimum  by  the  Maximum  Principle 
(a)   Differentiation  of  State  Variables. 

The  differentiations  of  the  state  variables  with  respect  to  the 
decision  variables  and  state  variables  are  given  below.   These  are 
used  to  determine  the  adjoint  variable  and  the  derivatives  of  the 
Hamiltonian  funcitons  in  the  following  sections. 
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(b>  Adjoint  Variables  8. 

Since 

=1  =  0.      -2  - 1> 

c3  =  0, 

we  can  write 

•J  =  o,   ,3  =  1, 

-J.o. 

3 
However,  since  x^  is  prefixed, 

_ 

3  , 
21  ^  Cl  • 

' 

3 
Then  H  becomes 
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(e)   Calculatior 
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Since  x1  and  q  /F  ar 
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values  or  x^  and 
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from  the  following  procedures: 
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Step  S. 

3 
Calculate  z^  from  equation  (91)  and  equations  (68)  and  (78), 

2    2    2    111 
and  calculate  z^ ,  z2 ,  z3,  z^  ,  z2 ,  z3  from  equations  (92) 

through  (97). 

Step  6. 

Calculate  «»'   ,H*    '«*   >  »2   *  H*   from  equations  (98) 

^    '63    >8~   >e|.  -.'Bi 
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then  the  assumed  S^  are  the  optimal  values;  otherwise  go  to 

the  next  step. 

Step  8. 

If  xi   ^s  greater  than  that  computed  in  the  preceding  itera- 

tion, then  one  half  of  the  original  aQ^  is  used;  otherwise 

the  original  aQ  ™  is  used. 

Step  9. 

The  new  set  of  decision  variables  (0")new  is  obtained  by 

<Onew  =  C^old  -*»J                              <103> 

When 

— „•  >  0    use  (-)  sign 
*»i 
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to  step  2  and  repeat  the  computation  until  the  optimum  is  obtained. 
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references  (10)  and  (12).   A  computer  flow  chart  based  on 

the  procedure  we  have  described  in  part  (e)  is 

given  Fig.  7. 

Table  1 .   Numerical  Values  for  the 

Constants 

Symbols 

Explanation 

Numerical  values 

a  or  T£ 

Steam  temperature 

274. 4°F 

B 

Coefficient  of  the  Clausieus-Clapeyron 

equation 

1.79  x  10   lbf/ft 

Cc 

Unit  cost  of  cooling  water 

4.4875  x  10"7  $/lb 

Cf 

Concentration  of  sea  water  feed 

0.035  wt.  fraction 

Cht 

Unit  cost  of  brine  heater 

3.76  x  10"5  S/ft2 

°P 

Heat  capacity  of  sea  water 

1.0  Btu/lb°F 

Cs 

Unit  cost  of  steam 

2.5  x  lO"4  $/lb 

Ccd 

Unit  cost  of  condensing  area 

2.397  x  10"5  $/ft2 

Sp 

Unit  cost  of  pump  and  pumping  power 

2.903  x  10"9  $/ft-lb 

Pc 

Unit  cost  of  feed  pretreatment 

1.1795  x  10"6  $/lb 

Nn 

No.  of  stages  in  n-th  effect 

23,  23,  22 

U 

Overall  heat  transfer  coefficient 

510  Btu/hr.  ft2  F 

A 

Latent  heat  of  flash  brine 

1000  Btu/lb 

As 

Latent  heat  of  steam  at  274.4  F 

and  45  psia 

928.9  Btu/lb 

R 

Ideal  gas  constant 

0.1104  Btu/lb. °F 

e 

Density  of  brine 

62.5  lb/ft3 

0 
*3 

Temperature  of  brine  entering  the 

first  effect 

250°F 
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5-3.   Parametric  Search  for  Overall  Optimum 

A  set  of  grid  points  of  x^  and  q  /F  is  shown  in  Table  2.   For  each 
grid  point  the  scheme  we  have  described  in  the  last  section  was  used 
to  seek  the  optimum.   Since  corresponding  to  a  given  xj  we  have  a  set 

of  qs/F>  we  can  use  a  graphical  method  to  find  an  overall  optimum  for 

3 
a  given  x  .   These  data  are  given  in  Table  3  and  the  corresponding 

figures  are  shown  in  Figs.  8  through  15.   The  optimal  policies  for 

3 
each  x-^  are  plotted  in  Fig.  16,  from  which  the  optimal  condition  of  the 

whole  system  was  obtained. 

3 
In  Table  3,  for  given  xj  and  qs/F,  the  optimal  values  of  rn  and 

(T^)n,  and  the  water  cost,  C,  are  tabulated.   For  example,  for  x1  =  0.05 

and  qs/F  =  16,  the  optimal  policy  is 

x1   =   0.99    r2  =  1.30    r3  =  1.72 

(Tf )x  =  198°F    (Tf )2  =  151°F    (Tf)3   =  102°F 

C  =  0.2907$/1000  gal. 

In  Figs.  8  through  15,  these  optimal  policies  are  plotted  against 

3 
qs/F  for  a  fixed  x± . 

From  Fig.  16,  the  overall  minimum  water  production  cost  is  found 

to  be  0.2855$/1000  gal.,  when  the  system  is  operated  under  the 

following  conditions: 

salt  concentration  of  the  flashing  brine 

3 
leaving  the  third  effect,  x^   =  0.065, 

the  ratio  of  heat  load  to  seawater  feed  a.s/F  =  27, 

recycle  ratio  in  the  first  effect,  rj.  =  2.14, 

recycle  ratio  in  the  second  effect  r2  =  2.88, 

recycle  ratio  in  the  third  effect,  r3  =  3.86, 
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j  r.»3.6l,r2=5.38,r,»7.60 

m|                                         ■    |  (?,) 3-189, (?,) 2=i-2, {t..;5=io;.      ! 

j   0=0.2909    $/1000    gal. 

:   ri=i.jo,r£  =  -..o9,r3  =  o.9T                 j  r.=3.o3,    r2«5  .  52,r3=7  .59 

-5      ;"f):=i5^,;,:-fi2=ii.T,(T::;)i=io3     (2j)j  =  159,  (i,);  =  1-2,  (?_.;■;  =105 
S-0.2962    $/1000    gal.                          |   0=0.293-    S/1000    gal. 

s 
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[\4                                 1    "'"   '                             i 

30  |  ;?J:=192,(T..),  ^^.('rp;  =23 

'   C  =  0.  3021    S/1000    /ral. 

- 

.3  ~-,-3.>G,r  =5.G5,r~=7.-2 

3J,  i                                                -]  (sji1-i9i,ter)2«oAo,-(sjj -100 

J  0=0.2936  A/1CGC   fjal.                                         1 

1   ",-.;. .50,  !--•*. 5, <,i-_  =  c. jo 

*                ■    -*                                  1 
35  ]   (?;,).  =19i-,(T,V  ■l^V.d*).   =  102   1 
0=0.2395    S/1000    gal. 

|  r. =3 .29  ,r2*4 . o5 ,r_»6. 50 

36  j   (T..X  =  192,(2,),  =lit31(Tl,)J,=lC3   j 
;               !   0=0.2539    S/1000    gal. 

'   -"1=2.33,r,  =  s.60,r_  =  6.'*S 
„  1                          '                   i                               j 

J'  i  (T'..)n  -195, (t;)  9-lH,(.J->  -103.  1                                                               ! 

j               .   0=0.2335    $/1000    gal.                          j 

- 

;  r1  =  3.33,r2=N,.6'.,r,  =  6.59                    ^-3.61,^-5.17.^-7.35 

38  |   nC.,).  =19^>(Tf)2=U*l.,(Tf):Ji=10li       (2f)l-191>(T-)2-liH,(?    $,  -102 

;   0=0.2533    S/IOCG    gal.                            j  0=0.2913    S/-000    gal. 

- 

1       39  i 

r1=3.62,r.2=5.20,r,  =  T.-2 

(T.,).  -191, (Tf )    -11*1, (I;),  =103 
0=0.2905    5/1000    Bal. 

]  r,=3.6o,r   =5.31,r,=?.5o                      [ 

'        fc0|                                                  '                      (»f)1"l89.l»,)j'l*i.(«;^"3 

1   0=0.2907    3/1000    gal. 

j  r.»3.6l,r2=5.38,r,»7.60 

m|                                         ■    |  (?,) 3-189, (?,) 2=i-2, {t..;5=io;.      ! 

j   0=0.2909    $/1000    gal. 

:   ri=i.jo,r£  =  -..o9,r3  =  o.9T                 j  r.=3.o3,    r2«5  .  52,r3=7  .59 

-5      ;"f):=i5^,;,:-fi2=ii.T,(T::;)i=io3     (2j)j  =  159,  (i,);  =  1-2,  (?_.;■;  =105 
S-0.2962    $/1000    gal.                          |   0=0.293-    S/1000    gal. 
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THE  OPTIMUM  POLICIES  FOR  VARICUO  x?  and  cj? 

s>^ 

0.12 

« 

r.=j.Su,r,=3.>v,r,-| .  yj   • 

(T.J,  =137, (T.J,  =137,  (T./),  =55 

-  -       -  •*      1 
0=0.3055  S/1000  gal. 

|                             r1--..0rTro-S^5-,r.,  =  S.&5 

3T   i 

i  (T  )  -188, (T;),  "136, (tA  =55 

i                         i-l       -  ^       -  ; 
, !  0  =  0.3073  S/iCCG  gal. 

1.0 

23                                      j  . 
(T:).  =107, (T.J  =138, (Tj,  =102. 

r  1     .  r  1.      ,        t  ?     . 

0=0.25^1  S/1000  gal.      /        < 

• 

| r, =3. 07, r, =5. 75, r, =8. 31 
111 

i  (T.)]  =lSo)(-1^,)2=135,(T..)5=10.3 

\   C=0 . 2532  S/1000  gal. 

.  '  ■ 

1 

' 

:  r.=i.So,r=5.£3,r.i  =  6.5G 
1*2   1  ~                 - 

|  (^).=l37,(T..)5=135,(Tf)./=103 
I  0  =  0.2525  $/100d  gal.      ■ 

r.=3.6o,r  =5.32, r  =6.55 
1  x         i.                     i 

"3   ..(T:;).=l85,(T-)9  =  133>(TJ)J  =  10ii 

! 0=0.2531  S/1000  gal. 

r1-lt.ll,r2-S.16,r3=9.35 

(Tf)  =137,  (T.,)  =136,(  2f),  =102 
0=0.2953  S/1000  gal. 

;r1  =  4.12,  r,=3.22,  1-^  =  5.35 
|  C=>0.  29U8  S/1000  gal. 

11. 

1      I 

"5  ; 

r1=!,.ll,,r,  =  6.22>r,  =  5.-.l 

(T^.)-=1oc;)vTj)?  =  137,(T_,  }i=10^ 
-*          -          ^ 

0=0. 25^5  S/1000  gal. 

!   -o   j                              "  '■■' 

r1«'».15,r2»6.25,r31B9.50 

(TJ).=l37,(TJ)0  =  137,(T..;  j  =  10, 
•0  =  0.2952'  5/1000  gal. 

:ri=2.53,r£  =  6.0i,r3  =  a.Sl. 

1*8   |.(5p.»l88,tT J  «ll»0,{Tf)  <L07  i 

0  =  0. 255o  S/1000  .gal. 

50                               (Tp1_-»l68,(Tp2i»133,(Ti.)5=106 

|0=0.2570  0/1000  gal. 
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THE  OPTIMUM  POLICIES  FOR  VARICUO  x?  and  cj? 

s>^ 

0.12 

« 

r.=j.Su,r,=3.>v,r,-| .  yj   • 

(T.J,  =137, (T.J,  =137,  (T./),  =55 

-  -       -  •*      1 
0=0.3055  S/1000  gal. 

|                             r1--..0rTro-S^5-,r.,  =  S.&5 

3T   i 

i  (T  )  -188, (T;),  "136, (tA  =55 

i                         i-l       -  ^       -  ; 
, !  0  =  0.3073  S/iCCG  gal. 

1.0 

23                                      j  . 
(T:).  =107, (T.J  =138, (Tj,  =102. 

r  1     .  r  1.      ,        t  ?     . 

0=0.25^1  S/1000  gal.      /        < 

• 

| r, =3. 07, r, =5. 75, r, =8. 31 
111 

i  (T.)]  =lSo)(-1^,)2=135,(T..)5=10.3 

\   C=0 . 2532  S/1000  gal. 

.  '  ■ 

1 

' 

:  r.=i.So,r=5.£3,r.i  =  6.5G 
1*2   1  ~                 - 

|  (^).=l37,(T..)5=135,(Tf)./=103 
I  0  =  0.2525  $/100d  gal.      ■ 

r.=3.6o,r  =5.32, r  =6.55 
1  x         i.                     i 

"3   ..(T:;).=l85,(T-)9  =  133>(TJ)J  =  10ii 

! 0=0.2531  S/1000  gal. 

r1-lt.ll,r2-S.16,r3=9.35 

(Tf)  =137,  (T.,)  =136,(  2f),  =102 
0=0.2953  S/1000  gal. 

;r1  =  4.12,  r,=3.22,  1-^  =  5.35 
|  C=>0.  29U8  S/1000  gal. 

11. 

1      I 

"5  ; 

r1=!,.ll,,r,  =  6.22>r,  =  5.-.l 

(T^.)-=1oc;)vTj)?  =  137,(T_,  }i=10^ 
-*          -          ^ 

0=0. 25^5  S/1000  gal. 

!   -o   j                              "  '■■' 

r1«'».15,r2»6.25,r31B9.50 

(TJ).=l37,(TJ)0  =  137,(T..;  j  =  10, 
•0  =  0.2952'  5/1000  gal. 

:ri=2.53,r£  =  6.0i,r3  =  a.Sl. 

1*8   |.(5p.»l88,tT J  «ll»0,{Tf)  <L07  i 

0  =  0. 255o  S/1000  .gal. 

50                               (Tp1_-»l68,(Tp2i»133,(Ti.)5=106 

|0=0.2570  0/1000  gal. 
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temperature  of  the  flashing  brine  leaving  the 

first  effect,  (Tf)l  =  196°F, 

temperature  of  the  flashing  brine  leaving  the 

second  effect,  (Tf)2  =  148°F, 

temperature  of  the  flashing  brine  leaving  the 

third  effect,  (Tf >3  =  103°F. 

5-4.   Overall  Optimum  by  the  Simplex  Method 

It  has  been  shown  that  the  discrete  version  of  the  maximum  prin- 
ciple can  be  applied  to  find  the  optimal  condition  for  a  sub-optimization 
problem  with  a  set  of  given  values  for  x^  and  qs/F.   If  a  multi- 
dimensional search  technique  is  combined  with  the  sub-optimization 
procedure  using  the  maximum  principle  to  minimize  the  objective 
function  depending  on  the  two  variables,  x   and  q  /F,  which  are  fixed 
in  each  sub-optimization  step,  the  overall  optimal  policy  of  the 
system  can  be  obtained  in  a  straight-forward  manner  by  using  one 
complete  computer  program. 

A  number  of  multi-dimensional  search  techniques  are  available, 
such  as  Powell's  method  (13),  Box's  method  (16),  Smith's  method  (15), 
etc.   Powells  method  is  known  to  be  an  efficient  method  for  finding 
the  minimum  of  an  objective  function  or  simply  a  function  without 
calculating  its  derivatives.   However,  Nelder  and  Mead  (14)  have  re- 
cently developed  the  so-called  "simplex  method"  which  is  reported  to 
perform  more  efficiently  than  Powell's  method.   It  is  said  that,  for 
two-dimensional  search  problems,  the  efficiency  of  the  simplex  method 
far  exceeds  that  of  Powell's  method.   The  simplex  method  was  used  in 
this  study. 

The  general  concept  of  this  method  for  the  minimization  of  a  function  oj 
n  variables  is  to  set  up  a  simplex  of  (n+1)  vertices,  that  is,  to  select  (n+1) 
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points  in  the  space  of  a   variables  and  calculate  values  or  the  function  at 
the  selected  points.   Then,  by  comparing  the  calculated  values  of  the  function 
among  themselves,  the  vertex  with  the  highest  value  (i.e.  the  worst  point  in 
minimization)  is  replaced  by  another  point  with  a  lower  value  of  the  function, 
which  is  determined  according  to  certain  operations  to  be  described  later. 
The  simplex  method  forces  the  function  to  approach  the  minimum  by,    at  each 
stage  of  operation,  discarding  the  worst  point  of  a  simplex  and  adapting  a 
better  point  to  form  a  new  simplex.  This  procedure  is  repeated  until  z:.& 
minimum  point  is  achieved. 

For  the  problem  in  hand,  the  simplex  can  be  represented  by  a  triangle 
as  shown  in  Fig.  17.   ?1,  P2>  and  ?3  are  the  points  in  the  two  dimensional 
space  of  x1  and  q  /F,  which  define  the  current  "simplex." 
We  define 

yn  =  the  value  of  the  objective  function  or  the  water  cost  x,  at 

the  point,  P  , 
r    '  n 

Pi  ■  the  vertex  with  the  lowest  value  of  the  objective  function  (y  ) 

in  the  simplex, 
P3  =  the  vertex  with  the  highest  value  of  the  objective  function  (y,) 

in  the  simplex, 
?,  ■  the  vertex  at  which  the  corresponding  value  of  the  objective 

function  (y_)  lies  between  (y  )  and  (y  ) , 
?^  ■  the  centroid  of  the  vertices,  P  and  P  ,  with  the  value  of  the 

objective  function  (y , ) . 

The  operations,  through  which  a  new  point  with  a  lower  value  of  the 
objective  function  is  found,  are  reflection,  expansion  and  contraction.   The 
reflection  of  the  worst  point,  P.,  with  respect  to  centroid,  ?   is  denoted 
by  ?j  and  its  co-ordinates  are  defined  by  the  relation 
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P,  -  P.  +  a'(P,  - 
0    4       4 

V 

(104) 

whe 

re  a'  is  a 

positive  constant,  the  reflection  coefficient.  Thus  P 

is  on 

the 

line  joining  P,  and  P., 
J       4 

on  the  far  side  of  P,  fron  ?  with  ?  ?  = 
4       3       5  4 

a'?.?,. 
J  4 

The  refle 

cted  point  P 

nay  be  expanded  to  P,  by  the  relation 

P6  = 

p4  +  Y'(P5  - 

V 

(105) 

The 

expansion 

coefficient  y 

' ,  which  is  greater  than  unity,  is  the  rat: 

:.o   of 

the 

distance  P 

~.    to  P.P,  . 
o  4     J  4 

The  contr 

action  of  the 

worst  point,  ?,,  with  respect  to  the  centroid, 

P4' 

is  represe: 

nted  byP'  and 

defined  by  the  relation 

P6  = 

?4  +  3'(P3  - 

V 

(106) 

where  S'  is  a  positive  number  between  0  and  1  and  is  the  ratio  of  the 

distance 

pT?T 

o  4 

■  to  P3P4. 

The  values  of 

these  coefficients  considered  best  by  N'el 

der  and 

Xead 

(14)  are 

a'  = 

1,      3'  = 

2>          and  y'    -2. 

The 

details  of 

the  procedure 

for  using  the  method  are  described  as  fol 

lows: 

First,  P 

is  reflected 

to  P.,  and  if  y  lies  between  y1    and  y„,  t 

r.en  ?; 

is  r 

eplaced  by 

Pj  and  we  start  the  procedure  again  with  a  new  simplex. 

if  y5  <  y: 

,  that  is,  if 

the  reflection  has  produced  a  now  mini 

,'  Chen 

we  e; 

spand  P.  to 

P£.   If  y,  < 
o      J 6 

y-p  we  replace  P,  'ay   ?,  and  restart  the  p: 

roccij . 

Zuz   . 

if  H  >  *v 

then  we  have 

a  failed  expansion,  and  we  replace  ?  by  ?.  before 

restarting. 
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r.en  ?; 

is  r 

eplaced  by 

Pj  and  we  start  the  procedure  again  with  a  new  simplex. 

if  y5  <  y: 

,  that  is,  if 

the  reflection  has  produced  a  now  mini 

,'  Chen 

we  e; 

spand  P.  to 

P£.   If  y,  < 
o      J 6 

y-p  we  replace  P,  'ay   ?,  and  restart  the  p: 

roccij . 

Zuz   . 

if  H  >  *v 

then  we  have 

a  failed  expansion,  and  we  replace  ?  by  ?.  before 

restarting. 
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If,  after  reflection,  we  find  that  y5  >  y^  and  y_  >  y2 ,  then  we 
define  a  new  P,  to  be  either  the  old  Po  or  Pr,  depending  on  whichever 
has  the  lower  yn  value,  and  then  contract  P-  to  PJ.   We  then  accept  Pi 
for  P«  and  restart  the  procedure,  unless  yl>y~y    that  is,  the  contracted 

point  is  worse  than  P,.   For  such  a  failed  contraction,  we  replace  P2 

(P,  +  P,)      (P,  +  P,  )       . 
and  P3  by  and  —3 ±—  respectively  and  restart  the  process. 

A  flow  diagram  of  the  method  is  given  in  Fig.  18,  and  a  complete 
computer  program  for  the  simplex  method  together  with  the  sub-optimization 
program  by  means  of  the  discrete  form  of  the  maximum  principle  is  given 
in  Table  A3  of  the  Appendix. 

The  optimal  water  production  cost  obtained  by  using  this  method  is 
$0. 2855/1000  gal.  and  the  corresponding  optimal  operating  conditions 
are  as  follows: 

salt  concentration  of  the  flashing  brine 

3 
leaving  the  third  effect,  x^  =  0.065, 

the  ratio  of  heat  load  to  seawater  feed       °,S/F  =  27> 

recycle  ratio  in  the  first  effect,  r-^  =  2.139, 

recycle  ratio  in  the  second  effect  r2  =  2.877, 

recycle  ratio  in  the  third  effect  r   =  3.861, 

temperature  of  the  flashing  brine  leaving 

the  first  effect,  (Tf)l  =  195.9  F, 

temperature  of  the  flashing  brine  leaving 

the  second  effect,  ^Tf^2  =  l47-7°F> 

temperature  of  the  flashing  brine  leaving 

the  third  effect,  <Tf'3  =  103°F* 

The  same  procedures  used  here  for  the  two  dimensional  search  can  be  extended 
to  the  n-dimensional  problem  (14).   The  worst  point  of  a  simplex  with  (n+1) 
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vertices  is  reflected,  expanded  or  contracted  in  the  same  manner  with  respect 
to  the  centroid  of  the  remaining  n  vertices  until  the  minimum  point  is  attair.ee. 

5-5.   Comparison  of  Results  from  the  Use  of  the  Two  Search  Techniques 

The  optimal  policies  from  the  two  search  techniques  are  tabulated  in  Table 
4.   Both  in  the  sub-optimization  stage  and  in  the  search  stage,  a  criterion  is 
adopted  to  test  if  the  -minimum  point  of  the  objective  function  is  attained  and 
ho  further  iteration  is  needed.   Theoretically  when  the  objective  functions 

attain  their  minimum  points,  it  is  necessary  that  the  values  of  the  deriva- 

. ,    .    -    .     3K      ,     .     .  . 
tives  or  the  corresponding  namiitoman  functions,  tt»  m  tne  sub-optimization 

o  a 

step  and  the  "standard  error"  defined  by 


Ik  (yi " y^  ' 3 


where  y  ,  i  =  1,  2,  3,  and  y,  are  the  values  of  the  objective  function  at  tt.e 

vertices  and  centroid  of  the  simplex  respectively,  in  the  search  step  are  both 

equal  to  zero.   In  the  actual  calculation,  these  values  are  compared  to  some 

pre-set  values  or  criteria  and  the  iteration  stops,  when  they  fall  below  such 

criteria.   In  the  computer  code  developed,  the  criterion  for  tt  is  desigr.atec 

o  o 

by  ER  and  that  for  the  "standard  error"  by  ERROR. 

The  numerical  results  in  column  (a)  of  Table  4  from  the  parametric  search 

-4 
are  obtained  by  setting  ER  =  1  x  10   ;  in  column  (b)  from  the  simplex  method 

by  setting  ER  =  1  x  10   and  ERROR  =  1  x  10   ;  in  column  (c)  from  the  simplex 

method  by  setting  ER  -=  0.5  x  10   and  ERROR  =  0.5  x  10"  . 

From  the  closeness  of  the  numerical  values  of  the  various  operating 

variables  between  columns  (a)  and  (b) ,  it  can  be  concluded  that  the  parametric 

search  and  the  simplex  method  both  lead  to  the  same  optimum  point  if  the  same 

criterion  is  used.   However,  for  the  parametric  search,  tedious  exhaustive 

numerical  and/or  graphical  search  procedures  must  be  carried  out.   But  by  using 
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Table  4.   Optimal  Policy,  Operating  Conditions,  and  Cost 
from  the  Two  Search  Techniques 


Item 

Symbol 

Unit 

;  Parametric 
Search 

(a) 

Sirapl 

2x  .Method 

(b) 

!.) 

Maximum  allowable 
error 

ER 

lxlO"4 

lxl  0"4 

!                  41 
;  O.SxlO-4  i 

ERROR 

lxlO"4 

- 

_4 
0.5x10 

exit  brine  cone, 
of  the  3rd  effect 

3 
*1 

wt . f rac . 

0.065 

0.065 

0.06475 

ratio  of  heat  load 
to  seawater  feed 

qs/F 

Btu/lb 

27.0 

27.0 

26.64 

recycle  ratio  in 
the  1st  effect 

rl 

2.14 

2.139 

2.125 

recycle  ratio  in 
the  2nd  effect 

r2 

2.88 

2.877 

recycle  ratio  in 
the  3rd  effect 

r 
r3 

3.86 

3.861 

3.819 

exit  brine  temp, 
of  the  1st  effect 

(Tfh 

°F 

196.0 

195.9 

195.85 

exit  brine  temp, 
of  the  2nd  effect 

<Tf)2 

°F 

148.0 

147.7 

147.72 

exit  brine  temp, 
of  the  3rd  effect 

(Tf)3 

°F 

103.0 

103.0 

103.0 

water  production 
cost 

3 
*2 

$/1000 

qal . 

0.2855 

0.2855 

0.2855 
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Che  simplex  method,  the  same  results  are  obtained  in  a  straight  forward  manner 
without  employing  such   tedious  procedures.  For  example,  if  criteria  for 
stopping  the  iteration  are  changed  from  one  set  of  value; to  another,  the  only 
thing  that  should  be  done  for  the  simplex  method  is  just  to  change  the  input 
data  for  the  criteria  to  the  corresponding  new  value,   in  contrast  to  this, 
the  parametric  search  requires  rather  labor,di/us  procedures,-  such  as  the  pre- 
paration of  a  new  set  of  nine  figures  similar  to  Figs.  8  through  16. 

However,  the  great  advantage  of  the  parametric  searcn  lies  os-  tr.e  ract 
that  it  gives  insight  into  the  MEXS  system  and  gives  detailed  information  about 
the  influences  of  the  variables  under  search,  namely  x^  and  q  /F,  or.  the  water 
production  cost  and  the  other  optimal  policies.  From  Figs.  8  through  15,  it 

is  seen  that  both  the  values  of  the  recycle  ratio  r  and  the  brine  temperature 

n 

(T.)  vary  linearly  and  slightly  with  the  value  of  q  /F.   3ut  the  water  pro- 
I  n  s 

duction-  cost  changes  considerably  with  the  value  of  q  /?,     From  Fig.  16,  it 

is  seen  again  that  the  values  of  (T^)   are  nearly  constant;  however,  they  vary 

3 

slightly  and  linearly  with  the  value  of  x.  .   On  the  other  hand,  tr.e  optimal 

recycle  ratio  r  and  the  optimal  ratio  of  the  heat  load  to  seawater  feed  q  /F 

3 
and  the  water  production  cost  vary  greatly  with  the  value  ox  x. . 

The  results  can  then  be  summarized  as  follows: 

(1)  The  optimum  temperature  of  the  flashing  brine  leaving  each  effect,  (T.)  , 
varies  only  slightly  with  the  values  of  q  /F  and  x. . 

(2)  The  optimal  recycle  ratio  in  each  effect,  r  ,  depends  or.  x^  but  varies 
only  slightly  with  q  /F. 

(3)  The  optimal  ratio  of  the  heat  load  to  seawater  feed,  q  /F,  varies  sicr.i- 

ficantly  with  x^. 

1 

(4)  The  water  production  cost  varies  greatly  with  the  values  of  x,  and  c_/F. 
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CHAPTER  6 
CONCLUSION 

In  this  study  a  detailed  analysis  of  a  MEMS  process  has  been  made 
and  quantitative  solutions  between  the  operating  variables  have  been 
obtained.   The  operating  variables  of  the  system  are  the  heat  input  to 
the  brine  heater,  the  recycle  ratio  in  each  effect,  the  brine  concen- 
tration and  temperature  leaving  each  effect,  and  the  number  of  stages 
in  each  effect. 

A  mathematical  model  of  the  MEMS  system  containing  these  operating 
variables  have  been  developed.   The  quantitative  relations  of  the  model, 
which  contain  the  operating  variables,  are  used  to  set  up  cost  equations 
which  relate  the  performance  of  the  system  to  the  unit  cost  of  product 
water.   These  equations  are  then  employed  in  the  optimization  study 
by  means  of  the  discrete  maximum  principle.   The  parametric  search 
techniques  and  simplex  method  are  used  in  conjunction  with  the  maxi- 
mum principle  to  find  the  overall  optimal  condition. 

While  the  simplex  method  gives  rise  directly  to  the  optimum  point, 
the  parametric  search  gives  detailed  information  about  the  influences 
of  the  individual  parameters  on  the  water  cost  and  the  other  operating 

variables.   It  is  obvious  from  Fig.  16  that  the  brine  temperature 

3 
changes  little  as  we  change  x.  .   On  the  other  hand  the  recycle  ratio 

rn  and  qs/F  do  change  considerably. 

The  overall  optimal  operating  conditions  are  summarized  in 

Table  5,  and  the  corresponding  capital  and  operating  costs,  the  various 

cost  items,  and  their  contributions  to  the  overall  fresh  water  cost  on 

a  percentage  basis  are  given  in  Table  6. 
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Table  5.  The  Overall  Optimal  Operating 

Conditions 

Symbols 

It 

ems 

Numerical 

values 

(1) 

Concentrations  (C.) 
v  I'D 

CF 

concentration  of  sea  water  feed 

0.035  wt. 

fraction 

(cf>i' 

brine  cone. 

entering  the  1st  -effect 

0.0397  wt, 

.  fraction 

(CA 

brine  cone. 

leaving  the  1st  effect 

0.0419  wt. 

i  fraction 

<cf)2. 

brine  cone. 

entering  the  2nd  effect 

0.0488  wt. 

fraction 

<cf)2 

brine  cone. 

leaving  the  2nd  effect 

0.0512  wt. 

fraction 

(cf)3, 

brine  cone. 

entering  the  3rd  effect 

0.0622  wt. 

fraction 

(cf)3 

brine  cone. 

leaving  the  3rd  effect 

0.065  wt. 

fraction 

(2) 

Temperature  (Tf)n 

<Tfh 

brine  temp. 

leaving  the  1st  effect 

o 
196  F 

(Tf>2 

brine  temp. 

leaving  the  2nd  effect 

148°F 

(xf)3 

brine  temp. 

leaving  the  3rd  effect 

103°F 

<Tf>0 

brine  temp. 

entering  the  1st  effect 

250°F 

(T^)3 

sea  water  temperature 

85°F 

(3) 

Recycle 

Ratios  r 
n 

rl 

recycle  ratio  in  the  1st  effect 

2.14 

r2 

recycle  ratio  in  the  2nd  effect 

2.88 

r3 

recycle  ratio  in  the  3rd  effect 

3.36 
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Numerical 

values 

(1) 

Concentrations  (C.) 
v  I'D 

CF 

concentration  of  sea  water  feed 

0.035  wt. 

fraction 

(cf>i' 

brine  cone. 

entering  the  1st  -effect 

0.0397  wt, 

.  fraction 

(CA 

brine  cone. 

leaving  the  1st  effect 

0.0419  wt. 

i  fraction 

<cf)2. 

brine  cone. 

entering  the  2nd  effect 

0.0488  wt. 

fraction 

<cf)2 

brine  cone. 

leaving  the  2nd  effect 

0.0512  wt. 

fraction 

(cf)3, 

brine  cone. 

entering  the  3rd  effect 

0.0622  wt. 

fraction 

(cf)3 

brine  cone. 

leaving  the  3rd  effect 

0.065  wt. 

fraction 

(2) 

Temperature  (Tf)n 

<Tfh 

brine  temp. 

leaving  the  1st  effect 

o 
196  F 

(Tf>2 

brine  temp. 

leaving  the  2nd  effect 

148°F 

(xf)3 

brine  temp. 

leaving  the  3rd  effect 

103°F 

<Tf>0 

brine  temp. 

entering  the  1st  effect 

250°F 

(T^)3 

sea  water  temperature 

85°F 

(3) 

Recycle 
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n 
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r2 
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Table 

5. 

The  Overall  Optimal  Operating  Condit 

ions 

(Continued) 

Symbols 

I  terns 

Nume 

rical  Values 

(4)   Flow 

Rat 

es  of  Various  Streams* 

F 

flow  rate  of  sea  water  feed 

2170 

gal./hr. 

wlt 

brine  stream  entering  the  1st  effect 

6810 

gal./hr. 

(Dx 

brine  stream  leaving  the  1st  effect 

1810 

gal ./hr . 

% 

recycle  flow  rate  in  the  1st  effect 

4640 

gal./hr. 

<L)2. 

brine  stream  entering  the  2nd  effect 

702S 

gal./hr. 

(L)2 

brine  stream  leaving  the  2nd  effect 

1485 

gal./hr. 

R2 

recycle  flow  rate  in  the  2nd  effect 

5215 

gal./hr. 

CM3, 

brine  stream  entering  the  3rd  effect 

7210 

gal./hr. 

(M3 

brine  stream  leaving  the  3rd  effect 

1170 

gal./hr. 

R3 

recycle  flow  rate  in  the  3rd  effect 

5725 

gal./hr. 

R4 

cooling  water  flow  rate 

1220 

gal./hr. 

Wl 

water  production  in  the  1st  effect 

360  gal./hr. 

w2 

water  production  in  the  2nd  effect 

325  gal./hr. 

w3 

water  production  in  the  3rd  effect 

315  gal./hr. 

(5)   Heat 

Loads  in  the  Brine  Heater* 

% 

heat  load  in  the  brine  heater 

4.87xl05Btu/hr. 

qs/F 

ratio  of  qs  and  F 

27 

qs7As 

steam  consumption  in  the  brine  heater  524 

lbs/hr. 

2W/(qs/As) 

lbs.  of  fresh  water  produced  per 

lb.  of  steam  consumed 

16 

♦Remark: 

Basis,  1000  gallon/hr.  of  fresh  water  production 
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Table  6.   Capital  and  Operating  Cost  Allocation 


Symbols      Items 


cost  ($/1000  gal)      Percentage^) 


(1)   Capital  cost 


0.10805 


37.838 


E2 


brine  heater 


?" 


heat  transfer  area 
S  E"     pump 
E(j      outshell 


(2)   Operating  cost 


0.00125 
0.08470 
0.00130 
0.02080 

0.17752 


0.439 
29.661 
0.455 
7.283 

62.162 


E7 


feed  brine 


Eo      cooling  water 
_n 


0.13131 
0.03247 
0.00568 
0.00806 


45.979 
11.370 
1.991 
2.822 


(3)   Total  water  production 
cost 


0.28557 


100.000 


Remarks : 


1.  Basis:   1000  gallon  fresh  water  production  per  hour. 

2.  Feed  brine  cost:   $0,015/1000  gal.  of  sea  water. 

3.  Cooling  water  cost:   $0,005/1000  gal.  of  sea  water.  : 
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If  a  more  realistic  water  production  cost  is  required,  we  merely 
have  to  incorporate  the  updated  cost  oata  and  cost  equation  into  the 
scheme  and  apply  the  optimization  techniques  developed  in  this  study 
for  the  MEMS  system.   Obviously  this  optimization  technique  can  be 
applied  to  the  other  sequential  multistage  and  multide-cision  processes. 
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NOMENCLATURE 

a 

«  steam  temperature,  °R 

A 

2 
<  heat  transfer  area,  ft 

3' 

«  coefficient  in  the  Clausius-Ciapeyr-on  ec.uation,  llj./ft 

3 

-  (i  ♦  fc)  B' 

=3 

.  capital  cost  per  unit  Of  heat  transfer  area  in  the  brine 
S/fi2                      • 

heater, 

C 

0 

.  unit  cooling  water  cost,  S/lb 

C  . 

co. 

0 

e 

«  unit  power  cost,  S/hp 

<CA 

■  salt  concentration  of  flashing  brir.e  at  location  n,  vt.j» 

c._. 

■  capital  cost  per  unit  of  heat  transfer  area  in  tne  conic/ 
chamber,  s/ft 

.sin; 

c.t 

-  ycB,  S/ft2 

Cj 

-  capital  cost  per  horsepower  for  the  recycle  pumps,  S/np 

c? 

«*  salt  concentration  in  the  seawater  feed,  vt  .ft 

v? 

=  heat,  capacity  of  water  per  pound,  3tu/lb,"F. 

c 

p? 

-  CJ  ^  Cc/>Jp  •  »/«-» 

st 

•  unit  steam  cost,  S/lb 

*c 

>  unit  pretreatment  cost  for  seawater  feed,  s/l'o 

2 

=  various  cost  items 

p 

»  flow  rate  of  seawater  feed,  lb/hr 

pB 

«  the  iiamiltor.ian  function  at  stage  n 
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84 

H- 

1,  H-2,  and  II-}"  the  heating  sections  in  the  first,  seco.-.i  and 

third  effects  respectively 

HE 

-l.HR-2  and  KR-3  =  the  heat  recovery  sections  ir.  the  first,  second 

and  third  effects  respectively 

h 

c 

•  unit  enthalpy  po^   pound  of  condensate,  3t'u/lb 

kf 

»  unit  enthalpy  per  pound  of  flashing  brine,  Btu/lb 

h. 

=  unit  enthalpy  per  pound  of  non-flashing  brine,  Stu/lb 

"v 

■  unit  enthalpy  per  pound  of  water  vapor,  3tu/ib 

•Jl; 

,  J„  and  J   -  the  circulation  pumps  in  the  first,  seccr.c.  and  third 

effects  respectively 

L  »  flow  rate  of  the  flashing  brine,  ib/hr 

Hx, 

,  K,  and  X.,  =  the  .nixing  points  in  the  first,  second  ana  third 

effects  respeotiveiy 

i  H-  and  S_  -  the  number  of  stages  in  the  first,  second  and  t.'.ird 

effects  respectively 

?. 

-  the  vertex  of  a  simplex  with' the  lowest  function  value,  y. 

*2 

*  the  vertex  of  a  simplex  with  function  value  y0,  yn  y,.  y, 

^    -   ^  j 

"i 

=  the  vertex  of  a  simplex  with  function  value  y,  or  the  worst 

point  of  the  simplex 

?4 

=  the  centroid  between  points  ?_  and  ?. 

> 

«  the  point  obtained  after  reflection  of  ?., 

5 

-  the  point  obtained  after  expansion  of  ?. 

0 

-   the  point  obtained  after  contraction  of  ?, 

5 
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p° 

-  vapor  pressure  of  pure  water,  lb  ./ft" 

? 

=  vapor  pressure  of  the  aqueous  solution,  lb_/ft 

.  q 

■  heat  transfer  rate,  ^tu/hr 

«s 

■  heat  input  per  unit  time  in  the  Torino  heater,  Jitu/hr 

R 

•  ideal  gas  constant,  Btu/ib,°R 

R, . 

R,  and  R   -  the  flow  rates  of  recycle  brine  strea.7.  in  the  first, 

second  and  third  effects  respectively 

H-i 

,  R-2  and  R-3  ■  the  heat  rejection  sections  in  the  first,  second 

and  third  effects  respectively 

h 

■  flow  rate  of  cooling  water,  Ib/hr 

*  <\  j 

*«  and  r.  «  the  recycle  ratios  in  the  first,  second  and  third 

effects  respectively 

s  « 

the  objective  function,  S/lOOOgal. 

'  "c 

condensing  temperature,  °R 

c»f)a  ■ 

temperature  of  the  flashing  brine  at  location  n,  *H 

T   = 

temperature  of  the  non-flashing  brine,  °R 

s 

steam  temperature  at  the  brine  heater,  °R 

(*S)  = 

f     i  '    • 

(At)  « 

c    -  j  > 

D  - 

overall  heat  transfer  coefficient,  3tu/hr.  ft   °? 

V  - 

vapor  rate,  ib/hr 

"11 

V,  and  V   ■  the  rates  of  condensate  formation  in  the  first, 

second  and  third  effects  respectively,  *o/.ir 
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3£W  «  total  wotor  production  rai.e,  lb/hr 


o 


*,  -  CO-  ,  wt.JS 


n 
x2  »  accumulated  water  cost  in  the  first  r.-th  effects,  O/lGOGgal. 

*"   ■  (?  -)   i  °H 

3    f  n 

x2  -  '^f'O  '  :esPe"4i"re  01'  the  flashing  brine  entering  the  first 
effect,  °R 

y  -  the  function  value  or  water  cost  zz   point  P 

z['   =  the  adjoint  variable  in  association  with  state  variable  x? 

IRZSK   IZTi'ERS 

eC  -  T.  -  I     ,'P 

cC  =  the  average  of tf in  the  n-th  effect 
c<  .  reflection  coefficient,  1 

£  -  vapor  pressure  depression  of  brine  streaas,  ib^/ft 
fi'»  contraction  coefficient,  C.p 
V-  expansion  coefficient,  2 

>|.  «  fractional  increase  in  pulping  power  required  due  to  friction 
'?,  *  ?ur-P  efficiency 

■5.  *  r.  t  recycle  ratio  in  the  n-z'.-    effect 
52-  (T-)   ,  °R 


<-  '•   /-n-1   «n    n-1 
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\   «  latent  haat  of  vaporization  ,  I 
X.  -   latent  heat  of  steair*  at'  2Jk.k°2 
^    -  density  of  water,  ia/cu.ft. 
\>  •  capitalization  charge 


tu/lb 

and  Vj  ps 


Btu/ls 
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CHAPTER    1 
INTRODUCTION 

The  fugacity  of  the  solvent  in  a  solution  is  always  lower  than 
that  of  the  pure  solvent  under  the  sarae  pressure  as  the  solution  (17). 
If  there  is  a  semipermeable  membrane  present  between  the  solution  and 
pure  solvent,  which  will  pass  solvent  molecules  preferentially  or 
exclusively,  there  will  be  a  net  flow  of  pure  solvent  into  the  solu- 
tion. At  an  equilibrium  state  the  fugacities  of  the  solvent  in  the 
solution  and  in  the  pure  state  become  equal,  and  therefore  there  is  no 
net  solvent  flow  into  the  solution.   This  state  of  equilibrium  can  be 
brought  about  by  raising  the  pressure  of  the  solution  to  its  osmotic 
pressure  (18).   If  the  pressure  on  the  solution  is  raised  above  the 
equilibrium  osmotic  pressure  then  the  pure  solvent  will  flow  out  of 
the  solution.   This  is  the  reverse  of  the  osmotic  process  or  the  so- 
called  reverse  osmosis  process.   This  method  requires  no  phase  change, 
therefore,  it  has  an  inherent  advantage  over  distillation  and  freezing 
desalination  processes  from  the  point  of  view  of  energy  requirement. 
However,  progress  on  reverse  osmosis  depends  greatly  on  the  availa- 
bility of  a  suitable  semipermeable  membrane  which  can  stand  up  well 
over  time  to  the  required  pressures  and  still  give  an  appreciable  flow 
of  potable  water.   It  was  not  until  recently  that  such  a  synthetic 
membrane  was  developed  and  serious  consideration  was  given  to  the  use 
of  reverse  osmosis  in  desalination. 

Much  attention  has  been  drawn  to  the  improvement  of  membrane 
fabrication  techniques  (19),  but  only  little  attention  has  been 
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exclusively,  there  will  be  a  net  flow  of  pure  solvent  into  the  solu- 
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brought  about  by  raising  the  pressure  of  the  solution  to  its  osmotic 
pressure  (18).   If  the  pressure  on  the  solution  is  raised  above  the 
equilibrium  osmotic  pressure  then  the  pure  solvent  will  flow  out  of 
the  solution.   This  is  the  reverse  of  the  osmotic  process  or  the  so- 
called  reverse  osmosis  process.   This  method  requires  no  phase  change, 
therefore,  it  has  an  inherent  advantage  over  distillation  and  freezing 
desalination  processes  from  the  point  of  view  of  energy  requirement. 
However,  progress  on  reverse  osmosis  depends  greatly  on  the  availa- 
bility of  a  suitable  semipermeable  membrane  which  can  stand  up  well 
over  time  to  the  required  pressures  and  still  give  an  appreciable  flow 
of  potable  water.   It  was  not  until  recently  that  such  a  synthetic 
membrane  was  developed  and  serious  consideration  was  given  to  the  use 
of  reverse  osmosis  in  desalination. 

Much  attention  has  been  drawn  to  the  improvement  of  membrane 
fabrication  techniques  (19),  but  only  little  attention  has  been 
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directed  to  the  system  analysis  of  the  process  itself.   Merten,  et . 
al.,  has  given  an  extensive  cost  analysis  of  a  single  stage  system  (20). 
Investigators  at  Kansas  State  University  (21,  22,  23)  have  proposed  a 
completely  mixed  model  of  the  multi-stage  sequential  system  based  on 
the  assumption  of  uniform  salt  concentration  inside  the  osmosis  unit, 
and  the  optimization  of  this  model  has  been  carried  out  by  the  same 
group  (21,  22,  23).   In  the  present  study,  the  plug  flow  model  is 
proposed  by  taking  into  account  the  concentration  change  inside  the 
tubular  osmosis  unit.   The  proposed  plug  flow  model  of  a  multi-stage 
sequential  system  is  described  in  Chapter  2.   The  quantitative  re- 
lations between  operating  variables  are  derived  in  Chapter  3.   Various 
cost  functions  of  the  system  are  established  in  Chapter  4.   The  outline 
of  the  optimization  procedure  of  the  system  is  described  in  Chapter  5. 
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CHAPTER  2 
PROCESS  DESCRIPTION 

A  simplified  flow  diagram  of  a  N-stage  reverse  osmosis  system  is 
shown  in  Fig.  1.   The  principal  components  of  the  system  are  listed 
in  Fig.  1.   Each  stage  consists  of  a  membrane  separator  unit,  MS;  a 
high  pressure  pump,  J-^  ;  and  a  recirculation  pump,  J 2-      q,  R,  and  W 
represent  respectively  the  flow  rate  of  brine  stream,  recycle  brine, 
and  water  produced.   Superscript  n  is  used  to  indicate  the  quantity 
referred  to  the  n-th  stage.   However,  the  subscripts  i  and  e  refer 
respectively  to  the  inlet  and  exit  quantities  to  the  membrane  separ- 
ator.  In  addition,  the  blowdown  turbine  at  the  end  of  the  process  is 
represented  by  J  . 

Sea  water  is  first  brought  through  a  prefilter  and  is  introduced 
into  the  first  stage  as  the  brine  stream,  q°.   It  is  then  pumped  by 
the  high  pressure  pump,  J    to  an  operating  pressure  in  excess  of  its 
osmotic  pressure  and  then  mixed  with  recycle  brine,  R    ,  at  the  mixing 
point,  M  .   The  resulting  combined  stream,  q  ,  is  carried  through  the 
membrane  separator  unit,  MS  by  means  of  the  recirculation  pump,  J2. 
The  membrane  separator  is  a  shell -and-tube  arrangement.   The  solvent 
water  of  the  brine  stream  under  a  pressure  higher  than  its  osmotic 
pressure  migrates  across  the  semipermeable  membrane  tube  to  the  lower 
pressure  shell  side  of  the  separator.   The  collected  water  product 
from  the  shell  side  of  the  first  stage,  W  ,  is  then  introduced  and 
stored  in  the  fresh-water  reservoir. 

The  exit  brine  stream  of  the  first  membrane  separator,  q  ,  is 
divided  into  two  streams.   One  stream,  q  ,  is  fed  into  the  second  stage 
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and  the  other  stream,  R  ,  is  recirculated  by  the  pump  J  ,  and  then 
mixed  with  the  sea-water  feed  at  the  mixing  point  M  .   The  operation 

of  the  subsequent  stages  is  similar  to  that  of  the  first  stage  except 

N 
at  the  last  stage,  where  the  brine  stream,  q  ,  is  allowed  to  blow 

down  through  a  recovery  turbine  before  it  is  rejected  as  waste. 

The  system  configuration  and  operation  described  here  are  not 

necessarily  optimal.   Several  versions  of  the  model  are  presented  in 

chapter  3.   The  best  model  can  only  be  decided  from  an  optimization 

study  of  each. 
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CHAPTER  3 
PROCESS  ANALYSIS 

A  fairly  general  model  of  a  sequential  reverse  osmosis  desali- 
nation system  will  be  set  up  first.   The  first  several  sections  are 
devoted  to  the  derivation  of  the  system  equations  of  "such  a  model. 
Three  simplified  versions  of  such  a  model  are  then  proposed  in  the 
last  section. 

A  schematic  representation  of  this  model  is  shown  in  Fig.  1  and 
the  n-th  stage  of  this  model  is  depicted  in  Fig.  2.   Each  stage  con- 
sists of  a  membrane  separator  unit,  a  recirculation  pump  and  a  high 
pressure  pump.   The  last  stage,  however,  includes,  in  addition,  a 
blowdown  turbine . 

The  flow  rates  of  the  brine  stream,  q,  recycle  brine,  R,  and  water 
production,  W,  and  the  superscript  and  subscript  representations  have 
been  defined  in  chapter  2.   Definitions  of  several  other  symbols  em- 
ployed in  the  derivation  are  listed  below: 

xn  =  the  mass  fraction  of  salt  in  the  brine  stream  leaving  the 

n-th  stage, 
x-  =  the  mass  fraction  of  salt  in  the  brine  stream  entering  the 

membrane  separator  of  the  n-th  stage, 

n 
xe  =  the  mass  fraction  of  salt  in  the  brine  stream  leaving  the 

membrane  separator  of  the  n-th  stage, 
q°  =  the  mass  flow  rate  of  the  sea  water  feed  (lb  /hr), 
x°  =  the  mass  fraction  of  salt  in  the  sea  water  feed, 
r  =  the  recycle  ratio  in  the  n-th  stage  defined  as  the  ratio  of 

Rn  and  q""1, 
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N  =   the  total  number  of  stages  in  the  sequence  i 

3f 

the 

process, 

W^  =  the  total  mass  flow  rate  of  the  fresh  water 

produced 

from 

the  whole  system  (lb /hr),  i.e. 

N 

wf  =  z   w" 

n=l 

p  =  operating  pressure  at  the  n-th  stage  (psi), 

AP  =  the  pressure  difference  across  the  membrane 

of 

the 

n- 

•th  stage 

(psi), 

n                                         2 

S   =  the  membrane  area  of  the  n-th  stage  (ft  ). 

3-1.   The  Fresh  Water  Production  Rate  W  and  Wf. 

The  material  balance  around  the  n-th  stage  is 

n-1    .n          n 
q     =  W   +  q            n=l,  2,  ...,N. 

(1) 

The  material  balance  for  the  process  as  a  whole  is 

q°  =  Wf  ♦  qN 

(2) 

A  salt  material  balance  for  the  first  n-  stages  gives 

0   0 
q"  =3_2_              n  =  1,  2,  ...,  N. 

x" 

(3) 

Substituting  equations  (3)  into  equation  (1).  yields 

n    0   0    l   -  l 

(4) 

q  X   '/-1    xn) 

Substituting  equation  (3)  into  equation  (2)  yields 

0      *° 
Wf  =  q  (1  -  ~) 

XT 

(5) 
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3-2.   The  Volumetric  Flux  of  Water  Through  the  Membrane,  F 

The  volumetric  flux  of  water,  F,  through  a  membrane  of  constant 
permeability  has  been  reported  (21,  22,  23)  as 

p  =      K(AP  -  12,000  x) 

s ~  (°) 


1  +  3.05  x  10° 


(Sc)1/3Da   Re7/8 


where 

ft3 

F  =  water  flux,  (^S )  . 

'  ft^-hr   ' 


K  =  the  membrane  constant,  (2         ), 

lft  -hr-psi  " 

AP  =  the  pressure  difference  across  the  membrane  (psi), 

Sc  =  Schmidt  number, 

d  =  diameter  of  the  membrane  tube  (ft), 

Re  =  Reynolds  number, 

x  =  mass  fraction  of  salt  in  the  brine  stream, 

2 
Da  =  daffusivity  of  NaCl  in  water  (cm  /sec). 

This  equation  can  be  written  as 
KAP  +  bx 


V 

1   ♦ 

r,      7/8 

Re 

whe 

re 

b 

= 

-12, 

000  K 

c 

= 

3.05 

x  105  - 

1 

;sc 

Kd 

(7) 
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3-3.   Inlet  and  Exit  Brine  Concentrations  of  MS  ,  x-^  and  xe. 

From  the  steady  state  material  balance  for  an  infinitesimal  ele- 
ment of  the  membrane  tube  as  shown  in  Fig.  3  we  obtain 

dq  =  -dW  =  -FpdS  (8) 

salt  material  balance  for  the  same  element  gives 

xq  =  ( q  +  dq )  ( x  +  dx ) 

=  sq  +  qdx  +  xdq  +  dxdq 

if  it  is  assumed  that  no  salt  passes  through  the  membrane. 

Neglecting  the  term  dxdq  yields 


*S  =  -  *1  (9) 

X       q 


Integration  of  equation  (9)  from  the  inlet  of  the  tube  to  an  arbitrary 
point  along  the  tube  gives 


4  "    x  (10) 

where  the  i  subscript  represents  the  quantity  of  the  inlet  stream. 

Substituting  equations  (8)  and  (10)  into  equation  (9)  yields 

"i  qi  %   =  FpdS  (11) 

■XT 

Substituting  equation  (7)  into  equation  (11)  yields 


(Re)7/3 


x2    (K&P  +  ox) 


9dS 
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Integrating  the  above  equation  from  x^  to  xe  and  from  0  to  S  gives 

f      £  b     ln  xe(KAP  +   bXj) 

Xiqi   (K(Re)7/eiP   "  K2AP2>    xi   (KAP  +  bxe) 

+   -A-  U-  ".!.>}  "PS  (12) 

K»P   x±    xe   J       .    _ 

In  carrying  out  this  integration,  the  values  of  AP  and  Re  were  assumed 

constant  in  the  range  of  the  whole  tube. 

After  changing  the  notation  from  x^  to  x^,  q^  to  q^,  aP  to  iPn, 

Re  to  Re  ,  and  xe  to  xe,  we  have  the  following  equation  for  the  n-th 

stage, 

r                                                                     n    n       n 
"  "    /                     ,  ,   x0(X4P   +  b  x,-  ) 
*i  1i   {     (      c      -     b    )  in  J^ i ' 

K(Ren)7/8APn   K2(iPn)2     xn(KAPn  +  bx£  ) 


_i_  (i_  -  i_)[-  PSn  (13) 

K4Pn  xn   x£ 


3-4.  Outlet  Brine  Concentrations  Between  Stages,  xn  and  xn   . 
In  Fig.  2  a  material  balance  around  the  mixing  point  Mn  is 

n         n  n-1 
q±  =  (1  +  r  )q  (14) 

salt  material  balance  around  point  Mn  is 


..n-1  „n-l 


+  xnrnqn"X   =  xn   qn  .(15) 
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Integrating  the  above  equation  from  x^  to  xe  and  from  0  to  S  gives 

f      £  b     ln  xe(KAP  +   bXj) 

Xiqi   (K(Re)7/eiP   "  K2AP2>    xi   (KAP  +  bxe) 

+   -A-  U-  ".!.>}  "PS  (12) 

K»P   x±    xe   J       .    _ 

In  carrying  out  this  integration,  the  values  of  AP  and  Re  were  assumed 

constant  in  the  range  of  the  whole  tube. 
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Substituting  equations  (14),  (16)  and  (3)  into  equation  (13)  and  noting 

that  Xg  =  x  ,  we  then  have 

(xn_1  ♦  rnxn)  It            c           b    ) 

[   K(Ren)7/8APn   K2(APn)2 

"f„  r,"/,   n,   ^,  n-1    n  h  ,)        -                   7 
x  [jUP  (1+r  )  +  b(x    +  r  x   )J   +    x          (        x   +  rn   _  x  J 

(xn-l  +  rnxn)(KAPn  +  bxn)        KAPn   x""1  +  r"x"   xn  J 

=  -£.-n"1   (£)                                          (17) 
q° 

3-5.   Reynolds  Number  Re  and  the  Recycle  Ratio  r  . 

The  cross-sectional  area  through  which  the  brine  stream  passes  at 

the  n-th  stage,  An  is  given  by 

^n  _  mnTt(d)2                                          rial 

4 

where  m  is  the  number  of  tubes  in  the  n-th  stage. 

The  fluid  velocity  inside  the  tubes  of  the  n-th  stage  is 

»"  =  2L                                             (19) 

An<? 

_   4qn"1(l  +  rn) 

m'hcd2/^ 

The  Reynolds'  number  is  defined  by 

Re"  =  du"P 
A 

=  4q""1(l  +  r") 

mn7td/i                                      (20) 

. 
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The  membrane  area  Sn  is  given .by 

Sn  =  mnd)tL 

where  L  is  the  length  of  the  separator  unit. 

Substituting  the  above  equation  and  equation  (3)  into  equation 
(20)  yields 

Ren  _   4q°x°d  (L)  (1  +  r") 

M  d  s^x11"1  (21) 

As  is  mentioned  before,  the  Reynolds  number  Re  is  assumed  to  be 
constant  inside  a  stage  in  the  derivation  of  equations  (12)  and  (13). 
From  equations  (19)  and  (20.)  one  can  see  that  Ren  is  defined  as  the 
value  of  Re  at  the  inlet  of  the  membrane  separator  in  the  n-th  stage. 
If  higher  accuracy  is  required  or  percentage  conversion  of  brine  to 
water  in  any  stage  becomes  very  high,  some  other  representation  of  Ren 
such  as  the  average  value  of  Re  between  the  inlet  and  outlet  in  the 
stage. must  be  made. 

3-6.   Energy  Requirement  for  the  High-pressure  Pump  J1   in  the  n-th 

n 
Stage,  E]_ 

n 
The  pumping  work  E^  is  primarily  used  to  increase  the  pressure 

n-1    n 
from  P    to  P  .   Since  the  velocity  difference  between  the  two 

successive  stages  is  small,  the  kinetic  energy  losses  and  friction 
losses  can  be  included  in  the  pump  efficiency.   Thus  the  power  require- 
ment for  the  high-pressure  pump  at  the  n-th  stage  can  be  written  as 


-n  _  1   *  If     P"  -  P"-1    n-1 
1  "  7m  TP        ?  q 
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where  f  „,  fjp   and  >ff   are  the  mechanical  and  pump  efficiency  and  the 
friction  loss  factor. 

Substituting  equation  (3)  into  the  above  equation  and  noting 


that 

Pn   .  p""1   =    (p»   .  pOj    .    (pn-l    _■  p0j    .  Apn      _Apn-l 
we  obtain 


n        1    +  »jf    iPn   -  4Pn-l      qo  xo 


Thus,  the  energy  requirement  for  the  high-pressure  pump  at  the 
n-th  stage  per  unit  water  production  can  be  given  in  terms  of  the  brine 
concentration  as 

3-7.   Energy  Requirement  for  the  Recirculation  Pump  j"  in  the  n-th 
Stage,  E2 


The  energy  required,  e£ ,  includes  the  energy  of  circulating 
9.   r  lbm/hr  of  the  recycle  brine  and  that  of  the  q""1  flow  work. 
The  friction  loss  comes  largely  from  the  fluid  flowing  in  the  membrane 
separator  unit.   This  lost  work  based  on  unit  time  is 

4  =  4f  i£ii  <y  qn  1  +  *t 

2ac    d     l     la    1p  (24) 

where  f  is  the  friction  factor. 
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From  rearranging  equation  (19),  the  flow  rate  of  the  brine  stream 
entering  the  membrane  separator  in  the  n-th  stage  can  be  written  as 
follows 

qn  =  ran  itdf,  u"f  (25) 

4 

For  turbulent  flow  the  friction  factor  can  be  approximated  by 

f  =  0.046 


(Ren)0-2 


Substituting  the  above  equation  and  equation  (25)  into  equation  (24) 
yields 

1  +  V 
E"  =  0.023  |£  _?_  l_M_)3      (Ren)2-8  Sn  (26) 

where  S  =  mnL*d,  the  membrane  area  in  the  n-th  stage. 

The  energy  requirement  per  unit  water  production  is 


*2  =  0.023  Lli£  ±  ^   (^)2-8  -J- 


V     IP      9c        d?  q°(l  .  £) 


(27) 


3-8.   Energy  Recovery  at  the  Reject  Brine  Turbine,  E, 

The  energy  recovery  from  depressurizing  the  high-pressure  brine 
stream  from  r     to  P   (discharge  pressure)  is  given  by 

E3  =  Mp(1-V^-T^   * 
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Noting  that  P^  -  P°  =^P^  and  substituting  equation  (3)  into  the 
above  equation,  we  obtain 


E3  -  W1   -%)     ^f  ^f  <28> 

1         V 


The  energy  recovery  per  unit  water  production  can  be  written  as 

E3  =  *  V        (1  "  fr>   ^     x°  (29) 


3-9.   Simplified  Models 

Three  simplified  models  A,  B,  and  C  are  presented  here.   Model  A 
is  essentially  similar  to  the  general  model  except  that  it  has  the 
same  membrane  area  in  each  stage  (i.e.  s"  =  S).   Model  B  is  a  sequen- 
tial system  without  recirculation  pump  in  each  stage  (i.e.  r   =  0)  and 
is  depicted  in  Fig.  4.   Model  C  is  a  sequential  system  with  only  one 
pump  in  the  first  stage  and  is  shown  in  Fig.  5.- 

The  model  A  has  been  suggested  by  the  fact  that  it  is  often 
economical  to  use  an  identical  unit  at  each  stage  for  a  multistage 
system. 

For  a  sequential  multi-stage  system,  if  the  cost  function  is  in 
the  linear  form,  the  system  with  recycle  operation  is  often  an  optimal 
configuration  (23,  28).   However,  if  the  cost  function  is  in  the  non- 
linear form,   this  may  not  be  true.  This  has  given  rise  to  model  B. 
It  also  appears  that  it  may  not  be  necessary  to  use  a  high-pressure 
pump  in  each  stage  but  just  to  let  the  pressure  of  the  brine  stream 
decrease  successively  in  each  stage.   Therefore,  the  model  C,  which 
is  simpler  than  the  model  B  is  proposed. 
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Reasons  for  proposing  different  models  may  not  be  entirely  valid. 
They  can  only  be  verified  through  an  optimization  study  of  each  model 
using  both  linear  and  non-linear  representations  of  the  cost  function. 

The  energy  requirements  for  each  model  are  derived  below.  The 

j  .  n    n 

second  subscripts  a,  b,  and  c  in  the  various  energy  terms,  E-  ,  E-h, 

n 
and  E^c  are  used  to  represent  the  various  energy  terms  of  the  models 

A,  B,  and  C  respectively.   There  is  no  such  additional  subscript 

attached  to  the  general  model. 

Model  A. 

The  basic  assumption  of  this  model  is  to  use  the  same  number  of 
tubes  in  each,  stage.   Since 

m  =m,    n  =  l,2,  N, 

Sn  =   mnditJL  =  md/tL  =  S    ,         n   =   1,    2,    N  (30) 

n   _  mn;t(d)2        mtt(d)2  ,  A,    n   =   1,    2 N  (31) 

4  4 

After  changing  Sn  to  S,  equations  (17),  (21),  (26),  and  (27) 
become 

(x""1    ♦   rnxn)    /  (   £ "  5 


(x""1    +   rV)    /  (   £_ 

I        K(Ren)7/8apn 


K~(AP")~ 

to   xn(K^Pn(l+rn)    +  b(xn-X    +   rnxn))     +    _J_       1    +   r"  i_    1 

('xn_1      +   rnxn)(KAPn  >  bx")  "    +    K4Pn    ^   x""1    +   rnxn   "  x°   '\ 

-  4     ^    ("%)  (33) 

X  q 
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(x""1    +   rV)    /  (   £_ 

I        K(Ren)7/8apn 


K~(AP")~ 

to   xn(K^Pn(l+rn)    +  b(xn-X    +   rnxn))     +    _J_       1    +   r"  i_    1 

('xn_1      +   rnxn)(KAPn  >  bx")  "    +    K4Pn    ^   x""1    +   rnxn   "  x°   '\ 

-  4     ^    ("%)  (33) 

X  q 
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and 

D„n   _    4q°x°d       (L)       (1    +    rnl 
M              d          s   xn-l 

(33) 

and 

E"      =    0.023    Z    +^f        ?        (^)3     (R-")2-8      S      ■ 
2a                      V?P          Sc        df 

(34) 

and 

«L    -«.«»      *    +    fc      ?        f^l*      fSe"!2-8 

(35) 

Wf                          ft$           9c        d?                                q°(l    -  >£) 

xN 

respectively. 

Equations    (22),    (23),    (28),    and    (29)    are   still   valid  for 

this 

model.      Therefore, 

n     _   ^n     _   1   +  ?f     iP"   .  APn-l     qox.o 

(36) 

*          MP                     ?                x""1 

Zia  =  fl   =   l   *  If     APn   _  apn-l                 xo 

(37) 

Wf         Wf       fa  ?p                   o               x""^!    -  x£   , 
5                                     xN 

E3a  =   E3  =    Uh    (1    -ft)      £^     ^ 

1              ?          xN 

(33) 

E3a       E3                                         ipN         x° 
Wf           Wf           Vp    U        V         c.        x-N-x0 

<39) 
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Model  B 

As  is  shown  in  Fig.  4,  only  a  high-pressure  pump  is  used  in  each 
stage  instead  of  a  high  pressure  pump  and  a  recirculation  pump.   The 
velocity  of  the  brine  stream  in  each  stage  may  be  adjusted  by  using  a 
different  membrane  area  in  each  stage. 

After  dropping  the  recycle  ratio  rn  in  equations  (17)  and  (21) 
we  obtain 


b     ln  xn(KAPn  +  bx""1) 


K(Ren)7/8aPn   K2(APn)2    xn_1(KAPn  +  bxn) 


+  K^Pn(Xn-1"  x"  J 


-%  I ■  TT  )  <4°> 

X     q 


and 


-  n  _  4q°x°d   /L%    1 

Re        >T"   {d]  s^=^  <41> 

Equations  (28)  and  (29)  are  still  correct  for  this  model.   Therefore, 
we  have 

E3b  =  E3=WP  <X  "'&)  IT^f  («) 

J    x 

f3b  =  E3   B  »  M      (1  -  /7.)  £E_  * (43v 


For  this  model,  the  energy  requirement  for  the  high  pressure  pump 

n 
In  the  n-th  stage  E-^   includes  not  only  the  energy  used  to  increase 

the  pressure  from  P     to  P  but  also  that  of  the  pumping  head  to 
overcome  friction  in  the  n-th  stage.   Therefore,  the  following  expres- 
sions are  adequate 
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E"b_ 
wf 

„n    „n 
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En 
2b  = 

Wf 

0 

(47) 

Model 

C 

Model  C  is 

shown  in  Fig 

.  5. 

In 

this  model  only  one 

high  pressure 

pump 

is  used  in 

the  first  stage, 

i.e. 

,  no  pumps  are  used 

in  the  re- 

maining  stages. 

Since  pressure  c 

hanges  along  the  tube  in  each  stage, 

an  exact 

solution  involves  a 

complicated  intergration.   An  approximate 

solut 

ion 

can  be 

obtained  if  we  assume 

constant  pressure  inside  each 

stage 

but  change 

s  from  stage 

to  s 

tage 

abruptly.   Therefore  equations 

(40) 

and 

(41)  oi 

model  B  are 

adequate 

here,  but  the  pressures  between 

the  n 

-th 

and  the 

(n-l)-th  stage  can  be  related  by  the  xo] 

.lowing  relation: 

(26) 
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CHAPTER  IV 
ECONOMIC. ANALYSIS 

Cost  of  the  plant  may  be  divided  into  the  two  major  parts:   the 

capital  cost  and  the  operating  cost.   The  capital  cost  consists  of 

X 
three  items: 

(a)  Pump  cost, 

(b)  Turbine  cost, 

(c)  Membrane  area  cost. 

The  operating  cost  includes  four  items: 

(a)  Power  cost  for  the  high  pressure  pump, 

(b)  Power  cost  for  the  circulation  pump, 

(c)  Energy  recovery  from  the  reject  turbine, 

(d)  Feed  brine  cost. 

Other  costs  such  as  labor  cost,  insurance  cost,  etc.,  are  not 
considered  here  as  they  have  little  effect  on  the  water  cost  when  the 
operating  conditions  are  changed. 

The  symbols  which  represent  the  cost  items  mentioned  above  are 
listed  below.   The  first  subscripts  are  referred  to  the  various  cost 
items  and  the  second  the  various  models. 

C.  ,  C    Cih''^l  "c^a  capital  costs  of  the  pumps  for  the  general. 

.model,  and  models  A,  B,  and  C,  respectively. 

^2'  ^2a*  ("2b*  ^2-  =  tne  caPi::a^  costs  of  the  turbine  for  t'r.^   genera 

model,  and  models  A,  3,  and  C,  respectively. 
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C-i>  C,  ,  C   ,  C,  =  the  feed  brine  costs  for  the  general  model,  and 
/    /a   Jo        /c  °  ' 

models  A,  E,  and  Cj  respectively. 
C  ,  C   ,  C  ,  ,  C   -  the  total  water  costs  per  unit  mass  of  production 
for  the  general  modal,  and  models  A,  3,  and  C, 
respectively. 

4-1.   The  Capital  Cost 

The  annual  capitalization  charge  for  the  equipment  items  is  taken 
to  be  0.074  of  the  initial  cost  per  year,  as  recommended  in  z'r^   Office 
of  Saline  Water  Report  (12) .   An  assumption  of  a  load  factor  of  330-on- 
stream  day  per  year  gives  a  capitalization  charge,  y,  of  9.4xl0~  of  the 
initial  cost  per  hour  on  stream. 

The  power  rule  relating  the  capital  cost  and  the  equipment  capacity 
as  suggested  by  Chilton  (27)  is 

C  -  k  Ta  (53) 
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where 

C  =  the  capital  cost  of  the  equipment 
X  =  the  capacity  of  the  equipment 
k  =  a  proportionality  constant 
a  =  a  positive  constant  less  than  1.0. 

This  rule  is  used  where  the  capital  costs  of  the  equipment  are 
concerned. 

(a)   Pump  and  Turbine  Cost,  C   and  C, 

The  proportionality  constants  in  the  power  rule  for  the  high 
pressure  pump,  recirculation  pump,  and  turbine  are  represented  re- 
spectively by  k1,  k2,  and  kt.   According  to  the  power  rule  mentioned 
above,  the  pump  and  the  turbine  costs  are  given  respectively  by 


N    f  k  (E")  l      *   k  (E")  2 
cx  =  ty   Z.  — 2 — § — 


(59) 


and 


(60) 


where  a-^ ,  a2,  and  a3  are  the  power  rule  coefficients  for  the  high 
pressure  pump,  the  circulation  pump,  and  turbine,  respectively. 

(B)   Membrane  Separator  Cost,  C, 


The  mass  of  the  membrane  separator  for  the  n-th  stage  is  given 
by  (20) 
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,  ,n    <?  d      . n  ,  n  ,,  ,     0.54   0.189  , 

Ws=^_  5  AP  (1.62  +  t7^ttW~     ) 


AP 
and 


w"  =   Jm_     S"  Ap"      (1.62  +  0.54  ,  0.1S9   /   °"m  j  (62) 

%.         q°(l  -  £)  L/°     L/D  V  AP" 

N  .  . 

x 

where 

n 

W  =  the  mass  of  the  membrane  separator  for  the  n-th  stage  (lb  ), 
s  *  m 

<?    =   the  density  of  the  material  of  the  construction  (lb  /ft3), 
'in  *   m 

<Tm  =  the  allowable  stress  of  the  material  of  construction  (psi), 
L/D  =  the  overall  length-to-diameter  ratio  of  the  membrane 
separator 

The  proportionality  constant  in  the  power  rule  for  membrane 
separator  is  represented  by  kg.   The  membrane  separator  cost,  C,,  can 
then  be  written  as 

N    (<)    ^ 
C3  =  ksS      w—         ($/lbm)  (63) 

n=l     ™f 

where  a.  is  the  power  rule  constant  for  the  membrane  separator. 

Since  we  have  assumed  certain  constant  values  for  L,  d,  and  ("g) ,  this 

gives  rise  to  an  inequality  constraint,  (— )D^i>Sn  that  must  be  satis- 

d 
n 
fied  in  the  selection  of  S  . 

4-2.   The  Operating  Cost 

The  unit  electrical  power  cost  is  represented  by  Ce,  $/psi-ft  . 
The  power  cost  for  the  high  pressure  pump  C4,  for  the  circulation  pump 
C  ,  and  the  energy  cost  recovered  from  the  reject  turbine  C.  are  given 
by 
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E 
N      h 


W       ($/lbm)  (64) 


and 


and 


N   Pn 
C5  =  Ce  X   ^2        ($/lbm)  (65) 

n=l   W. 
C6  =  Ce   -1  ($/lbm)  (66) 

wf 

respectively. 

The  cost  of  the  brine  feed  per  unit  water  production  can  be  given 
by 

C7  =  CF   c£     ($/lbm)  (67) 


Substituting  equation  (5)  into  the  above  equation  yields 

N 
c7  =  CF      x  (&8) 


The  total  water  cost  per  unit  water  production  is  the  sura  of  the 
various  cost  items,  i.e., 


C-t  =  Ci  +  C2  +  C3  +   C4  +  C5  -  C6  +  C7 
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C-t  =  Ci  +  C2  +  C3  +   C4  +  C5  -  C6  +  C7 
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4-4.   Water  Costs  for  Models  A,  B,  and  C 

The  derivation  of  the  cost  equations  for  models  A,  B,  and  C  is 
similar  to  that  for  the  general  model.   The  cost  equations  for  the 
various  models  are  listed  below: 


(A)  Model  A  a  a 

N  f  k  (En  )  1    +  k  (En  )  2 

r         ■    hi  c  I       1   la       2   2a 

cla  ■  f  £ .  (70) 


C. 


n=l 


2a  "  fkt   Hj—  (71) 


n   "4 
N      (W   ) 

C3a  "  f*s  =F     JJ—  (72) 

n=l        I 


where 


C"   ^  -"   ^*^  ^j^>  (-) 

C4a  =  Ce   X      il2.  (74) 

n=l    Wf 

N     E" 
C5a  -  C   Z     ^a  (75) 

n=l    W^ 


C6a  =  Ce   fsa  (76) 


C?a  =  °F   ^~^?  (?7) 


cta  "  cla  +  c2a  +  C3a  +  c4a  +  c5a  "  c6a  +  C7a  (78) 
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(B)  Model  B 

r              1 1   r 

(75) 

Clb   K.k  C 

n=l    W, 

a3 

c    -   V 

(£0) 

C2b   *t      „_ 

N   (V/*)*4 

C-,  ■  tlik  £ 

3b     s  . 

n=l    W 

'  (81) 

N    S?. 

r     r         *° 

(82) 

4b    e   ,  -  W . 
n=l    f 

■    C5b  "  ° 

(£3) 

s„ 

C     " 

6b   e   wf 

(84) 

s 

r       _        X 

(85) 

C7b   CF    v 

><      0 
X   -  X 

cb    lb    2b    3b    4o 

-  c,,  +  C-, 
6b    7b 

(86) 

(C)   Model  C 

. 

1  al 

(87) 

Clc  "  *  V 

W- 

i 
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(C)   Model  C 

. 

1  al 

(87) 
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(E,J 


C2c  *  *kt 


(88) 


N    .  (Wn)  4 


C3c  -  iK     £ 


n-1     W. 


(£9) 


c;  .=  c 


lc 


(SO) 


C.  =  0 
5c 


(SI) 


6c    e 


"3c 


(52) 


C,  -  C„ 

7c    F 


N    o 

X   -  X 


(S3) 


tc    lc    2c    3c    4c    oc    7c 


(94) 
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CHAPTER  5 
OPTIMIZATION 

Equations  (17)  and  (21)  show  the  relations  among  the  various 
quantities  of  the  n-th  stage  such  as  the  brine  concentration  x  , 
the  recycle  ratio  r  ,  the  pressure  difference  across  the  membrane 
AP  ,  the  Reynolds  number  Re",  and  the  membrane  area  Sn.    The  water 
production  cost  is  a  function  of  these  5N  variables  and  is  given  by 
equation  (69).   However,  the  2N  relationships  given  by  equations  (17) 

and  (21)  reduce  the  number  of  independent  variables  from  5N  to  3N. 

N 
Since  the  brine  concentration  leaving  the  last  stage,  x  ,  must  be  pre- 
fixed in  order  to  calculate  the  energy  requirements  in  each  stage, 
the  total  number  of  independent  variables  becomes  3N-1. 

A  discrete  version  of  the  maximum  principle  is  powerful  for  search- 
ing the  optimum  condition  of  a  multi-stage  multi-decision  process. 
For  the  process  in  hand,  which  is  an  N-stage  3-decision  process,  the 
performance  equations  are  summarized  in  section  5-1 .   The  derivatives 
of  the  state  variables  are  listed  in  section  5-2.   The  adjoint  vari- 
ables and  derivatives  of  the  Hamiltonian  functions  are  determined  re- 
spectively in  sections  5-3  and  5-4.   The  computing  procedures  are  de- 
scribed in  section  5-5. 

Similarly,  the  same  procedures  can  be  applied  to  the  simplified 
models. 

5-1.   Performance  Equations 

Equation  (17)  can  be  rewritten  with  the  aid  of  equation  (5)  as 
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,    n-1       n    n.f.                0                    .              b                         X^K^f  l+XD)    +    b|MrV )) 

1                        |l(Re)7/8Apn        k^p")2'    *"         (x""1    -    rnxn)(KAPn   ♦  bx") 

+     1             1  ♦  r»              l   ' 

?                  ° 

Wfx°       (           ^n)    x          S         (95) 

KAP      lxn-l    +   rnxn        ^ 

a  =  1 ,    . . . ,  n 

Rearrangement  of  equation    (21)    gives 

fi{X   -  x^x*1"1    Ren  Sn      ' 

T  "                                       X'                                                                        -     '                                         r,      -      1                                IT               f  1  fi  1 

4  x°   d(L/d)    Wf 

Substituting   the   various   equations    into   equation    (69)    yields 

Ct,±       l{kl(1+;?f    .AP--AP"- 

0           i'   ai 
(l-Xljx1^1 

X 

•  ■*-  f.  fc  i  *• <-v" s")  1 

+  w7  M  c1  -  v                 0    J 

N 

X 

♦  Jk    1      f5£     s"   iPn    (i.62   +   0.54   +   0.18S    f^.  }]  4 
Wf      n=1    U„                       (                   L/D           L/D    -JipnJ 

N      c    1    t1!,      ,nn         .„n-l                  o                                       l    +   >;      o              , 
,    -       j-     [              If     AP      -  AP                      x                    +nnoi                'f    V     (/«)3 

6   n=l    L    V(p                ?                   xn-l(i    .  X°                             ?m  ft>   3C      a? 

x° 

(Ren)2.8gJ        .Ce?m^p      (1    .^, 

APN           x° 

9       N.  _      o      ^      *"  ' 

_>        X              X 

x° 

N           o 

X        "    X 

(97) 
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The  3N  decision  variables  and  2N  state  variables  are  defined  as 
follows: 

gn  =  Ren  n  =  1,  ..  .  ,  N                   (98) 

g"  =APn  n  =  1,  ...,  N  .                 (99) 


|£  =  Sn  n  =  1,  ...,  N  (100) 


x"  =  x"  n  =  1,  ...,  N  (101) 

N 


(102) 


After    such   transformations   equations    (95)    and    (96)    become 


r  xnf  Kfln(l+rn)    +  bCx"-1*!1^)] 
n-1           n  n     f                 c                            b         „        lL     V  ;  K    2  2'-> 

(x,  +    r    X,  )   (   ( ; -   lYr,  

1  (   \f^n\7/s  nn        „2,nr\2  <    n-1    .      n      nWIr.n       .    .    n, 

\  K(fix)       82     K  (82)  vx!         r    xi)(K62  iJ 

1  1    +   rn  I1)  '  x?        n-l  n 

+  (  "  — )     =   B    (1  -  -J.)    x"  g"         <103> 


K$l        xj"1    ♦   rnx£         xl      ]  xN 

where 


B      =  — I —  n   =    1,     . ..,    N 


1     x?wf 


and 


rn   =   B2    (1    -  — J   ^     gn      xn- 


1      n   =   1,    ...,    N  (104) 
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where 


B,   = 


M 


4   x°    d(L/D)Wf 


x2   is   defined  as   follows : 


cn       Qn-1      a.    • 

n  n-1  (  °2    ~  y2  I  x  r  0   0  1 

*2    =   *2         +   33  I  34      xn-l(1.xoJ  *  B^B-t    »)*•»     J  "  J 


xN 
1 


[e"  f>2 


(B?    ♦   B^^f4 


fi"  -  fi"'1 

°*       vi  n  2.8     .  n 

+  B9     17TT To   +  B10(Si)  S3 


(105) 


n   =    1,    . . .,    N-1 


»N       .N-1 


N           N-1         „     [  „         u2        °2 
X-,    =   x"         +   B„      B„      

2  2  3  I      4        n 


,al 


*1     (1    -  2) 

N 
1 


-)    ^hii!)m»! 


B6[fia    6*    (b7*b8(B^J 


i.   *4 


Bio(6^)2-8  6^ 


N- 
xl 


+   &1A      B 


[-  ArJ 


N  ■ 


14      Bll 


1  1 


-   Bi2 


+   Cc- 


*-■;     *?-*? 


(106) 
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where 


\i>k 
B,   =  i-± 

3      wt- 

1    ♦'*. 


f        xl    wf 


B5    =    0.023   ———     j~      if-) 


t*. 


B8    =   0.189 


Mp      f 

B10  -   Ce  BS/Wf 


x      W 

Bll  Ma  %  (»  "  ^>     -^ 
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=  <j^2 

B14 
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Equations  (105)  and  (106)  contain  the  term  Q2~  >  that  is,  it  has 

.  .    .  n 

decision  in  memory.   To  avoid  this,  x,  is  introduced  as  (11) 

n   .  n 
x3=02.'  n  =  l,...,N  (107) 


Then  equations  (105)  and  (106)  can  be  rewritten  as 

k      xi   (*  -  !ir 

xN 

1 

+  s6  Co?  e^  <b7  *  H8(e^)-)]a4  ♦  b9  ^  "X3Ko 

xl      f1    "  7> 
xl 


♦   B1Q(6")2-8  63  n   =   1,     ...,    N  (108) 


n     N-i        i      e2  -  x3-1   >*i         r     rn  2*8  flNla2 

X2   =   X2         +  B3    [  B4      N2  J  )         ♦   B13  (85(61,  8  3J 

+  b6  (e»  e»  (b7  *  vis)*))-*  *  b9  e^ ;  r  o 

*? 

N  flN 

N  2.8     N  f  82        \a3 


N   2.8      N  r  »2         \a3  ' 

*  b10  (Si)     e  3 +  Bi4    Bn  s — -J     -  \2  ir-z- 

xl  -  V  *i-*x\ 

1 

+    CB-    ^     -    X° 

F     1  1  (109) 
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Now,  the  optimization  problem  may  be  formulated  as  this: 

Find  a  set 

of  decision  variables  Q^,   S2>    and  S3  (n  =  1i  2>  •••  N) 

to  minimize  the 

water  cost  x  with  xv  prefixed. 

5-2.   Derivatives  of  State  Variables 

For  convenience  the  symbols  g  and  h?  are  used  to  represent  the 

various  combinations  of  the  state  variables  x?,  decision  variables 

@",  and  constants  Bn  as  defined  before.   The  two  symbols  are  listed 

respectively  in 

Tables  1  and  2. 

(1)   x? 

2x1 

n 
g 

43           n  =  1,  2,  ...,  N-l             (HO) 

<  = 

gn 

41 

N 
jx1 

N 

f*a                                                    (in) 

gN 
42 

!i. 

n 
2M.                          n  =  1,  2,  ...,  N-l             (112) 

»es 

y41 

>xx   = 

< 

N 
°44                                            (113) 

N 
942 

K_ 

n 
945           n  =  1,  2,  ...,  N-l             (114) 

^3 

gn 

41 

2  x 

N 
945                                                (115) 

J42 
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f*a                                                    (in) 

gN 
42 

!i. 

n 
2M.                          n  =  1,  2,  ...,  N-l             (112) 

»es 

y41 

>xx   = 

< 

N 
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N 
942 

K_ 

n 
945           n  =  1,  2,  ...,  N-l             (114) 

^3 

gn 

41 

2  x 

N 
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J42 
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n          n 
?Xl   =   g46                           n   =   1,    2,     .  .  . 

->nn          n 

,    N-1 

(116) 

N               N 
2  *l             946 

■ 

(117) 

1              a42 

?X1        =   0                           n   =   1,    2,     ... 

,    N 

(118) 

-,   n"l 
7X2 

9x" 

7    1        =   0                           n  =    1,    2,     ... 

„    n-1 
2x2 

,    N 

(119) 

(2) 

n 
X2        . 

n 

p            n          n 
— -  =      h7   +  h8               n   =   1 ,    2 ,    . . . 

,    N-1 

(120) 

__2        -        h?        +        hg         +        h16 

5  6? 

(121) 

-    n 

"^=h17      +h?8      +19      n   =   l>    2,    ■ 
,y2 

. ..,    N-1 
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3  xN           N           N              N              N             N 
S.  "   h17+  h18    +   h19   +   h20  +   h21 

(123) 
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n          n 
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,    N-1 
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N               N 
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■ 
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,    N 
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9x" 
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„    n-1 
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,    N 
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n 
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n 

,  n           ,  n           ,  n 
"22   +   «23    +      24 

n   =   1, 

2,     . 

..,    N-1 

(124) 

■>     N 
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>6* 

N            N             N             N 
n22  +  n23   +   h24   +  h25 
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3x2 

"    h26 

n  =   1, 
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..,    N-1 
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n 
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n 

0 
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,.,    N 

(130) 

^x3    = 

a82 

1 

n   =   1, 
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...    N 
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n 

,  n           ,  n           ,  n 
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2,     . 
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>6* 

N            N             N             N 
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?x. 
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3x, 
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^x. 
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Table   1 .      Symbol    Representation  of  g 

n          n-1            n  n 
91   =  xl          +   r  xl 

gn  - 

K6n2  (6n)7/8 

»  -        b 

93        K2(@2V 

g"  =  x^[K82    (1   +   rn)   +  b9l] 

%  =  %  +  bx? 

g"   =       1          (1    *   rn         1  ) 
K  fln             a"                n 

K  82          gx          xx 

o 
n                ,           xl 
97  -  H    (I   -  ^) 

1 

x° 

98    =  B2    (1    "  4   ) 
Xl 

n 

n.n            n                 94            n 
99  =  (g2     -  93)     m            +  g6 

is 

- 

n     _      n.n   n-1 
910   "    98S3X1 
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- 

n     _      n.n   n-1 
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n 

Hi  ' 

n 

912  = 

n 

913  " 

n 

914  " 

9n     = 
15 

g"      = 
16 

n 

917  " 

n 

918  = 

n 

g19  = 

n 
920 

Table   1.      Symbol    Repress 

n   n  n 
Xl99910 
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■ntation   of  gn    (Continued) 

132 

n,    n          n,    n 
-      91<92    "   93>xl        "   n 
n                    95910 

■  94    ■ 

,    n            n.      n  n 

(g2     -  g3)  xl9l0 

n 
910 
...          /i           ni      n          n 

rngn 

n    ,    n          n             n 
91    <92   -  93)      ,94   ^       an, 

n                      n             1       ' 
g4                 *! 

-(rf  -  g") 

J         ,    n   n        .     n. 

(y-i      +  by    ) 
95 

n     - 
91             1             (1    +"t°)      r11  ' 

K8^  (xn)2           (gn)2 

n,    n            n                     n 
yx(g2  -  2g,    )             g4 

—     ln  

fin                    gngn 

b2                                   ly5 

n 
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n 

913  " 
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15 

g"      = 
16 

n 
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n 
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Table  1.   Symbol  Representation  of  gn  (Continued) 


n,  n    n. 
n      9, (9-  -  g,)      n       n 

gn  =  .  _l_2 3_  K  "  (1  +  rn) 

21         n        * 


K  gn(gn  -  gn) 
n       ~V  J2   u3; 
9c 


"22         g. 


n  n 
9,9A 


n     "lJ6 
923  =  "Tn" 


n     n  n  n-1 

924  ■  98  ex  xx 


9"  =  gnxn_1  -  x"gngn 

25    y7  1      1  9J24 


n   n-     n 
n    9i(g2   "  93)   n  n  n 
926  =  n xl92495 


n   _  ,  n     n.   n  n 

927  -  (g2  -  g3)  xl9 


24 


n 
928 

= 

n 
924 

gn 
y29 

= 

93M3 

n  n  I.   * 


(1  +  r  ) 


g"   =  1  +  xn  gn 
a30         1  y29 


gn  =  gnfln  -  gngn   +  (gn      -  gn)  gn 
U31    J7H3    y9J30    ^2  J3J  93 


"3'  "30 
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Table  1.   Symbol  Representation  of  gn  (Continued) 
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924 

gn 
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93M3 
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Table  1.   Symbol  Representation  of  Qn    (Continued) 


n ,  n    n_   n 
J32 


n_  .  gi(g2  -  03)  *i  (Ken  gng  +  bgnQ) 


9, 


n  n    ,,  ,   n.  n 
32991  "  <X  +  r  >930 


33        „»n  n 


„n   _  .n  n-1  x 
934  "  fl3  Xl 


n  2 
(x,  ) 


935  "   B29l  9 


34 


936  "  Bl934 


gn  =  gngn  xn  -  gn 
37    9  35  1     36 


n   n    n 

938  "  -7T^-    Xl  95  93S 


n   _    ,  n    n.   n  n 

g39  -  -  (92  -  g3)  xl935 


n   n         n,   n 
gn   .   935  8j  -  <X  +  r  j  xl 
40  ,san    n 

Ke2gx 


g    =  gn  +  g"  +  gn  +  qn 
a41    y16   y17    yl8    J19 
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Table  1.   Symbol  Representation  of  Qn    (Continued) 
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n  n    ,,  ,   n.  n 
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33        „»n  n 
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Table  1.   Symbol  Representation  of  gn  (Continued) 


«n_n.n.n     n     n 
g   ■  g   +  g   +g   +o   +  q 
42    y41    J37    y38    J39    J40 


n     n     n     n     n     n 

343  =  Sj!  +  g12  +  g13  +  g14  *  g15 


gn   =  gn   +  gn   +  gn   +  gn 
44     20     21     22     23 


n     n     n     n     n 
945  "  925  +'  926  +  927  +  928 


gn  =  gn  ♦  gn  +  gn 

46     31    y32    y33 
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Table  1.   Symbol  Representation  of  gn  (Continued) 
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Table  2.   Symbol  Representation  of  h 

fl"  -fl""1 
^n   02   °3 

1   x^d-l) 
N 
Xl 

*%  =  (e^)2,8e^ 

4  =63  air* 

h4  =  Bs +  B7  <©2>3a 

*.  •; 

N     0 
Xl  "  xl 

N 
hn  =     xl 

6     N     0 
xl  "  xl 

a                                  n   n  1  H     n  a2_1 
h?  =  2.8a2B13Bs93  (9")  -   (B5hg) 

b»-8.a^{ij)*.ae; 

0 

(x°  -  x°)2 

„n    n-1 

h10"   2     3    .  hn 
iu     n-1       9 

xl 
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Table  2.   Symbol  Kepresentation  of  h   (Continued) 

hll  =  aiB3  (Vl^^loV 

h12  =  VlO 

h^  =  *3B14B11  t8!!^3  ^9 
o 

xl 

1 

n«4  .  C^   -  *l2f  3 

xl 

.  n   _   n     n    ,  n     ,  n 
h15  "  hll  *  h12  +  h13  +  h14 

n     n   ->xn 
h16  =  h15  !S 

w£ 

a  -1 
_  „  ,  n  ,„  ,  n .  l 
n   .  a1B3B4h6  (B^) 

hl7          n-X 
*1 

n  n  4    n      n 
„n   _  a4B6(h3h4)     (h4  +  B7h3) 

2(6^ 

n-l 
xl 
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Table  2.   Symbol  Kepresentation  of  h   (Continued) 
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*1 
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2(6^ 

n-l 
xl 
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Table  2.   Symbol  Representation  of  hn  (Continued) 


a  -1 


h2i  =-sri-T(a3Bi4Bii  <Buh5)"3 


hn   =  a  B    (B  h")*2"   B   (fln)2-8 
22     2  13    5  2       5  ^el' 


h23  "  *4B6  ^V/^    ^H 


n  =  n    /nn,2.8 


h"  =  B    (nn) 
24    10  "V 


n     n  3x1 
h25  "  his  -r- 


*26  =  -  ^  [*9      ♦  BlB3B4  (B^)*1"1] 


n  =  .  n   3x1 


h27  =  h15  T^T 
?X1 


h28 


n-1 


"29  "  "28  ]a9   "  "l°3B4  lV" 


a,-l 
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Table  2.   Symbol  Representation  of  hn  (Continued) 
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a,-l 
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5-3.   Adjoint  Variables 

N 


(a)  Adjoint  Variables 
Since  xvj  is  the  total  cost  function,  we  have 


and  we  can  write  (11) 


N 
However,  since  x^    is  prefixed, 


H* 


Then  tr  becomes 


1  1     2 


Cdfferentiating  HN  with  respect  to  63  yields 


N      N 
gHN  =  ZJJ  »£  +  3x'^ 


Setting  =  0  yields 

?6* 


. 

>x» 

4  = 

,•» 

N 

>e» 


(136) 


****.  -5-1.      -5-0.  (137) 


N 

*!  *  ex.  (13a) 


(139) 
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(b)   Adjoint  Variables  z. 

From  the  definition  of  adjoint  variables  (11)  and  the  known 

N 
values  of  zij  and  the  derivative  of  state  variables  in  section  5-2,  we 

obtain  the  following  expression  for  zn 


1 

-r1 

-r1 

z'n-1 

1 

-i 

,   n-1 
pX2 

+ 

n 
z2 

-2 

2 

-      n-1 
?X2 

n-1 

Z' 

-■5 

>"l 

+ 

n 
z2 

n 
3x2        . 

3 

>*r 

-    n-1 
*X3 

n  =  1,  2,  ...,  N       (140) 


n  =  1,  2,  ...,  N     (141) 


n  3x 
'-2 


n-1 


n  =  1,  2,  ... ,  N   (142) 


5-4.   Derivatives  of  Hamiltonians 

From  the  definition  of  the  Hamiltonian  (11)  and  the  known  value 
of  the  derivative  of  the  state  variables  in  Section  5-2  and  zn  in 
Section  5-3,  we  have 

36"      <      38?       3B? 
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5-4.   Derivatives  of  Hamiltonians 

From  the  definition  of  the  Hamiltonian  (11)  and  the  known  value 
of  the  derivative  of  the  state  variables  in  Section  5-2  and  zn  in 
Section  5-3,  we  have 
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5-5.   Computing  Procedures 

A  suggested 

computational  procedure  to  seek  the  optimal  decisions 

9n  for  a  fixed  x  J 

is  as  follows : 
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Step  6. 

Calculate  OH  ,  ?Hn,  3H  from  equations  (143)  through  (145). 

>e"   nu2   263 

Step  7. 

n 
If  ^H  are  zero  or  less  than  the  allowable  errors  preassigned, 

*" 

then  the  assumed  Q;  are  the  optimal  "values,  otherwise  go  to 

the  next  step. 

Step  8. 

N 
If  x2  is  greater  than  that  computed  in  the  preceding  itera- 

tion, then  one  half  of  the  original  A  9  .  is  used;  otherwise 

the  original  A{)^  is  used. 

Step  9. 

The  new  set  of  decision  (fi.)  new  is  obtained  by 

<>new«  <e?>oldi*8i                   (146) 
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n 
liL  >  0    use  (  +  )  sign 

< 

Then  return  to  step  2  and  repeat  the  computation  until  the  opti- 

mum is  obtained. 

The 

computational  procedure  described  here  is  similar  to  that 

in  Section  5-2  of  PART  ONE.   From  the  experience  in  PART  ONE,  it  is 

believed 

that  the  same  procedure  can  be  applied  to  determination  of 
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the  optimal  condition  of  the  process.   However,  the  numerical  result 
is  yet  to  be  determined  by  the  actual  computer  computation. 
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NOMENCLATURE 

a    =   a  positive  exponent  of  power  rule  for  capital  cost  of  equip- 
ment; a1 ,  a  ,  a,,  and  a  are  the  exponents  for  the  high  pres- 
sure pump,  the  recirculation  pump,  the  turbine,  and  the  mem- 
brane separation  unit,  respectively. 

A    =   cross-section  area  of  the  membrane  separator  unit  in  the  n-th 
stage,  ft  . 

5     K-d  ,3^      /,2,     .2 

c    =   3.05  x  10  constant,  ft  -ft-sec/ft  -hr-psi-cm 

(Sc)1/3  Da 

Ce  =  electrical -power  cost,  $/psi-ft 

Cp  =  the  unit  cost  of  brine  feed,  $/lbm 

Ct  =  the  total  water  cost  per  unit  water  production,  $/lb 

Cn  =  the  various  cost  items  (n=l,  ...,  7)   $/lbm 

d  =  the  diameter  of  the  membrane  tube ,  ft 

Da  =  molecular  diffusivity  of  salt,  cm  /sec. 

E.    =  the  pumping  work  of  the  high  pressure  pump  at  the  n-th  stage 

for  the  general  model;  models  A,  B,  and  C  are  represented  re- 

n    n        n         3 
spectively  by  Ela,  E^,  and  Elc,  psi-ft  /hr. 

n 
Eo    ■   the  pumping  work  of  the  recirculation  pump  at  the  n-th  stage 

for  the  general  model;  models  A,  B,  and  C  are  denoted  by 

n    n        n 
Exa,  E2b>  and  E2C ,  respectively,  psi-ft  /hr 
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B.j       =  the  energy  recovery  from  the  blowdown  turbine  at  the  end  of 
the  process  for  the  general  model;  models  A,  B,  and  C  are 
denoted  by  E3a,  E3b,  and  E3c,  respectively,  psi-ft  /hr. 

t         =     Fanning  friction  factor. 

F    =  the  volumetric  flux  of  water  through  the  membrane,  ft3/ft2-hr. 

H    =  Hamiltonian  functions  of  the  n-th  stage. 

J±        =  the  high  pressure  pump  at  the  n-th  stage. 

n 
J2    =   the  recirculation  pump  at  the  n-th  stage. 

J3   =  the  reject  turbine  at  the  last  stage. 

k    =  the  proportionality  constant  in  the  power  rule  cost  expression 
for  the  equipment  and  kx ,  k2,  kt  and  ks  such  proportionality 
constant  for  the  high  pressure  pump,  the  recirculation  pump, 
the  turbine,  and  the  membrane  separator,  respectively. 

K    =   the  membrane  constant,  ft3/ft2-hr-psi 

L/D  =  the  overall  length-to-diameter  ratio  of  the  membrane  separator. 

M    =  the  mixing  point  in  the  n-th  stage. 

n 
m    =  the  total  number  of  tubes  in  the  membrane  separator  unit  of 

the  n-th  stage. 
MS    =   the  membrane  separator  unit  at  the  n-th  stage 
N    =  the  total  number  of  stages  in  the  system 
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E3 

= 

the  energy  recovery  from  the  blowdown  turbine  at  the  end  of 
the  process  for  the  general  model;  models  A,  B,  and  C  are 
denoted  by  E3a,  E3b,  and  E3c,  respectively,  psi-ft  /hr. 

t 

= 

Fanning  friction  factor. 

F 

= 

the  volumetric  flux  of  water  through  the  membrane,  ft  /ft2-hr. 

Hn 

= 

Hamiltonian  functions  of  the  n-th  stage. 

7n 

= 

the  high  pressure  pump  at  the  n-th  stage. 

n 
J2 

= 

the  recirculation  pump  at  the  n-th  stage. 

J3 

= 

the  reject  turbine  at  the  last  stage,. 

k 

the  proportionality  constant  in  the  power  rule  cost  expression 
for  the  equipment  and  kx ,  k2,  kt  and  ks  such  proportionality 
constant  for  the  high  pressure  pump,  the  recirculation  pump, 
the  turbine,  and  the  membrane  separator,  respectively. 

K 

= 

the  membrane  constant,  ft3/ft2-hr-psi 

L/D 

= 

the  overall  length-to-diameter  ratio  of  the  membrane  separator. 

M" 

= 

the  mixing  point  in  the  n-th  stage. 

n 

m 

= 

the  total  number  of  tubes  in  the  membrane  separator  unit  of 
the  n-th  stage. 

MSn 

= 

the  membrane  separator  unit  at  the  n-th  stage 

N 

the  total  number  of  stages  in  the  system 
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p    =  pressure  within  membrane  separator  chamber  of  the  n-th  stage, 
psi. 

P    =  atmosphere  pressure,  14.7  psi. 

AP   =  pressure  difference  across  the  membrane  at  the  n-th  stage, 
psi. 

q    =  mass  flow  rate  of  brine  solution  discharged  from  the  n-th 
stage,  lbm/hr. 

n 
q^    =   mass  flow  rate  the  brine  entering  the  membrane  separator  of 

the  n-th  stage,  lbm/hr. 

n 
qe    -  mass  flow  rate  of  the  brine  leaving  the  membrane  separator  of 

the  n-th  stage,  lbj^/hr. 

qo  =  mass  flow  rate  of  brine  feed,  lbm/hr. 

R  =  the  flow  rate  of  recycle  stream,  lbm/hr. 

r  =  the  recycle  ratio  at  the  n-th  stage. 

Re  =  Reynolds  number  at  the  n-th  stage. 

S  =  membrane  area  at  the  n-th  stage,  ft 

Sc  ■  Schmidt  number. 

T  =  the  capacity  of  the  equipment 

W  =  flow  of  fresh-water  produced  from  the  n-th  stage,  lb_/hr. 
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wf 

= 

the  total  water  production  from  the  system,  lbE/hr. 

wS 

= 

the  mass  of  the  shell-and-tube  membrane  separator  unit  of 

the 

n-th  stage  for  the  general  model;  models  A,  E,  and  C  are  de- 

noted by  Wsa,  Wsb,  and  Wsc,  respectively,  lbm. 

n 

X 

• 

the  mass  fraction  of  salt  component  in  the  brine  solution 
leaving  the  n-th  stage. 

n 

X. 

1 

= 

the  mass  fraction  of  salt  component  in  the  brine  solution 
entering  the  membrane  separator  of  the  n-th  stage. 

n 

= 

the  mass  fraction  of  salt  component  in  the  brine  solution 
leaving  the  membrane  separator  of  the  n-th  stage. 

n 
xi 

= 

n 

X 

n 
*2 

= 

accumulated  water  cost  at  the  first  n  stages  and  x  =  C  . 

x3 

= 

Qn2 

n 

zi 

adjoint' variables  in  association  with  stage  variables  x  . 

i 

Greek  Letters 

?* 

= 

loss  factor 

?» 

= 

mechanical  efficiency 

7p 

= 

pump  efficiency 

•; 

fiS 
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jU.  =  viscosity  of  the  brine  solution,  lbf/in  . 

9  =  density  of  brine  solution,  lbm/ft  . 

7m  =  density  of  material  of  construction,  lb  /ft  . 

<^m  =  allowable  stress  of  materials  of  construction,  psi. 

\lf  =  capitalization  charge  of  initial  cost  per  hour  in  stream, 
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APPENDIX 

She  PORTRAK  program  symbols,  their  corresponding  mathematical 
notations  ana  their  explanations  are  summarized  in  Sable  A-l.  She 
FCRSHAX  computer  program  (i)  for  the  discrete  maximum  principle  is 
presented  in  Sable  A-2;  the  FC8SRAM  computer  program  (2}  for  the 
simplex  method  is  given  in  Sable  A-J.   The  input  data  and  sample  ot 
results  for  the  computer  program  (l)  are  presented  in  Sable  A-4. 
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Table  A-l-     PROGRAM  SYMBOLS  AND  EXPLANATION 

Program 

Explanation 

K&thematic 

Syabols 

Sysicls 

A 

a*  earn  temperature,  27/T.it°? 

—  7      03.   -  _ 

Ut{A) 

Condensing  area  cost  in  the  n-th  effect 

i> 

B   .  ' 

Coefficient  of  the  Ciausius-Clapeyror*  equation 
including  the  friction  loss,  1.79x10°  lb  ./ft 

3 

-7  / 

CC 

Cooling  water  unit  cost,  5.5673x10  5/lb 

V 

cqs 

c   ,2 
Unit  cost  of  condensing  area,  2. J57xlO  8/" 

c£ 

CH2 

Unit  cost  of  brine  heater,  3.76x10  $/ft 

c>t 

C? 

Heat  capacity  of  water,  1 .  03tu/ib, '? 

p 

CPP 

Unit  pumping  cost,  2.50JxlO  S/ft-lb 

c 

CST' 

Unit  steam' cost,  2.5x10"  S/ib 

c  _ 

St 

-5 

Construction  cost,  2.CSxiO~  S/lOOQgal. 

~6 

2 

Density  of  water,  62.5  Ib/ou.ft. 

? 

Do  (a) 

Increment  cf  decision  variables 

a9 

2?. 

Kaximum  allowable  error  for  --- 

S3 

28E08 

Maximum  allowable  error  for  standard  error 

. 

function  ye ' 

2?:HCS 

-'■*-. 

A 

:-:?s 

Latent  heat  of  steam  at  27^. h°?   and  V5  psia, 
526.5  Btu/lb 

As  . 

?c 

Unit  coot  of  pretrcatmenx,  1:796x10  S/lb 

p 
C 
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Program 

Explanation 

iitio 

Symbols 

"J  •"■"-■- 

s 

PH(a) 

Pump  coat  and  pumping  cost  ir.  the  n-th  effec 

-- 

t"   > 

i 

'ilhe  ratio  of  heat  load  to  seawater  feed 

- 

«./* 

(JP(i>) 

The  value  of  Q  at  vertex  P 

«-s/* 

H 

Ideal  gas  constant,  0.1104  Btu/llj,'!! 

a 

8(n) 

Xuober  of  stages  ir.  the  '.n-th  effect 

-v 

T2ST 

Average  a, in  the  n-th  effect  i 

7.' 

1  5 

atandarc.  error  function  y  '  ,  <■  /      \2/-, 

3   -'iTi^i-^'  /; 

TY(n) 

Direction. of  decision  increment 

Tim, 

Derivative  of  Kamiltor.ian  H  with  respect  to 

n  1 

~ ..  ^ 

E1Y11 

"1 

n  -  vv  1 

derivative  of  Eamiltohiar.  H'L  with  respect  to 

°3 

°  '• 

2»" 

T2YI2, 

2 
Derivative  of  Kasiiitonian  H  with  respect  to 

^7 

.:-:2 

H2Y12 

3»r 

C2YJ2, 

..2   .  . 

derivative  ox  Eaniltoman  :-.  with  respec  to 

aa 

2 

E2YJ2 

vi> 

Derivative  of  Eamiltosian'E  with  respect  to 

fj3 

*H3 
3©| 

u 

Overall  heat  transfer  coefficient, 
510  Btu/nr.ft2."? 

" 

V 

Unit  positive  number,  1 

1 
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Tabic  A-i   (Continued) 

Program 
Symbols 

Explanation 

gmoot 

&tic 

V 

Overall  water  production  rate,  ^000  ga* . /nr . 
or  S.^xiO"5  ib/hr 

-'', 

W?3(n) 

Water  production  rate  in  the  n-th  effect 

V, 

vc 
X3(<0 

Water  cost , S/iOCCgal . 

4 

WCI(r.) 

Function  value  or  water  'Cost  at  vertex  ? 

y, 

Xl(n*l) 

Outlet  concentration  in  the  n-th  effect,  C^W- 

xi 

»<») 

(C,),  at  vertex  ? 
v  f'3           n 

X." 

• 

3»(n+l) 

Temperature  gradient  in  the  n-th  effect, a?" 

--"' 

X2(iifl) 

Accumulated  water  production  cost  in  the  first 
n-th  effects,  S/iOGGgai." 

X2 

X3(a*D 

xin(n-i) 

XlYJ(a-l) 

Brine  temperature  in  the  n-th  effect,  (T#/_ 
derivative  of  x*.  with  respect  to  6" 

,  n  . . ,  v  '      .^n 
jenvative  01  x.  with  respect  91 

5  v. 

*sl 

5X".* 

WYl(a-l) 

derivative  of  AT'*  with  respect  to  01* 

■3  6', 

SSXl(n-2) 

derivative  of  AT  '  with  respective  to  x'.' 

IX'.'"' 

' 

XZYi(n-'l) 

Derivative  of  x"  with  respective  to  o'' 

:'x'i 
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Tabic  A-i   (Continued) 

Program 
Symbols 
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'   Table 

i'.-l    (Continued) 
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pp  Ojjrr  sjd 

; 

Explanation 

Katheraatic 

Symbols 

iy^bolL 

X2Yi(n. 

■1) 

Derivative  of  x" 

with  respoct  to.j). 

25 

'a  «• 

xaxi(«- 

■2) 

Derivative  of  x„ 

with  respect  to  x" 

y  X, 
3  X '-' 

X2X>(n-2) 

.  Derivative  of  x' 

with  respect  to  x"7' 

;  *2 

3  xl" 

?a 

Coefficient  of  reflection 

*' 

• 

T3 

Coefficient  of  contraction 

*8 

Coefficient  of  expansion 

y ' 

«(») 

RgCycie  ratio  in 

the  n-th  effect,  rn 

6*-' 

•  Y3W 

Temperature,  drop 

in,  the'n-th  effect 

S2  "  S2"X 

Zll 

Adjoint  variable, 

zl 

-1 

Z12 

Adjoint  variable, 

2 

z. 

zT 

- 

Z12 

Adjoint  variable, 

A    ■■•' 

n  > 

' 

ZJ1 

Adjoint  variable, 

2y 

J 

1 
zt 

Z'/i 

Adjoint  variable, 

2 

.> 

' 

'   Table 

i'.-l    (Continued) 

156 

pp  Ojjrr  sjd 

; 

Explanation 

Katheraatic 

Symbols 

iy^bolL 

X2Yi(n. 

■1) 

Derivative  of  x" 

with  respoct  to.j). 

25 

'a  «• 

xaxi(«- 

■2) 

Derivative  of  x„ 

with  respect  to  x" 

y  X, 
3  X '-' 

X2X>(n-2) 

.  Derivative  of  x' 

with  respect  to  x"7' 

;  *2 

3  xl" 

?a 

Coefficient  of  reflection 

*' 

• 

T3 

Coefficient  of  contraction 

*8 

Coefficient  of  expansion 

y ' 

«(») 

RgCycie  ratio  in 

the  n-th  effect,  rn 

6*-' 

•  Y3W 

Temperature,  drop 

in,  the'n-th  effect 

S2  "  S2"X 

Zll 

Adjoint  variable, 

zl 

-1 

Z12 

Adjoint  variable, 

2 

z. 

zT 

- 

Z12 

Adjoint  variable, 

A    ■■•' 

n  > 

' 

ZJ1 

Adjoint  variable, 

2y 

J 

1 
zt 

Z'/i 

Adjoint  variable, 

2 

.> 

' 

157 


»ble  A-2 


Co&putar   ?rc" 


9    r! 


.ASH    DISTILLATION    3Y    MAX  J  MUM    PRINCIPLE 
[MENS  I  ON    Y1(4)-.Y3{4>»X1<4)  .X2i41.X3(4)  »DT( 
;  MENS  ION    X1Y3I3)  .DTY1(3).,X2Y1<3)  »X2Y3(31  »D 
INtNiiON    CPRI4I  .0S(  S)  »TV<5)  »PUI4)  >AR(«)  .T-£ 
;SKAT(1CF7.1) 
;S.v,AT  C7E1G.31 
JUKAT  15512.51 
"6512.5). 

?XAT(2X»37HTHS  FOLLOWING  ARE  CJ'TIMUN  CL/TP 
IMAT(2X,6HY1 12)»E12.5»7X,6HY1(3)=E12.5,7X 
5.x  AT ! 2X ,6HY3< 2) "512 . 5 .7X . 6HY3 ( 3 1 =£12.5 ,7X 
!MA7(2X»6HX1(1)»£12.S>7X.6HX1 i21=E12.5»7X 
2X»6.HDT(2)»512.5»7X»6H0T(31  "H12.5»7X 
5.v,AT(2X.6HX3(  1)  "512.5  .7X.6HX3  (2  )  =E12.5,7X 
;;*KA7(  2X.6KX2I2)  "512.5  i7X»6HX2(31  »E12.5>7X 
2X»6HH1Y11»E12.S»7X.6HH2Y12=E12.5,7X 
2X  .6KH1  Y31»E12  .5  >7X  .,6HH2Y32=£12 . 5  »11 
2X»5HTER.MS>20X,7HCCST!$S  >7X.10HP£RCE 
2  X ,  5  HST5AM , 15  X , £12 . 5.5X , £12 . 5 ) 
2  X  »  6  HH  £  A  T  EK  » 1 4 X  t  £ 1 2 . 5  »  5 X  >  E  i  2 • 5 1 
2X , IShCCNDcNSING  AREA»5X»Ei2.5»5Xi£l 
2X  .7.HPUMP I N3..13X  .£1 2 .  5 » 5X , £  12 . 5  1 
2X, 12HCONSTRUC7ION»8X.£12.5»5X,EI2.5 
2X , 12MPR5TREATMENT. 8X  »£12 .5  »5X . £12 .5 
2X  ,5.-;7ES.v,5»:2'X.  11H.1  -ST    Err  EC7  ,6X  >  1 1- 

/: 

2X. IChTEMP.    DR0P>6X»£12 . 5»5X«E12.5  >5 

2X , 11KWATER    PROD. >iX,E12-5>iX,E12.5, 
2X.26HREAD    NEW    OPTIMIZATION    DATA; 

2X,  i-x:  ;^,;=e:  2.  i.::x.2.-j.o=_  i2.  =  ,7x,  or, 

2X,E;2.i>2X,E:2.i,2X>E12.i,2X,E12.3, 


„.NAM  I 

-  (-    ■  '    »    T 

2RXAT 

;.-<mat 


P<  C  ^  ->   i    . 

:<ea3<  : 

S  c  A  3  C  1 
S  E  A  Ci  i  1 


^2 
;X?  i 


^2  x: ; 


BY    KIAI 

4 1 ,  x  i  y  ; 

7X1(2) 


U  ;  *jA  ,  A 
.-6HY1  [4 
>6KY3( 4 
>6HX1  !  i 
>  6  H  D  T  !  4 
*  6  r.  X  5  i  i 
,6.-;X2  i  4 
,6HM3Y1 
X»2HW=£ 
NTASEt/ 


;  kc 

(2x: 


S14 
1  2  1 


X2x: : 2; 


:2 .  a 
L2.5 
.2  .  5 
.2.5 
.2.5 

.2.5; 


2.5,/ ; 


»/! 

2  NC 


X  » E 1 2  •  5  ; 
5X.E12.: 


./; 


.,/; 


2X,;OnCC,\.D.  CCST,6X, 


,5»5X»£12.5»5X,£i2.5; 
2X,  iOHPUMP.  CCST,6X,E:2.5,"5X,Ei2.5,5X,E12.51 
2X.3H3.P.  £L. >SX.£12.5»5X.E12.5.5X»£12.5»/1 
4 )  X 1  (  4 1  >Q*WC»  T£  ( l'l  i  TE  ( 2 1  »"7£  (  3 1 
3!Y1(2)»Y1(3)»Y1(4)»Y3(2).Y3<3) 
4)DS(1)  >OS(2i  »D'S(31  »»S<4)  >DS(5>  >ER 
1 )  U » A  »ri7  »hTS » X3  { 1 1  »  D  •  5  f  2  1  » S  (  3 1  »S  { 4 1  » V 

2  (R»XH  1 1  »CST  tW»CHT  »CCO»CPP 

3  I 3>C5  >  PC >c?  >  cc 


CP-Y3! I 1/HT1 

x-xi  [  i-li !/< i.i-Yi ; :  ;-y:  ; : ; 


'XI 


157 


»ble  A-2 


Co&putar   ?rc" 


9    r! 


.ASH    DISTILLATION    3Y    MAX  J  MUM    PRINCIPLE 
[MENS  I  ON    Y1(4)-.Y3{4>»X1<4)  .X2i41.X3(4)  »DT( 
;  MENS  ION    X1Y3I3)  .DTY1(3).,X2Y1<3)  »X2Y3(31  »D 
INtNiiON    CPRI4I  .0S(  S)  »TV<5)  »PUI4)  >AR(«)  .T-£ 
;SKAT(1CF7.1) 
;S.v,AT  C7E1G.31 
JUKAT  15512.51 
"6512.5). 

?XAT(2X»37HTHS  FOLLOWING  ARE  CJ'TIMUN  CL/TP 
IMAT(2X,6HY1 12)»E12.5»7X,6HY1(3)=E12.5,7X 
5.x  AT ! 2X ,6HY3< 2) "512 . 5 .7X . 6HY3 ( 3 1 =£12.5 ,7X 
!MA7(2X»6HX1(1)»£12.S>7X.6HX1 i21=E12.5»7X 
2X»6.HDT(2)»512.5»7X»6H0T(31  "H12.5»7X 
5.v,AT(2X.6HX3(  1)  "512.5  .7X.6HX3  (2  )  =E12.5,7X 
;;*KA7(  2X.6KX2I2)  "512.5  i7X»6HX2(31  »E12.5>7X 
2X»6HH1Y11»E12.S»7X.6HH2Y12=E12.5,7X 
2X  .6KH1  Y31»E12  .5  >7X  .,6HH2Y32=£12 . 5  »11 
2X»5HTER.MS>20X,7HCCST!$S  >7X.10HP£RCE 
2  X ,  5  HST5AM , 15  X , £12 . 5.5X , £12 . 5 ) 
2  X  »  6  HH  £  A  T  EK  » 1 4 X  t  £ 1 2 . 5  »  5 X  >  E  i  2 • 5 1 
2X , IShCCNDcNSING  AREA»5X»Ei2.5»5Xi£l 
2X  .7.HPUMP I N3..13X  .£1 2 .  5 » 5X , £  12 . 5  1 
2X, 12HCONSTRUC7ION»8X.£12.5»5X,EI2.5 
2X , 12MPR5TREATMENT. 8X  »£12 .5  »5X . £12 .5 
2X  ,5.-;7ES.v,5»:2'X.  11H.1  -ST    Err  EC7  ,6X  >  1 1- 

/: 

2X. IChTEMP.    DR0P>6X»£12 . 5»5X«E12.5  >5 

2X , 11KWATER    PROD. >iX,E12-5>iX,E12.5, 
2X.26HREAD    NEW    OPTIMIZATION    DATA; 

2X,  i-x:  ;^,;=e:  2.  i.::x.2.-j.o=_  i2.  =  ,7x,  or, 

2X,E;2.i>2X,E:2.i,2X>E12.i,2X,E12.3, 


„.NAM  I 

-  (-    ■  '    »    T 

2RXAT 

;.-<mat 


P<  C  ^  ->   i    . 

:<ea3<  : 

S  c  A  3  C  1 
S  E  A  Ci  i  1 


^2 
;X?  i 


^2  x: ; 


BY    KIAI 

4 1 ,  x  i  y  ; 

7X1(2) 


U  ;  *jA  ,  A 
.-6HY1  [4 
>6KY3( 4 
>6HX1  !  i 
>  6  H  D  T  !  4 
*  6  r.  X  5  i  i 
,6.-;X2  i  4 
,6HM3Y1 
X»2HW=£ 
NTASEt/ 


;  kc 

(2x: 


S14 
1  2  1 


X2x: : 2; 


:2 .  a 
L2.5 
.2  .  5 
.2.5 
.2.5 

.2.5; 


2.5,/ ; 


»/! 

2  NC 


X  » E 1 2  •  5  ; 
5X.E12.: 


./; 


.,/; 


2X,;OnCC,\.D.  CCST,6X, 


,5»5X»£12.5»5X,£i2.5; 
2X,  iOHPUMP.  CCST,6X,E:2.5,"5X,Ei2.5,5X,E12.51 
2X.3H3.P.  £L. >SX.£12.5»5X.E12.5.5X»£12.5»/1 
4 )  X 1  (  4 1  >Q*WC»  T£  ( l'l  i  TE  ( 2 1  »"7£  (  3 1 
3!Y1(2)»Y1(3)»Y1(4)»Y3(2).Y3<3) 
4)DS(1)  >OS(2i  »D'S(31  »»S<4)  >DS(5>  >ER 
1 )  U » A  »ri7  »hTS » X3  { 1 1  »  D  •  5  f  2  1  » S  (  3 1  »S  { 4 1  » V 

2  (R»XH  1 1  »CST  tW»CHT  »CCO»CPP 

3  I 3>C5  >  PC >c?  >  cc 


CP-Y3! I 1/HT1 

x-xi  [  i-li !/< i.i-Yi ; :  ;-y:  ; : ; 


'XI 


1SS 


Table  A-2      (Continued) 


l*))/(XlCil*( 1.+Y1I4] ) 1 1/CP 


-xi  (  n/xii :-:  j  i  )/<  I'.tyk  :  ;*x;<  i;/xi: :-: 


42  XI 1 !)■< X*X1 l 1-11 !/( l.+Yli 1 1-Y1 
Y3(41=w7*A'uOGt  (XI  (  3)  +Y1  (' 
OC  43  1=2,4 

x  3 1 : )  ■  x  3  ( I  - 1 1  f  Y  3  ( :  i 

4  3  DT  !  :  i  =  I 0/C?t  Tt  i  .'-1  i  *• 

wA=w/-i  :.-x:i :  i  /xi  <4i  i 

C2»Cri"*0*WA/  !0*<A-X3  H)+0.§*i5T(2)  J  )  ■ 

c;«cst*g»wa/hts 

X2; :  ;=C2->-CiTc: 

DC  44  1-2.4 

w?j  ; :-:  i»wa*xi  ti)«(i./xi(:-i')-i./xi( in 

X»W?R( 2-l)»HT/U 

AS; : )=CC0~X/iD?l tl-TEI I~li-t-Y3(I>/(2.*S(I ) ) ) 

Y*SXP(-riT/(R*X3l  1-1)  )  )-£X?'(-riT/  (R*(X3i  1-1  >+Y3<  I)  )  ! 


P'u( :  >=c??~x:  i 
X2t : i=x2; :-: i 


<XI( 


,/ixi;  i-D*(i.+Yi(  i; ;  i 

1  ;  I  l*Xlf  I  1,1.-/  [HT*X1  i  i-i)  i 

)  tY.  1  (  I )  ) 

>*(0T( I)-TE(I-1)*0.5*Y3(I )/S(] 


x  s  3  «  x : 

X52=X= 


Yl( : 1#B*Y*WA/ 

;-:  i*aR(  :  ;-<-pj: ; ) 

=  X1  i  i  )  ■*  t  X  1  t  I )— XI ( I— 1 
*— C?*X3  (  I  1  *  C  X 1  £  I  - 1  )  •*■ 

i «-DT( : )/(xi ( :-: ;/x: ; 
x3:=x:  i :  i ~x;y:  ii-i)/hxk;i« 
x53=x:  ; ',  !*&TYi!i-i>*(l./x3  ; :-:  j-i./xk  : )  ) 
X32*x33/; ;ot; : >-7£( :-: ;tG.5*y3> : i/S( I ) ;**2; 

DT  i  1  ',  =  CC0'"'H7',>*WA'*  £  X B  I  —  X62  /  /Li 

ya:  =  £x?  ;-.-;r/  (S'Xj; i-n  i  i-exP(-HT/(R*(X3U-l)*Y3( : ; i j i 

Yl ; i )*C??*X1 ( i l*8»YAl*WA/< }*X1 (1-1)1 

X2YH  :-i  t-OTii  i+vi  1 1 ) 

. (1)*X1Y3<  i-l)/( (XI (I >**2)*<D7U )-TE(I-l)+0.5*Y3(U/S( 

,; i)*; i./xii I-i j-i./xi  ti j) 

s3/(2.c*s< : ; "i dt : : >-TE( i-i  i+o.5*Y3U i/S( I ) >»2; 

XA*ri7«EXP(-ri7/!R*{X3(  I'-l)+Y3<  I  1  )  )  )  /(R*(X3  <  1-1  )-t-Y3  I  I )  )**2! 
y.-c?P»Xl(i!*Yl<  I  >»Q*WA*XA/(X1  (I-1)«0) 
44    X?Y3( 1-1) =CCO*HT*WA* <XS1-XS2)/U*Y 
5"=X3  (3  1-Y3  (41-545. 
CT»!0-WAtTE( 3)*W)/BT 

cs«;?<:-cci*watCC*ct 

X2(4>=X2(4)tC6 

AR.'l  )»AR(21tAR(3)+AR(41  ■    _■ 

PU( A>«?Ut2  >-fPU(3>*PU(4) 

■/.'PR !  4  i  =W?R  I  1  i  rWPR (  2  1  tW?R  .' 3  ) 

PCliBCl1r100./X2(40 

PC2-C2"'..G0./'X2;4i 

PCS»AS(11*10G./X2(4) 

PC4»Pu ( 1 i *100«/X2 (4) 

P  C  5  -  C  5  *  '.  0  0  .  /  X  2  i  4  ) 

P  C  6  *  C  6  *  1 0  G  •  /  X  2  £  4  i 


;  i ; 


;  I ; 


1SS 


Table  A-2      (Continued) 


l*))/(XlCil*( 1.+Y1I4] ) 1 1/CP 


-xi  (  n/xii :-:  j  i  )/<  I'.tyk  :  ;*x;<  i;/xi: :-: 


42  XI 1 !)■< X*X1 l 1-11 !/( l.+Yli 1 1-Y1 
Y3(41=w7*A'uOGt  (XI  (  3)  +Y1  (' 
OC  43  1=2,4 

x  3 1 : )  ■  x  3  ( I  - 1 1  f  Y  3  ( :  i 

4  3  DT  !  :  i  =  I 0/C?t  Tt  i  .'-1  i  *• 

wA=w/-i  :.-x:i :  i  /xi  <4i  i 

C2»Cri"*0*WA/  !0*<A-X3  H)+0.§*i5T(2)  J  )  ■ 

c;«cst*g»wa/hts 

X2; :  ;=C2->-CiTc: 

DC  44  1-2.4 

w?j  ; :-:  i»wa*xi  ti)«(i./xi(:-i')-i./xi( in 

X»W?R( 2-l)»HT/U 

AS; : )=CC0~X/iD?l tl-TEI I~li-t-Y3(I>/(2.*S(I ) ) ) 

Y*SXP(-riT/(R*X3l  1-1)  )  )-£X?'(-riT/  (R*(X3i  1-1  >+Y3<  I)  )  ! 


P'u( :  >=c??~x:  i 
X2t : i=x2; :-: i 


<XI( 


,/ixi;  i-D*(i.+Yi(  i; ;  i 

1  ;  I  l*Xlf  I  1,1.-/  [HT*X1  i  i-i)  i 

)  tY.  1  (  I )  ) 

>*(0T( I)-TE(I-1)*0.5*Y3(I )/S(] 


x  s  3  «  x : 

X52=X= 


Yl( : 1#B*Y*WA/ 

;-:  i*aR(  :  ;-<-pj: ; ) 

=  X1  i  i  )  ■*  t  X  1  t  I )— XI ( I— 1 
*— C?*X3  (  I  1  *  C  X 1  £  I  - 1  )  •*■ 

i «-DT( : )/(xi ( :-: ;/x: ; 
x3:=x:  i :  i ~x;y:  ii-i)/hxk;i« 
x53=x:  ; ',  !*&TYi!i-i>*(l./x3  ; :-:  j-i./xk  : )  ) 
X32*x33/; ;ot; : >-7£( :-: ;tG.5*y3> : i/S( I ) ;**2; 

DT  i  1  ',  =  CC0'"'H7',>*WA'*  £  X B  I  —  X62  /  /Li 

ya:  =  £x?  ;-.-;r/  (S'Xj; i-n  i  i-exP(-HT/(R*(X3U-l)*Y3( : ; i j i 

Yl ; i )*C??*X1 ( i l*8»YAl*WA/< }*X1 (1-1)1 

X2YH  :-i  t-OTii  i+vi  1 1 ) 

. (1)*X1Y3<  i-l)/( (XI (I >**2)*<D7U )-TE(I-l)+0.5*Y3(U/S( 

,; i)*; i./xii I-i j-i./xi  ti j) 

s3/(2.c*s< : ; "i dt : : >-TE( i-i  i+o.5*Y3U i/S( I ) >»2; 

XA*ri7«EXP(-ri7/!R*{X3(  I'-l)+Y3<  I  1  )  )  )  /(R*(X3  <  1-1  )-t-Y3  I  I )  )**2! 
y.-c?P»Xl(i!*Yl<  I  >»Q*WA*XA/(X1  (I-1)«0) 
44    X?Y3( 1-1) =CCO*HT*WA* <XS1-XS2)/U*Y 
5"=X3  (3  1-Y3  (41-545. 
CT»!0-WAtTE( 3)*W)/BT 

cs«;?<:-cci*watCC*ct 

X2(4>=X2(4)tC6 

AR.'l  )»AR(21tAR(3)+AR(41  ■    _■ 

PU( A>«?Ut2  >-fPU(3>*PU(4) 

■/.'PR !  4  i  =W?R  I  1  i  rWPR (  2  1  tW?R  .' 3  ) 

PCliBCl1r100./X2(40 

PC2-C2"'..G0./'X2;4i 

PCS»AS(11*10G./X2(4) 

PC4»Pu ( 1 i *100«/X2 (4) 

P  C  5  -  C  5  *  '.  0  0  .  /  X  2  i  4  ) 

P  C  6  *  C  6  *  1 0  G  •  /  X  2  £  4  i 


;  i ; 


;  I ; 
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-2      (Continued) 


701 =X3!  1  )-X3,'2  ) 
T02«X3i"2) -X3(3) 
T03*X3  ;  3  !  -X3  ;  4 ) 

X31«-CCD*W*X1(  1  )*X1Y1  (3  )*V*HT*Xl  ( 1  )*(l-./Xl  (  3)-:l-/.Xl  (*")  1  /U 

X32-X31/;  ;  (XK4  )-xi  ( 1 )  )**2 ) * ;  DT  i  4 i -7£ ; 3  ;  to. 5"Y3  (4)  /  s  ;4 ;  ;  ; 
X5A--Ci3?-»x:  (  1)*Y1  (4)  *B*X1!  1  )*X1Y1 !  3)  *W/  I (XI (4J-X1 ( 1) ) **2) 

X2Y1  13)  =  X2YI  i3)+XB2-rX5A"YA:/(C'i-X'.  I  3)  ) 

PT=-wa*X1 tlJ/(X'l!4>*<Xl(4i-Xl<l))) 

a-=?t-x:y: (3) 

X2Y1 [31-X.2Y1 (3  It (#C-CC)*ATfCC*Q*AT/3T 

X2Y1 ;1)=X2Y1<1  )-CHT*0#W.A*DTYl  f  1  ;  /  1  U*7£i  I  )*(A-X3<  1 i+0.5*OT(2 )  )< 

X3i=xi(l>*il.0/xr(3)-i.0/xr(4) )*HT/(0Tt4)-TE(3)+o.5*Y3{4)/st4 

■X32=-(W*X1!1  >*X1Y3(3>*V/((X1(4)-Xlfl))**2i  )  ■-X3;''CC3/LJ 
X=V.'*xi !1)*X1Y3(3)»YA1*B*X1 (1 }*Yi(4)/( (X1{4)-X1(1> )**2) 
X2Y3(3)=X2Y3(.3)>XB2-CPP*X7(D*X1(3)  ) 
X2Y3(3)"X2Y3(:3J+(PC-CCl*PT*xiY3(3)+CC*(Q*P7*XlY3(3i/aT-CT/3T) 

do  ;;  >j,s 


.  0  ■*■  X " 


D7XlfK)aX0/(X 

X  =  (  - 1  .  /  »'  X  1  £  I  — 
X  =  X  /  ;  3  7 
YA-OTXl 
YA-YA/i 
XA«5X?i- 

Y=-CPP*> 


7£(  :-l  )+X2<  I )*Y1(] ) ) 
(  J-il'*(Xli  1-1  >+Yl<  I)*X1(1.)  )) 
)**2)*1 ./<X1(I-1)*X1(I)))' 
:-1Jt0.5*.Y3(  I  l/.Si  I  )  ) 

*(  :./xi(  :-i  i-i./xi  in) 

i : >-te< :-: >tQ.5*Y3< : i/s; ; i )**2 ; 

/<R*X3( I-l > ) )-EXP(-HT/(R*<X3( I-l>+Y3( I)  >  )  ) 
li»Yl(J  l*3*XA*WA/.(D*Xl  (  I-Di^i- 
^X!  ;  ;<  :  =Y-rCCi3*r;7,fWA*Xl  ( 1  )*  <  X-YA  )  /o 
Xi-£X?  ;  -.-,7/  ( R*X3t  I-l  i  1  )/(X3(  I-rl)**2) 

YA«cXP(-HT/{R*<X3(T-1)+Y3(  I  )>  1  )/(<  X3<  I-l  )+Y3(.i  )  1**2  1 
X 2X3  IK)  "i.O*CPP*Xl  f  1  )*Y1<  I  1  *:WA--<5*.-t7*IX5-YA;  /  (R*&*X1  !  I- 

X£i«xi < i ) *KT* ( 1 .c/xi ( 3J-1.G/X1 14) ) / (DT(4)"-T£(3)+0.5*Y3 

i£52=W»Xl ( 1)*X1(4)/(X1(3)*(X1(4)-X1(1) )**2! 

Y-»-CPP*XK  1  )*Ylt4)*B*XB2*XA/(0*Xl'(3)  ) 

X2X1  (21sX2XK2)+Y-CC0*xa2*XBl/U 

X2X1 (2;=X2X1 (2Sf(PC-CC>*PT*Xl(4>/Xl  (  3  ):*CC*Q*PT*X1  (4)/< 

X2X3(2)=X2X3(21-CC*CT/BT  '    ■ 

213=-X2Y3 !31 /X1Y3C3) 

Z:2=2l3»XU4)/XH3)i-X2Xl<2] 

232«X2X3(2) 

■2:i=Zi2*Xl  (3)  /Xl  (21+X2X1  (1 ) 

Z3:»X2X3< : ;vz 
r.iYiw:i*xiYi<i>*x2Yi<  n 

r  "  Y  3 '  -  v  1 '  ^  X  1  Y  3 
«2Y12"2:2*X1Y: 


x:  ;s;- 


2)-rX2Yl  (2)  ■ 


r.2Y32»212*XiY3I2)»X2Y3(2)+Z32 
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-2      (Continued) 


701 =X3!  1  )-X3,'2  ) 
T02«X3i"2) -X3(3) 
T03*X3  ;  3  !  -X3  ;  4 ) 

X31«-CCD*W*X1(  1  )*X1Y1  (3  )*V*HT*Xl  ( 1  )*(l-./Xl  (  3)-:l-/.Xl  (*")  1  /U 

X32-X31/;  ;  (XK4  )-xi  ( 1 )  )**2 ) * ;  DT  i  4 i -7£ ; 3  ;  to. 5"Y3  (4)  /  s  ;4 ;  ;  ; 
X5A--Ci3?-»x:  (  1)*Y1  (4)  *B*X1!  1  )*X1Y1 !  3)  *W/  I (XI (4J-X1 ( 1) ) **2) 

X2Y1  13)  =  X2YI  i3)+XB2-rX5A"YA:/(C'i-X'.  I  3)  ) 

PT=-wa*X1 tlJ/(X'l!4>*<Xl(4i-Xl<l))) 

a-=?t-x:y: (3) 

X2Y1 [31-X.2Y1 (3  It (#C-CC)*ATfCC*Q*AT/3T 

X2Y1 ;1)=X2Y1<1  )-CHT*0#W.A*DTYl  f  1  ;  /  1  U*7£i  I  )*(A-X3<  1 i+0.5*OT(2 )  )< 

X3i=xi(l>*il.0/xr(3)-i.0/xr(4) )*HT/(0Tt4)-TE(3)+o.5*Y3{4)/st4 

■X32=-(W*X1!1  >*X1Y3(3>*V/((X1(4)-Xlfl))**2i  )  ■-X3;''CC3/LJ 
X=V.'*xi !1)*X1Y3(3)»YA1*B*X1 (1 }*Yi(4)/( (X1{4)-X1(1> )**2) 
X2Y3(3)=X2Y3(.3)>XB2-CPP*X7(D*X1(3)  ) 
X2Y3(3)"X2Y3(:3J+(PC-CCl*PT*xiY3(3)+CC*(Q*P7*XlY3(3i/aT-CT/3T) 

do  ;;  >j,s 


.  0  ■*■  X " 


D7XlfK)aX0/(X 

X  =  (  - 1  .  /  »'  X  1  £  I  — 
X  =  X  /  ;  3  7 
YA-OTXl 
YA-YA/i 
XA«5X?i- 

Y=-CPP*> 


7£(  :-l  )+X2<  I )*Y1(] ) ) 
(  J-il'*(Xli  1-1  >+Yl<  I)*X1(1.)  )) 
)**2)*1 ./<X1(I-1)*X1(I)))' 
:-1Jt0.5*.Y3(  I  l/.Si  I  )  ) 

*(  :./xi(  :-i  i-i./xi  in) 

i : >-te< :-: >tQ.5*Y3< : i/s; ; i )**2 ; 

/<R*X3( I-l > ) )-EXP(-HT/(R*<X3( I-l>+Y3( I)  >  )  ) 
li»Yl(J  l*3*XA*WA/.(D*Xl  (  I-Di^i- 
^X!  ;  ;<  :  =Y-rCCi3*r;7,fWA*Xl  ( 1  )*  <  X-YA  )  /o 
Xi-£X?  ;  -.-,7/  ( R*X3t  I-l  i  1  )/(X3(  I-rl)**2) 

YA«cXP(-HT/{R*<X3(T-1)+Y3(  I  )>  1  )/(<  X3<  I-l  )+Y3(.i  )  1**2  1 
X 2X3  IK)  "i.O*CPP*Xl  f  1  )*Y1<  I  1  *:WA--<5*.-t7*IX5-YA;  /  (R*&*X1  !  I- 

X£i«xi < i ) *KT* ( 1 .c/xi ( 3J-1.G/X1 14) ) / (DT(4)"-T£(3)+0.5*Y3 

i£52=W»Xl ( 1)*X1(4)/(X1(3)*(X1(4)-X1(1) )**2! 

Y-»-CPP*XK  1  )*Ylt4)*B*XB2*XA/(0*Xl'(3)  ) 

X2X1  (21sX2XK2)+Y-CC0*xa2*XBl/U 

X2X1 (2;=X2X1 (2Sf(PC-CC>*PT*Xl(4>/Xl  (  3  ):*CC*Q*PT*X1  (4)/< 

X2X3(2)=X2X3(21-CC*CT/BT  '    ■ 

213=-X2Y3 !31 /X1Y3C3) 

Z:2=2l3»XU4)/XH3)i-X2Xl<2] 

232«X2X3(2) 

■2:i=Zi2*Xl  (3)  /Xl  (21+X2X1  (1 ) 

Z3:»X2X3< : ;vz 
r.iYiw:i*xiYi<i>*x2Yi<  n 

r  "  Y  3 '  -  v  1 '  ^  X  1  Y  3 
«2Y12"2:2*X1Y: 


x:  ;s;- 


2)-rX2Yl  (2)  ■ 


r.2Y32»212*XiY3I2)»X2Y3(2)+Z32 


160   , 

Table   A-2        (Continued) 

H3Y13=Zl'3#XiYl  (3)tX2Y1  (3) 

I flM» 1 ) 66 ,66 ,67 

6  6    TEC 1)"X.01*( <XX(1)»Y1(2)*X1I2) )/U.*yit2) >+Xl(2 

'   /  \J  a  G  6 

T£<2!«1.0075*<  (Xlt21+Yl<3>#Xl<3)  1/U.+YK31  !->-x: 

3  )  ;  /G.G69A 

,TE(31=W.32-r!  (Xl-<3»+Yl(41*Xl(*)>/U.i-Yl(4))TXl(4 

1 /G • 063 G 

y,  =  M  •*•  l 

GO    TO    Hi                             • 

67    Ir<ASS<HiYn)-£R>47,47,51 

47    IF (ABS(HiY31)-£R)4S»48,51 

4  3     : .-  ;  A3S  IH2Y12)  -ER  )  49  ,  49  ,  5  1 

49    IF<AI3S(H2Y32)-£R)5G,  50,51 

5  0    I  F(A35<H3Y13)-gft  1-63,63, 51 

51     IF     irilYll)     137,137,136 

137   TVfl  )=-!.. 

13c     IF     ;H2Y12)     139,139,140 

139    TV(2>«-1. 

140     .,-     !M3Y13)     141,141,142 

, 

141  7v;3)=-i.                                                  .... 

142     I.-     IH1Y31)     143,143,144 

14  3    TV(  4  )■-!.. 

144    I.-     ir.2Y32)     145,145,146 

145    "V(5;--l. 

146     I  ri  WC-X2  i  4  1  )  36  ,36,  35 

36    DO    37    1=1,5 

37    06! I )=C.5*OS; I ) 

Yi;2S=Y117 

Yi;3;=Y127                        •     ' 

. 

Y1(4)«Y13T                                 "    ■ 

Y3;2;=Y317 

Y3,'3)=Y32T 

GO    70    36 

35'  WC«X2<4) 

36    Y117=Y1 12)     • 

Y127-YK3: 

Y13T=Y1!4) 

Y317=Y3(2) 

Y327=Y3(3) 

Yl(2i=YH2!-TV<l)*0'S(l) 

Y1(3)*Y1(3)-TV<2)*DS{2)  ■ 

Y3(2 )»Y3!2 I -TV! 4 ! »0S ! 4 ) 

Y3!3i»Y3{3).-TV(&)»DS(5)    ■ 

GO    70    71 

.6i'A.:<;7E:3>25!Xl(4),0,X2:4) 

160   , 

Table   A-2        (Continued) 

H3Y13=Zl'3#XiYl  (3)tX2Y1  (3) 

I flM» 1 ) 66 ,66 ,67 

6  6    TEC 1)"X.01*( <XX(1)»Y1(2)*X1I2) )/U.*yit2) >+Xl(2 

'   /  \J  a  G  6 

T£<2!«1.0075*<  (Xlt21+Yl<3>#Xl<3)  1/U.+YK31  !->-x: 

3  )  ;  /G.G69A 

,TE(31=W.32-r!  (Xl-<3»+Yl(41*Xl(*)>/U.i-Yl(4))TXl(4 

1 /G • 063 G 

y,  =  M  •*•  l 

GO    TO    Hi                             • 

67    Ir<ASS<HiYn)-£R>47,47,51 

47    IF (ABS(HiY31)-£R)4S»48,51 

4  3     : .-  ;  A3S  IH2Y12)  -ER  )  49  ,  49  ,  5  1 

49    IF<AI3S(H2Y32)-£R)5G,  50,51 

5  0    I  F(A35<H3Y13)-gft  1-63,63, 51 

51     IF     irilYll)     137,137,136 

137   TVfl  )=-!.. 

13c     IF     ;H2Y12)     139,139,140 

139    TV(2>«-1. 

140     .,-     !M3Y13)     141,141,142 

, 

141  7v;3)=-i.                                                  .... 

142     I.-     IH1Y31)     143,143,144 

14  3    TV(  4  )■-!.. 

144    I.-     ir.2Y32)     145,145,146 

145    "V(5;--l. 

146     I  ri  WC-X2  i  4  1  )  36  ,36,  35 

36    DO    37    1=1,5 

37    06! I )=C.5*OS; I ) 

Yi;2S=Y117 

Yi;3;=Y127                        •     ' 

. 

Y1(4)«Y13T                                 "    ■ 

Y3;2;=Y317 

Y3,'3)=Y32T 

GO    70    36 

35'  WC«X2<4) 

36    Y117=Y1 12)     • 

Y127-YK3: 

Y13T=Y1!4) 

Y317=Y3(2) 

Y327=Y3(3) 

Yl(2i=YH2!-TV<l)*0'S(l) 

Y1(3)*Y1(3)-TV<2)*DS{2)  ■ 

Y3(2 )»Y3!2 I -TV! 4 ! »0S ! 4 ) 

Y3!3i»Y3{3).-TV(&)»DS(5)    ■ 

GO    70    71 

.6i'A.:<;7E:3>25!Xl(4),0,X2:4) 

161 


iblo  A-2  (Continued) 


W  S  : 

ws  : 

WR  I 

6*  X3.; 

WS  I 

wr  : 

ws : 
ws: 

wr: 
wr  : 
wr  : 

ws: 
ws: 
ws : 

ws: 
ws: 
ws: 


ws : 

WS  I 

ws: 

ws: 
ws: 

GO 


r£; 3,6  jy:  (2  )  ,y:  (3  ;  ,v:  ;•-) 

r£(3.7)Y3(21  ,Y3.'3  !  >Y3(',) 
rei3>8JXlli ) ,X1<2) >X1 (3) 
*  £(3. 9  JOT  (2)  ,DT(  3)  »OT  (4) 

>4   :=2.4 

:  ) =  X3 (  ;  ) -v60. 

rE(3»10)X3(2.)»X3(3>  ,X3<4) 

*E(3»ll)X2!l)tX2t2) ,X2<3) 

■£(3>12)H1YU.H2Y12,H3Y13 

•£(3,13>H1Y31,H2Y3.WPR(4) 

'£(3»26}DS< 1 i  »OSf 2 ) , 0S( 3  > »DS(4) »DS(5) 

'E(3,i4J 

■£;3,i5ic:.pc: 

'£;3,:6iC2,?C2 
£(3il7>AR(  i;  ,?C3- 
'£(3,18  )?o'(  1 )  >?C4 
'E(3,:9)C5.PC5 
E(3»20)C6rPC6 
£(3,21) 

E(3»22.)T01.TD2  .7  33 
E(3.23>WPR(1) .WPR{2>.WPR(3)     . ' 
E(3.27)AR(2)  »'AR(3)  ,AR<4) 
E(3»28)PU("2)  «PU(3)  ,.PU{4) 
£;3,29>TE(1) ,T£;2) ,7£(3) 
£(3.241 
0    4C        '     • 


161 


iblo  A-2  (Continued) 


W  S  : 

ws  : 

WR  I 

6*  X3.; 

WS  I 

wr  : 

ws : 
ws: 

wr: 
wr  : 
wr  : 

ws: 
ws: 
ws : 

ws: 
ws: 
ws: 


ws : 

WS  I 

ws: 

ws: 
ws: 

GO 


r£; 3,6  jy:  (2  )  ,y:  (3  ;  ,v:  ;•-) 

r£(3.7)Y3(21  ,Y3.'3  !  >Y3(',) 
rei3>8JXlli ) ,X1<2) >X1 (3) 
*  £(3. 9  JOT  (2)  ,DT(  3)  »OT  (4) 

>4   :=2.4 

:  ) =  X3 (  ;  ) -v60. 

rE(3»10)X3(2.)»X3(3>  ,X3<4) 

*E(3»ll)X2!l)tX2t2) ,X2<3) 

■£(3>12)H1YU.H2Y12,H3Y13 

•£(3,13>H1Y31,H2Y3.WPR(4) 

'£(3»26}DS< 1 i  »OSf 2 ) , 0S( 3  > »DS(4) »DS(5) 

'E(3,i4J 

■£;3,i5ic:.pc: 

'£;3,:6iC2,?C2 
£(3il7>AR(  i;  ,?C3- 
'£(3,18  )?o'(  1 )  >?C4 
'E(3,:9)C5.PC5 
E(3»20)C6rPC6 
£(3,21) 

E(3»22.)T01.TD2  .7  33 
E(3.23>WPR(1) .WPR{2>.WPR(3)     . ' 
E(3.27)AR(2)  »'AR(3)  ,AR<4) 
E(3»28)PU("2)  «PU(3)  ,.PU{4) 
£;3,29>TE(1) ,T£;2) ,7£(3) 
£(3.241 
0    4C        '     • 


Table  A-}   Computer  Program  (2) 

FLASH  DISTILLATION  BY  SIMPLEX  METHOD  AND  MAXIMUM  PRINCIPLE 
DIMENSION  Yl{4)»Y3l4)»Xl(4).X2<4),X3(4i,DT<4),XlYl<3),S<4) 
DIMENSION  X1Y3<3).X2Y1(3>.DTY1(3).XIY3<3),DTX112),X2X3(2>.X2X1(2) 
DIMENSION  WPR(4),DS(5),TV(5>  .PUt4)  >AR(4)  ,TE<3)  ,QF(6)  ,XF[6)  ,'ACI(6) 
DIMENSION  OF(6),XF(6),WCI(6) 

1  FORMAT(10F7.1 ) 

2  FORMATI7E10.3) 

3  FORMATI5E12.5) 

4  FORMAT(6E12.5) 

5  FCRMAT(2X,37HTHE  FOLLOWING  ARE  OPTIMUN  OUTPUT  DATA./) 

6  F0RMAT<2X,6HY1(2)=E12.5,7X,6HY1<3).  =  E12.5,7X,6HY1(4)=E12.5) 

7  F0RMAT(2X,6HY3(2)=E12.5,7X,6HY3(3)=E12.5,7X,6HY3!4)=E12.5) 

8  FORMAT (2X,6HXK1)=E 12.5. 7X.6HX1 (2) =£12. '5,7X,6HX1(3)=E12.5) 

9  FORMAT<2X,6HDT<2)=E12.5.7X,6HDT(3)=E12.b,7X,6HDT<4)=-12.5> 

10  FORMAT<2X,6HX2  12)=E12.5,7X,6HX2(3)=E12.5,7X,6HX2(4)=E12.5) 

11  F0RMAT(2X,6HX3C1)=E12.5,7X,6HX3<2>=E12.5,7X,SHX3(3)=E12.5) 

12  FORMAT(2X,6HHlYll=E12.5,7X,6HH2Y12=E12.5,7X,6riH3Y13=^2.5) 

13  F0RMAT(2X,6HH1Y31=E12,5,7X,6HH2Y32=E12.5,7X,2HW=E12.5,/) 

14  FORMAT <2X,5HTERMS.20X,7HCGST(S) ,7X , lOHPERCENTAGr , / ) 

15  FORMAT<2X,5HSTEAM,15X,E12.5,5X,E12.5) 

16  FORMAT(2X,6HHEATER,14X,E12.5,5X,E12.5l 

17  FORMAT(2X,15HCCNDENSING  AREA , 5X , E12. 5  ,5X  .  El  2 . 5 ) 

18  FORMAT(2X,7HPUMPING,13X,E12.5,5X,E12.5) 

19  FORMAT(2X,12HCONSTRUCTION,8X,E12.5,5X,E12.5) 
2C  F0RMAT(2X,12HPRETREATMENT,8X,E12.b.5X,E12.5>/) 

21  FORMATI2X.5HTERMS.12X.11H1  ST  EFFECT ,6X , 1 1H2  ND  EFFECT ,6X . 11H3  RD 
1EFFECT,/) 

22  FORMAT(2X,10HTEMP.  DROP ,6X ,E 12. 5 , 5X ,E12 . 5 , 5X , E12 . 5 ) 
,23  FORMATI2X.11HWATER  PROD.  ,5X  ,  E12  .  5  ,  5X  ,  E12.  5  ,  5X  ,  El  2  .  5  ,/  ) 
24  F0RMAT(2X,26HREAD  MEW  OPTIMIZATION  DATA) 

2  5    FORMAT  (2X,6HX1(  4.)  =E  12.5.1  IX. 2HO  =  E12.  5,  7X.6HX2!4)=£I2.5) 

2  6    FORMAT(2X,E12.5,2X.E12.5,2X.E12.5,2X,E12.5,2X,E12.5,/) 

27  FORMATI2X.10HCOND.    COST ,6X , E12 . 5 ,5X  ,E12 . 5  ,5X , E12 . 5 ) 

28  FORMAT(2X,10HPUMP.     COST ,6X , El  2. 5 . 5X ,E12 . 5 , 5X , E 1 2 . 5 ) 

29  FORMAT12X.8HB.P.     EL. .8X.E12 . 5 ,5X .E12. 5 , 5X , E12 . 5 , / ) 
31    FORMAT  (2X.6HXK  4)  =E 12 . 5 ,2X , 2riQ=El 2.5 ,2X , 3HTN=E12 . 5 ) 

3  3    FORMATI2X,6HXF(I)=E12.5,2X,6HC!F(I)=E12.5,2X,7HWCI(I)=P12.5,/) 
34    FORMAT<2X,3HTN=F4.1) 

123  FORMAT  (2X.6HXK  4)  =E12 . 5 ,2X ,2HQ=E1 2 . 5 , 2X ,3HTX=E12 . 5 ) 

124  FORMAT(3E12.5) 

71    READ    4,TX,TN,TA.Tb.TR,T2 

125  RFAD    124.XK4)  ,0,22 

126  READ    3.WCTE11 ) ,TE(2) ,TE(3) ,ER 

READ    3.Y1(2).Y1(3),Y1(4),Y3(2),Y3(3) 

READ    3,(DS( I ) .1=1,5) 

READ    4,(TV( I ) ,1=1,5) ,2Z 

READ    l.U.A.HT.HTS.X3(l) ,D, S( 2 ) »S( 3 ) »S (4 ) > V 

READ    2.R.XK  1)  ,CST  .W  .CHT  .CCD  .CPP 

READ    4,B,C5.PC,CP,CC,ERR0R 

M  =  l 

L  =  l 

41  DO    42     1=2,3 
X=EXPF(-CP*Y3( I 1/riT) 

42  XI  (  I  )  =  (X*X1(  1-1)  )/(l.  +  Yl(I  )-YK I )*X) 

Y3(4)=HT»L0GF( <Xl<3)+Yl{4)*Xl(4))/(XK41*(l.+Yli4))))/CP 

DO    43    1=2,4 

X3(  I)=X3(I-1)+Y3(I  ) 

43  DTI  I )  =  (D/CP+TE(I-1)*(1.-X1 (1)/X1( 1-1)  ))/(l.+YK  I  ) *X 1 ( 1 ) /Xl ( I -1 1 ) 


Table  A-}   Computer  Program  (2) 

FLASH  DISTILLATION  BY  SIMPLEX  METHOD  AND  MAXIMUM  PRINCIPLE 
DIMENSION  Yl{4)»Y3l4)»Xl(4).X2<4),X3(4i,DT<4),XlYl<3),S<4) 
DIMENSION  X1Y3<3).X2Y1(3>.DTY1(3).XIY3<3),DTX112),X2X3(2>.X2X1(2) 
DIMENSION  WPR(4),DS(5),TV(5>  .PUt4)  >AR(4)  ,TE<3)  ,QF(6)  ,XF[6)  ,'ACI(6) 
DIMENSION  OF(6),XF(6),WCI(6) 

1  FORMAT(10F7.1 ) 

2  FORMATI7E10.3) 

3  FORMATI5E12.5) 

4  FORMAT(6E12.5) 

5  FCRMAT(2X,37HTHE  FOLLOWING  ARE  OPTIMUN  OUTPUT  DATA./) 

6  F0RMAT<2X,6HY1(2)=E12.5,7X,6HY1<3).  =  E12.5,7X,6HY1(4)=E12.5) 

7  F0RMAT(2X,6HY3(2)=E12.5,7X,6HY3(3)=E12.5,7X,6HY3!4)=E12.5) 

8  FORMAT (2X,6HXK1)=E 12.5. 7X.6HX1 (2) =£12. '5,7X,6HX1(3)=E12.5) 

9  FORMAT<2X,6HDT<2)=E12.5.7X,6HDT(3)=E12.b,7X,6HDT<4)=-12.5> 

10  FORMAT<2X,6HX2  12)=E12.5,7X,6HX2(3)=E12.5,7X,6HX2(4)=E12.5) 

11  F0RMAT(2X,6HX3C1)=E12.5,7X,6HX3<2>=E12.5,7X,SHX3(3)=E12.5) 

12  FORMAT(2X,6HHlYll=E12.5,7X,6HH2Y12=E12.5,7X,6riH3Y13=^2.5) 

13  F0RMAT(2X,6HH1Y31=E12,5,7X,6HH2Y32=E12.5,7X,2HW=E12.5,/) 

14  FORMAT <2X,5HTERMS.20X,7HCGST(S) ,7X , lOHPERCENTAGr , / ) 

15  FORMAT<2X,5HSTEAM,15X,E12.5,5X,E12.5) 

16  FORMAT(2X,6HHEATER,14X,E12.5,5X,E12.5l 

17  FORMAT(2X,15HCCNDENSING  AREA , 5X , E12. 5  ,5X  .  El  2 . 5 ) 

18  FORMAT(2X,7HPUMPING,13X,E12.5,5X,E12.5) 

19  FORMAT(2X,12HCONSTRUCTION,8X,E12.5,5X,E12.5) 
2C  F0RMAT(2X,12HPRETREATMENT,8X,E12.b.5X,E12.5>/) 

21  FORMATI2X.5HTERMS.12X.11H1  ST  EFFECT ,6X , 1 1H2  ND  EFFECT ,6X . 11H3  RD 
1EFFECT,/) 

22  FORMAT(2X,10HTEMP.  DROP ,6X ,E 12. 5 , 5X ,E12 . 5 , 5X , E12 . 5 ) 
,23  FORMATI2X.11HWATER  PROD.  ,5X  ,  E12  .  5  ,  5X  ,  E12.  5  ,  5X  ,  El  2  .  5  ,/  ) 
24  F0RMAT(2X,26HREAD  MEW  OPTIMIZATION  DATA) 

2  5    FORMAT  (2X,6HX1(  4.)  =E  12.5.1  IX. 2HO  =  E12.  5,  7X.6HX2!4)=£I2.5) 

2  6    FORMAT(2X,E12.5,2X.E12.5,2X.E12.5,2X,E12.5,2X,E12.5,/) 

27  FORMATI2X.10HCOND.    COST ,6X , E12 . 5 ,5X  ,E12 . 5  ,5X , E12 . 5 ) 

28  FORMAT(2X,10HPUMP.     COST ,6X , El  2. 5 . 5X ,E12 . 5 , 5X , E 1 2 . 5 ) 

29  FORMAT12X.8HB.P.     EL. .8X.E12 . 5 ,5X .E12. 5 , 5X , E12 . 5 , / ) 
31    FORMAT  (2X.6HXK  4)  =E 12 . 5 ,2X , 2riQ=El 2.5 ,2X , 3HTN=E12 . 5 ) 

3  3    FORMATI2X,6HXF(I)=E12.5,2X,6HC!F(I)=E12.5,2X,7HWCI(I)=P12.5,/) 
34    FORMAT<2X,3HTN=F4.1) 

123  FORMAT  (2X.6HXK  4)  =E12 . 5 ,2X ,2HQ=E1 2 . 5 , 2X ,3HTX=E12 . 5 ) 

124  FORMAT(3E12.5) 

71    READ    4,TX,TN,TA.Tb.TR,T2 

125  RFAD    124.XK4)  ,0,22 

126  READ    3.WCTE11 ) ,TE(2) ,TE(3) ,ER 

READ    3.Y1(2).Y1(3),Y1(4),Y3(2),Y3(3) 

READ    3,(DS( I ) .1=1,5) 

READ    4,(TV( I ) ,1=1,5) ,2Z 

READ    l.U.A.HT.HTS.X3(l) ,D, S( 2 ) »S( 3 ) »S (4 ) > V 

READ    2.R.XK  1)  ,CST  .W  .CHT  .CCD  .CPP 

READ    4,B,C5.PC,CP,CC,ERR0R 

M  =  l 

L  =  l 

41  DO    42     1=2,3 
X=EXPF(-CP*Y3( I 1/riT) 

42  XI  (  I  )  =  (X*X1(  1-1)  )/(l.  +  Yl(I  )-YK I )*X) 

Y3(4)=HT»L0GF( <Xl<3)+Yl{4)*Xl(4))/(XK41*(l.+Yli4))))/CP 

DO    43    1=2,4 

X3(  I)=X3(I-1)+Y3(I  ) 

43  DTI  I )  =  (D/CP+TE(I-1)*(1.-X1 (1)/X1( 1-1)  ))/(l.+YK  I  ) *X 1 ( 1 ) /Xl ( I -1 1 ) 


Tabic     A-3    (continued) 
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WA  =  W/(  l.-XK  1I/X1I4)  ) 

C2  =  CHT*CJ*WA/(U*(A-X4U  >+0.5*DT(2  1  )  ) 
C1  =  CST*C0*WA/HTS 
X3( 1)=C2+C5+C1 
r  2    44    1=2.4 

V.PRI  I-1)=WA*X1  (1)*<1./X1(I-1»-1./X1U  )  ) 
X=WPR< 1-1 )*HT/U 

AR(  I  )=CCD*X/(DT(I )-TE(  I-11+Y3I I >/<2.*S(I )  ) ) 
Y=FXPF(-HT/(R*X3< 1-1 ) )  )-EXPF(-HT/ ( R*(X3 I  I -1 ) +Y3 ( I )  )  )  ) 
PU( I >=CPP*X1 ( 1 )*Y1 t I )*B*Y*WA/(D*X1 ( 1-1 ) ) 
X2( I )=X2( 1-1 )+AR( I )+PUt I) 

XlYlf  1-1  )=X1  t  I  )*(X1(  I)-X1(  I-l)  >/(Xi(  J-11*(1.  +  YKI)  )  ) 
X1Y3(I-1)=-CP*X1(I)*(X1(I-1)+Y1(I)*X1(I ) )/(HT*Xi(I-D) 
DTY1(I-1)=-DT< I)/(X1U-1)/X1(1)+Y1(I)) 

XBl=Xl(l)*XlYKI-l)/( (XI (I )**2)*(DT(I)-TE(I-1)+0.5*Y3( I )/SiI ) )) 
XB3=Xl(l)*X2Yl(I-l)*(l./xi(I-li-l./XKI)) 
XB2=XB3/( <DT< I )-TE( I-1)+0.5*Y3( I >/S< I ) )**2) 
OT(l)=CCD*HT*WA*(XBl-XB2)/U 

YA1=EXPF(-HT/(R*X3( I -1 ) ) J-EXPF ( -HT/ ( R* ( X3 ( I -1 ) +Y3 ( I ) ) ) ) 
Y1(1)=CPP*X1(1)*B*YA1*WA/(D*X1( I-l)  ) 
X2Y1(I-1)=DT(1)+Y1(1 ) 

XB1=X1  (l)*XlY3(I-l)/(  (XI  (I  )*«2)*<X2(I  )-TE[  I-1)+0.5*Y3<  I  >/S(I)  !  ) 
XB3=X1(1)*(1./X1(I-1)-1./X1(I)) 

XB2=XB3/(2.0*S(I ) * ( DT ( I ) -TE ( I -1 ) +0 . 5*Y3 ( I )/S(I) )**2) 
XA=HT*EXPF(-HT/(R*(X3( I-l ) +Y3 ( I ) 1 ) ) / ( R* ( X3 ( I-l )+ Y3 ( I ) ) **2 ) 
44  X2Y3( I-1)=CCD*HT*WA*(XB1-XB21/U+Y 
BT=X2(3)+Y3(4)-545. 
CT=<Q*WA+TE<3)*W)/BT 
C6=(PC-CC)*WA+CC*CT 
X2(4)=X2(4)+C6 
AR( 1 )=AR( 2 )+AR(3)+AR (4) 
PU(1)=PU(2)+PU(3)+PU(4) 
WPR ( 4 ) =  WPR ( 1 ) +WPR ( 2 ) +WPR ( 3 ) 
PC1=C1*100./X2(4) 
•   PC2=C2*100./X2(4) 
PC3=AR( 1)*10C./X2(4) 
PC4=PU(1)*100./X2(4) 
PC5=C5*1U0./X2(4) 
PC6=C6*10G./X2(4) 
TD1=X3(1)-X3(2) 
TD2=X3(2)-X3(3) 
TD3=X3(3)-X3I4) 

XB1=-CCD*W*X1(1)*X1Y1(3)*V*HT*XI(1)*(1./X1(3)-1./X1(4) ) /U 
>  32  =  XBl/(  ( (XlU)-Xlll)  >**2)*(DT(4)-TE(3)+0.5*Y3(4)/S(4)  1  ) 
X5A  =  -CPP*XHl)*Yl(4)*B*Xl(l)*XlYK3)*W/( ( XI ( 4 ) -XI ( 1 ) ) **2 ) 
X2YK3)=X2Y1(3)+XB2  +  X5A*YA1/(D*X1(3)  ) 
PT=-WA*X1(1)/(X1(4)*(X1(4)-X1(1) ] ) 
AT=PT*X1Y1 (3) 

X2Y1(3)=X2Y1(3)+(PC-CC)*AT+CC*Q*AT/DT 
X2Y1(1)=X2Y1(1)-CMT*Q*WA*0TY1(1 ) / ( U*TE( 1 )*( A-X3 ( 1 ) +0.5*DT I  2 ) )**2) 

XBl-Xl(l)*(l.-u/Xl(3)-1.0/Xl(4))*HT/(DT(4)-Tt(3)  +  0.5*Y3(4)/M4>) 
XB2=-(w*XK 1)»X1Y3(3)*V/I (XI ( 4 ) -X 1 ! 1 ) )  **2  ) )*Xbl*CCD/U 
X=W*Xl(l)*XlY3(3)*YAl*B*Xl(l)*Yl(4)/( ( X 1 ( 4 ) -XI ( 1 ) ) **2 ] 
X2 Y3 ( 3 ) =X3Y2 ( 3 ) +XB2-CPP*X/ ( D*X1 ( 3 ) ) 

X2Y3(3>=X2Y3(3)+(PC-CC)*PT*XlY3(3)+CC*iQ*PT*XlY3i3)/«T-C-/bT) 
DO  45  1=3,4 

K  =  I-2 

XQ=1.0*X1(1)*(TE( I-l )+DT( I )#Y1 ( I ) ) 


Tabic     A-3    (continued) 
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WA  =  W/(  l.-XK  1I/X1I4)  ) 

C2  =  CHT*CJ*WA/(U*(A-X4U  >+0.5*DT(2  1  )  ) 
C1  =  CST*C0*WA/HTS 
X3( 1)=C2+C5+C1 
r  2    44    1=2.4 

V.PRI  I-1)=WA*X1  (1)*<1./X1(I-1»-1./X1U  )  ) 
X=WPR< 1-1 )*HT/U 

AR(  I  )=CCD*X/(DT(I )-TE(  I-11+Y3I I >/<2.*S(I )  ) ) 
Y=FXPF(-HT/(R*X3< 1-1 ) )  )-EXPF(-HT/ ( R*(X3 I  I -1 ) +Y3 ( I )  )  )  ) 
PU( I >=CPP*X1 ( 1 )*Y1 t I )*B*Y*WA/(D*X1 ( 1-1 ) ) 
X2( I )=X2( 1-1 )+AR( I )+PUt I) 

XlYlf  1-1  )=X1  t  I  )*(X1(  I)-X1(  I-l)  >/(Xi(  J-11*(1.  +  YKI)  )  ) 
X1Y3(I-1)=-CP*X1(I)*(X1(I-1)+Y1(I)*X1(I ) )/(HT*Xi(I-D) 
DTY1(I-1)=-DT< I)/(X1U-1)/X1(1)+Y1(I)) 

XBl=Xl(l)*XlYKI-l)/( (XI (I )**2)*(DT(I)-TE(I-1)+0.5*Y3( I )/SiI ) )) 
XB3=Xl(l)*X2Yl(I-l)*(l./xi(I-li-l./XKI)) 
XB2=XB3/( <DT< I )-TE( I-1)+0.5*Y3( I >/S< I ) )**2) 
OT(l)=CCD*HT*WA*(XBl-XB2)/U 

YA1=EXPF(-HT/(R*X3( I -1 ) ) J-EXPF ( -HT/ ( R* ( X3 ( I -1 ) +Y3 ( I ) ) ) ) 
Y1(1)=CPP*X1(1)*B*YA1*WA/(D*X1( I-l)  ) 
X2Y1(I-1)=DT(1)+Y1(1 ) 

XB1=X1  (l)*XlY3(I-l)/(  (XI  (I  )*«2)*<X2(I  )-TE[  I-1)+0.5*Y3<  I  >/S(I)  !  ) 
XB3=X1(1)*(1./X1(I-1)-1./X1(I)) 

XB2=XB3/(2.0*S(I ) * ( DT ( I ) -TE ( I -1 ) +0 . 5*Y3 ( I )/S(I) )**2) 
XA=HT*EXPF(-HT/(R*(X3( I-l ) +Y3 ( I ) 1 ) ) / ( R* ( X3 ( I-l )+ Y3 ( I ) ) **2 ) 
44  X2Y3( I-1)=CCD*HT*WA*(XB1-XB21/U+Y 
BT=X2(3)+Y3(4)-545. 
CT=<Q*WA+TE<3)*W)/BT 
C6=(PC-CC)*WA+CC*CT 
X2(4)=X2(4)+C6 
AR( 1 )=AR( 2 )+AR(3)+AR (4) 
PU(1)=PU(2)+PU(3)+PU(4) 
WPR ( 4 ) =  WPR ( 1 ) +WPR ( 2 ) +WPR ( 3 ) 
PC1=C1*100./X2(4) 
•   PC2=C2*100./X2(4) 
PC3=AR( 1)*10C./X2(4) 
PC4=PU(1)*100./X2(4) 
PC5=C5*1U0./X2(4) 
PC6=C6*10G./X2(4) 
TD1=X3(1)-X3(2) 
TD2=X3(2)-X3(3) 
TD3=X3(3)-X3I4) 

XB1=-CCD*W*X1(1)*X1Y1(3)*V*HT*XI(1)*(1./X1(3)-1./X1(4) ) /U 
>  32  =  XBl/(  ( (XlU)-Xlll)  >**2)*(DT(4)-TE(3)+0.5*Y3(4)/S(4)  1  ) 
X5A  =  -CPP*XHl)*Yl(4)*B*Xl(l)*XlYK3)*W/( ( XI ( 4 ) -XI ( 1 ) ) **2 ) 
X2YK3)=X2Y1(3)+XB2  +  X5A*YA1/(D*X1(3)  ) 
PT=-WA*X1(1)/(X1(4)*(X1(4)-X1(1) ] ) 
AT=PT*X1Y1 (3) 

X2Y1(3)=X2Y1(3)+(PC-CC)*AT+CC*Q*AT/DT 
X2Y1(1)=X2Y1(1)-CMT*Q*WA*0TY1(1 ) / ( U*TE( 1 )*( A-X3 ( 1 ) +0.5*DT I  2 ) )**2) 

XBl-Xl(l)*(l.-u/Xl(3)-1.0/Xl(4))*HT/(DT(4)-Tt(3)  +  0.5*Y3(4)/M4>) 
XB2=-(w*XK 1)»X1Y3(3)*V/I (XI ( 4 ) -X 1 ! 1 ) )  **2  ) )*Xbl*CCD/U 
X=W*Xl(l)*XlY3(3)*YAl*B*Xl(l)*Yl(4)/( ( X 1 ( 4 ) -XI ( 1 ) ) **2 ] 
X2 Y3 ( 3 ) =X3Y2 ( 3 ) +XB2-CPP*X/ ( D*X1 ( 3 ) ) 

X2Y3(3>=X2Y3(3)+(PC-CC)*PT*XlY3(3)+CC*iQ*PT*XlY3i3)/«T-C-/bT) 
DO  45  1=3,4 

K  =  I-2 

XQ=1.0*X1(1)*(TE( I-l )+DT( I )#Y1 ( I ) ) 


Table      A-3    (oontii 

DTX1(K)=XQ/(X1<  1-1  )*(X1  (  I-D+Yl  (  I  |*X1  (1  )  )  ) 

X=(-1./(X1(I-1 !**2)+l./(Xl ( 1-1 )*X1( 1 )) ) 

X=X/(DT( I )-TE( 1-1)+0.5*Y3( I )/S(I) ) 

YA  =  DTX1!K)*(  l./XK  I-D-l./Xl  (I  )  ) 

YA=YA/( (DT(I)-TE(I-l)+0.5*Y3<l)/S< I))**2) 

XA«EXPF(-HT/{R#X3(  I'-l)  )  J-EXPF <-HT/ (R*(X3<  I-D+Y3  ( I  )  )  > ) 

Y=-CPP*X1(1)*Y1( I )*B*XA*WA/(D*X1( I-l)**2) 

X2X1(K)=Y+CCD*HT*WA*X1 ( 1)*(X-YA)/U 

X5=EXPF(-HT/(R*X3( 1-1 ) ) ) / ( X3 ( I -1 ) **2 ) 

YA  =  EXPF(-HT/(R*(X3(  I -1  )+Y3  (  I  )  )  )  )  /  1  (  X3  U-l  )+Y3  (  I  )  )  **2  ) 

45 

X2X3(K)=1.0*CPP*X1( 1>*Y1(I )*WA*B*HT*(X5-YA)/(R*D*X1( 1-1) ) 
XB1=X1(1)*HT*(1.0/X1(3)-1.0/X1(4) ) / ( DJ ( 4 ) -TE ( 3 ) +0. 5*Y3 ( 4 ) /S ( 4 ) ) 
XB2=W*X1< 1)*X1(4)/(X1(3)*(X1(4)-X1(1) 1**2  1 
Y=-CPP*X1(1I*Y1(4)*6*XB2*XA/(D*X1(3) ! 
X2X1 (2 ) =X2X1 ( 2 )+Y-CCD*XB2*XSl/U 

X2X1(2)=X2X1(2>+(PC-CC)*PT*X1<4)/X1<3)+CC*Q*PT*X1<4)/(XH3)*BT) 

> 2X3(2)=X2X3(2)-CC*CT/BT 

i.l3  =  -X2Y3(3)/XlY3<3) 

Z12=Z13*X1  (4)/XK3)+X2XK2  1 

Z32=X2X3(2) 

Z11  =  Z12*XK3)/X1(2)+X2X1(1) 

Z31=X2X3(2) 

H1Y11=Z11*X1Y1(1)+X2Y1(1) 

H1Y31=Z11*X1Y3(1)+X2Y3(1)+Z31 

H2Y12=Z12*X1Y1(2)+X2Y1(2) 

H2Y32=Z12*X1Y3(2)+X2Y3(2)+Z32 

H3Y13=Z13*X1Y1(3J+X2Y1(3) 

IFIM-1 166,66,67 

66 

TE(1!=1.01+((X1(1)+Y1(2)*X 1(2)1/(1. +Y 1(2)1 +X 1(2)1/0. 06 
TE(2)=1. 0075+1 (X1(2)+Y1(3)*X 1(3)1/ (1. +Y 1(3) 1+XK3) 1/0.0694 

TE(3)=0.32  +  (  (X1(3)  +  Y1(4)*X1(4)  1/11.  +  YK4)  1+XK4)  1/0.0630 
M  =  M+1 

GC    TC    41 

67 

M  =  l 
IF(ABSF(H1YU)-ER)47,47,51 

47 

I  F(  ABSF  (  HI  Y3D-ER  148,48,51 

48 

I F(ABSF(H2Y12)-ER 149,49,51 

49 

I F( ABSF (H2Y32J-ER 150,50,51 

50 

IF(ABSF(M3Y13)-ER)68,68,51 

51 

IF    (H1Y11)     137,137,138 

137 

TV(1)=-1. 

138 

IF    (H2Y12)     139,139,140 

139 

TV(2)=-1. 

140 

IF    (H3Y131     141,141,142 

141 

TV(3)=-1. 

142 

IF    (H1Y31)     143,143,144 

143 

TV(4)=-1. 

144 

IF    (H2Y32)     145,145,146 

145 

TV(5)=-1. 

146 

IF(WC-X3(4)  136,36, -35 

36 

DC    37    1=1,5 

37 

DS( I )=u.5*DS( I ) 

Y1(2)=Y11T 

Y1(3)=Y12T 

Y1(4)=Y13T 

V3(2)=Y31T 

Y3(3)=Y32T 

GC    TC    38 

Table      A-3    (oontii 

DTX1(K)=XQ/(X1<  1-1  )*(X1  (  I-D+Yl  (  I  |*X1  (1  )  )  ) 

X=(-1./(X1(I-1 !**2)+l./(Xl ( 1-1 )*X1( 1 )) ) 

X=X/(DT( I )-TE( 1-1)+0.5*Y3( I )/S(I) ) 

YA  =  DTX1!K)*(  l./XK  I-D-l./Xl  (I  )  ) 

YA=YA/( (DT(I)-TE(I-l)+0.5*Y3<l)/S< I))**2) 

XA«EXPF(-HT/{R#X3(  I'-l)  )  J-EXPF <-HT/ (R*(X3<  I-D+Y3  ( I  )  )  > ) 

Y=-CPP*X1(1)*Y1( I )*B*XA*WA/(D*X1( I-l)**2) 

X2X1(K)=Y+CCD*HT*WA*X1 ( 1)*(X-YA)/U 

X5=EXPF(-HT/(R*X3( 1-1 ) ) ) / ( X3 ( I -1 ) **2 ) 

YA  =  EXPF(-HT/(R*(X3(  I -1  )+Y3  (  I  )  )  )  )  /  1  (  X3  U-l  )+Y3  (  I  )  )  **2  ) 

45 

X2X3(K)=1.0*CPP*X1( 1>*Y1(I )*WA*B*HT*(X5-YA)/(R*D*X1( 1-1) ) 
XB1=X1(1)*HT*(1.0/X1(3)-1.0/X1(4) ) / ( DJ ( 4 ) -TE ( 3 ) +0. 5*Y3 ( 4 ) /S ( 4 ) ) 
XB2=W*X1< 1)*X1(4)/(X1(3)*(X1(4)-X1(1) 1**2  1 
Y=-CPP*X1(1I*Y1(4)*6*XB2*XA/(D*X1(3) ! 
X2X1 (2 ) =X2X1 ( 2 )+Y-CCD*XB2*XSl/U 

X2X1(2)=X2X1(2>+(PC-CC)*PT*X1<4)/X1<3)+CC*Q*PT*X1<4)/(XH3)*BT) 

> 2X3(2)=X2X3(2)-CC*CT/BT 

i.l3  =  -X2Y3(3)/XlY3<3) 

Z12=Z13*X1  (4)/XK3)+X2XK2  1 

Z32=X2X3(2) 

Z11  =  Z12*XK3)/X1(2)+X2X1(1) 

Z31=X2X3(2) 

H1Y11=Z11*X1Y1(1)+X2Y1(1) 

H1Y31=Z11*X1Y3(1)+X2Y3(1)+Z31 

H2Y12=Z12*X1Y1(2)+X2Y1(2) 

H2Y32=Z12*X1Y3(2)+X2Y3(2)+Z32 

H3Y13=Z13*X1Y1(3J+X2Y1(3) 

IFIM-1 166,66,67 

66 

TE(1!=1.01+((X1(1)+Y1(2)*X 1(2)1/(1. +Y 1(2)1 +X 1(2)1/0. 06 
TE(2)=1. 0075+1 (X1(2)+Y1(3)*X 1(3)1/ (1. +Y 1(3) 1+XK3) 1/0.0694 

TE(3)=0.32  +  (  (X1(3)  +  Y1(4)*X1(4)  1/11.  +  YK4)  1+XK4)  1/0.0630 
M  =  M+1 

GC    TC    41 

67 

M  =  l 
IF(ABSF(H1YU)-ER)47,47,51 

47 

I  F(  ABSF  (  HI  Y3D-ER  148,48,51 

48 

I F(ABSF(H2Y12)-ER 149,49,51 

49 

I F( ABSF (H2Y32J-ER 150,50,51 

50 

IF(ABSF(M3Y13)-ER)68,68,51 

51 

IF    (H1Y11)     137,137,138 

137 

TV(1)=-1. 

138 

IF    (H2Y12)     139,139,140 

139 

TV(2)=-1. 

140 

IF    (H3Y131     141,141,142 

141 

TV(3)=-1. 

142 

IF    (H1Y31)     143,143,144 

143 

TV(4)=-1. 

144 

IF    (H2Y32)     145,145,146 

145 

TV(5)=-1. 

146 

IF(WC-X3(4)  136,36, -35 

36 

DC    37    1=1,5 

37 

DS( I )=u.5*DS( I ) 

Y1(2)=Y11T 

Y1(3)=Y12T 

Y1(4)=Y13T 

V3(2)=Y31T 

Y3(3)=Y32T 

GC    TC    38 

Table  A-3  (continued) 

165 

35  WC=X3(4) 

38  Y11T=Y1(2) 

Y12T=Y1(3) 

Y13T=Y1(4) 

Y31T=Y3(2) 

Y32T=Y3(3) 

Y1(2)=Y1(2)-TV(1)*DS(1) 

Y1(3]=Y1<  3>-TV(2)*DS(2) 

Yl  (  4)=YH4)-TV<3>*DSC3) 

Y3(2)=Y3(2)-TV(41*DS(4)                     '  " 

Y3(3)=Y3(3)-TV<5)*DS(5) 

DO  32  1=1,5 

32  TV( I )=1. 

. 

GO  TC  41 

68  IFITX-1. 169,69,63 

69  IFITN-1. 172.78.79 

72  IFtTZ-1. 1150, 151. 152 

150  QF(1)=Q 

XF(11=X1(4) 

WCI  (1)=X3<4) 

TZ  =  1. 

GO  TO  125 

151  0F(2)=O 

XF(2)=X1(4! 

WCI (2)=X3(4) 

TZ  =  2. 

GO  TO  125 

152  QF(3)=Q 

>  F(3)=X1(4) 

WCI  (31=X3(4) 

98  IFIWCI (l)-WCI (2) 173.73.74 

74  Q=QF(1 1 

X1(4)=XF( 1) 

X3(4)=WCI (1) 

OF( l)=OF(2) 

XF(1)=XF(2) 

WCI (1)=WCI (2) 

QF(2)=Q 

XF(2)=XK4) 

WCI (2)=X3(4) 

73  IFIWCI (21-WCI (3) 175,75.76 

76  Q=0F(3) 

XI (4)=XF( 3) 

X3(4)=WCI (3) 

QF(3)=GF(2) 

XF(3)=XF(2) 

WCI (3)=WCI (2) 

QF(2)=0 

XF(2)=X1(4) 

WCI (2)=X3(4) 

IFIWCI  (21-WCI  (  1')  177,77,75 

77  Q=GF(2) 

X1(4)=XF(2) 

X3(4)=WCI (2) 

0F(2)=0F(1) 

XF(2)=XF(  1) 

WCI  (2)=WCIU) 

0F(1)=Q 

Table  A-3  (continued) 

165 

35  WC=X3(4) 

38  Y11T=Y1(2) 

Y12T=Y1(3) 

Y13T=Y1(4) 

Y31T=Y3(2) 

Y32T=Y3(3) 

Y1(2)=Y1(2)-TV(1)*DS(1) 

Y1(3]=Y1<  3>-TV(2)*DS(2) 

Yl  (  4)=YH4)-TV<3>*DSC3) 

Y3(2)=Y3(2)-TV(41*DS(4)                     '  " 

Y3(3)=Y3(3)-TV<5)*DS(5) 

DO  32  1=1,5 

32  TV( I )=1. 

. 

GO  TC  41 

68  IFITX-1. 169,69,63 

69  IFITN-1. 172.78.79 

72  IFtTZ-1. 1150, 151. 152 

150  QF(1)=Q 

XF(11=X1(4) 

WCI  (1)=X3<4) 

TZ  =  1. 

GO  TO  125 

151  0F(2)=O 

XF(2)=X1(4! 

WCI (2)=X3(4) 

TZ  =  2. 

GO  TO  125 

152  QF(3)=Q 

>  F(3)=X1(4) 

WCI  (31=X3(4) 

98  IFIWCI (l)-WCI (2) 173.73.74 

74  Q=QF(1 1 

X1(4)=XF( 1) 

X3(4)=WCI (1) 

OF( l)=OF(2) 

XF(1)=XF(2) 

WCI (1)=WCI (2) 

QF(2)=Q 

XF(2)=XK4) 

WCI (2)=X3(4) 

73  IFIWCI (21-WCI (3) 175,75.76 

76  Q=0F(3) 

XI (4)=XF( 3) 

X3(4)=WCI (3) 

QF(3)=GF(2) 

XF(3)=XF(2) 

WCI (3)=WCI (2) 
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WCI  (2)=WCIU) 

0F(1)=Q 

Table  A- 3  (continued) 
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78 


79 
80 

81 


83 
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Table  A-3  (continued) 
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XF(3)=0.5*(XF(3)+XF(  1  )  ) 

Q  =  OFni 

X1(4)=XF(3) 

TN=5. 

PRINT  31.X1I4) ,Q,TN 

GO  TO  126 

94  WCI !3)=X3<4) 

QF<2)=ti.5*(QF(2)+QF<  1  )  ) 

XF(2)=u.5*(XF(2)+XF( 1  )  ) 

Q=QF(2) 

XI (4)=XF(2> 

TN  =  6. 

PRINT  31 .Xl(4) ,Q,TN 

TO  TO  126 

95  WCI (2)=X3(4) 

96  TEST=( ( (WCI t D-WCI (4) )**2.+<WCI (2)-WCI (4) )**2.+<WCI (3)' 

-WCI (4) )**2. 

1 ) /3. )»».S 

PUNCH  33.IXF1I ) ,QF( I ),WCI ( I ) ,1=1,3) 

IF(TEST-ERR0R)97,97,98 

97  TX=2. 

GO  TO  98 

131  Q=QF(1) 

X1(4)=XF( 1) 

PRINT  31.XK4)  ,Q,TX 

GO  TO  126 

63  PUNCH  5 

PUNCH  25,X1<4) ,0,X2(4) 

PUNCH  6.YK2)  ,Y1(3)  ,Y1(4) 

PUNCH  7.Y3I2) ,Y3(3) ,Y3(4) 

PUNCH  8.XK1)  ,X1(2)  ,X1(3) 

PUNCH  9,DT(2) ,DT(3) ,DT(4) 

DO  64  1=2,4 

64  X3( I )=X3(  I  1-460. 

PUNCH  10,X2(1) ,X2(2) ,X2(3) 

PUNCH  11.X3I 1) ,X3(2)  ,X3(3) 

PUNCH  12,H1Y11,H2Y12,H3Y13 

PUNCH  13,H1Y31,H2Y32,WPR(4> 

PUNCH  14 

PUNCH  15, CI, PCI 

PUNCH  16.C2.PC2 

PUNCH  17,AR( 1) ,PC3 

PUNCH  18,PU( 1) ,PC4 

PUNCH  19.C5.PC5 

PUNCH  20.C6.PC6 

PUNCH  21 

PUNCH  22,T01,TD2,TD3 

PUNCH  23,WPR( 1 ) ,WPR(2) ,WPR(3) 

PUNCH  27.ARI2) .AR(3) .AR(41 

PUNCH  28, PU(2),PU(3) ,PU(4) 

PUNCH  29,TE(1) ,TE(2) ,TE(3) 

PUNCH  24 

GO  TO  71 

TND 
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?a\)la  A-4       Input   <Jats  and  SanipXe   Output  ?. 

ooult:: 

Ti-!5  following  are   input  data 

xi  ui                o                   wc               recii 

.06  ,005-00.2;-      OOOS-rGl    1. 000005-00    2.303505-00 

T5;2;                T5:5,   ■ 

2.393605-00    2.231605-00 

*  *       '                    i  *  k  j  j                    t  1 1 4 )                    1 3  (  2  ) 

Y '-  {  '  ' 

.24029E*C:    0.327125-rOl    0 .44967-£*0X-0.543*37£;^O2.. 

0.430975^02 

^S<l>_           _    0S(2)                   35(3)                   DS{4j 
j.^0000c-01    5.000005-Cl    5.000005-01    i.0000C5~C0 

3s:5;                   ER 

.•u-;^&ui"oG    -.C0CG0-~°G4 

U                A              riT              ,-iTS           X4  ( X  i            0             -S  (  2  J 

3  '  3  j          S  f  4 ''                V 

510.0       734.4    1000.0      923.9      710.0         62.5   '      23.0 

23.0         22.0            1.0 

3                   a1(1!                 CST    '                  W                      CH7 
l»X04E-CX    0.0355-rCG    0.0255-02    3.3405tG3    0.3765-0 

4    2.3975-05    2.9035-09. 

5                           C5                  .         PC                           C? 

cc 

;.7900C£T09    2.0S0005-02    1.796005-06    1.0C00C5-00 

5.937505-07 

~"F    •rC_-Ow:\G    AR5  .OPTIMUM   Cu'TPoT    DATA 

.1;-./-    6.5v5G.C-£-C2                             G=    2  .700005-1-01 

X3:4;=    2.5337^5-01 

Yli2;-    2.13o89i-0u                   Yli3)=    2.876735-00 

Yi£4.J=    3.350915-00 

v~  ;2:>  =-3.4-402 5-rOl                   Y3(3}=-4.e0402E+01 

Y3!4!=-4. 4-5277  £tC1 

XI;  11  =    3.'5CO0C!5-02                   XI  (2)=    4.192245-02 

X )  v  3  )  =    5.123/35-02 

O.T:.?)=    3.6-1765-00                     Oj;3i=    3.055995-00 
XXC'-    1.5336:5-01                   X0;2;=    1.371205-01 

ST<4)=    7.526^95-00 
X£.:3;=    2.177e4£-01 

Xi;?i«    :.05950c+'w2                  *3;3!=    X.47919E+02 
hi?'.!"    7.2^.',005-'07                    :  i  2  Y 1 2  =    2.304405-05 

X5'4;=    1.033915-02 

"..•>'.-    6.393435-05                   H2Y32=    1.702305-05    ■ 

-  '    ■'-                                                        CCSTiSi                     ?5RC5.\TAG5 

■    -" '•                                             1.3-3035-01                  4.597995-rOl 

-rj-~-~<                                         1.253265-0?                  4.333335-01 

::'.".r-:i:.No  ar5a            s.47040E-02            2.966055-rOi 

"'    ''?:n'CS                               .     9.359365-03                 3.277345-00 

"■.:--'JCT:0.\                          2.030005-02                  7.233435-00 

s"-  =  T?<eATy£NT                        3.313255-02                 -.335975-rOl 

"-    v~                              1    ST    EFFECT               2    NO    EFFECT 

3    xZ    5  F  F  5  C  T 

'■'■■.    OXC-                     3.404025-rOl                  4.3-4023-01 

-r.45^7*7"  —  *"  ' 

•'■"—■"    PntOO).                 2.9o3g05t03                 2.  74i22£i-03 

2.612v7£i-03 

CCS.    COS"                    2.773965-02                  2.323115-02 

2.871335-02 

P-U'^?.    COST                    6.019495-03                  2.412395-03 

n     -  _  ^  .  ~  r  _  - 

i-p-     =1-.                          2.371445-00                  2.447^75-00 

• 
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In  PART  ONE  a  detailed  analysis  of  a  HEMS  prooos3  is  made  and 
a  mathematical  model  of  the  process  is  developed.   An  optimization 
study  of  such  a  model  is  carried  out  by  a  discrete  analog  of  the 
maximum  principle  in  conjunction  with  two  search  techniques:  the 
parametric  search  and  the  simplex  method.   Both  methods  lead  to 
the  same  optimal  results.   In  contrast  to  the  parametric  search, 
the  simplex  method  gives  rise  directly  to  the  optimum  point.  The 
parametric  search,  however,  gives  detailed  information  about  the 
influences  of  the  individual  parameters  on  the  water  cost  and  the 
other  operating  variables.   In  PART  TWO  a  general  mathematical 
model  of  a  sequential  multistage  reverse  osmosis  process  is  deve- 
loped. This  model  is  obtained  under  the  assumption  of  plug  flow 
inside  the  tubular  osmosis  unit  to  take  into  account  the  brine 
concentration  changes  along  the  membrane  tube.   Several  simplified 
versions  of  this  model  are  also  proposed. 
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influences  of  the  individual  parameters  on  the  water  cost  and  the 
other  operating  variables.   In  PART  TWO  a  general  mathematical 
model  of  a  sequential  multistage  reverse  osmosis  process  is  deve- 
loped. This  model  is  obtained  under  the  assumption  of  plug  flow 
inside  the  tubular  osmosis  unit  to  take  into  account  the  brine 
concentration  changes  along  the  membrane  tube.   Several  simplified 
versions  of  this  model  are  also  proposed. 


